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FOREWORD 


The  proceedings  of  the  Third  Semi-Annual  ANP  Radiation 
Effects  Symposium,  held  at  the  DInkler-Plaza  Hotel  In 
Atlanta,  Georgia,  October  28  through  30,  1958,  are  In 
six  volumes.  Each  of  the  first  five  volumes  presents  the 
unclassified  papers  from  one  of  the  five  sessions;  the  sixth 
volume  presents  classified  papers  from  all  five  sessions. 

Each  volume  contains  a  complete  table  of  contents  and 
an  Index  of  authors.  Volume  One  contains  a  list  of  tbe 
names  of  all  who  attended  the  Symposium. 
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PREcDICTING  THE  PERFORMANCE  OF  IRRADIATED  ELECTRONICS 
SYSTEMS  BY  SIMULATION  ON  THE  ANALOG  COMPUTER 

by 

V.  C.  Brown  and  N.  M.  Peterson 
CONVAIR 

A  Division  of  General  Dynamics  Corporation 
Fort  Worth,  Texas 


The  method  of  simulating  electronics  systems 
on  analog  computers  for  the  purpose  of  studying  the 
effects  of  radiation  is  described  and  a  comparison 
made  between  the  results  obtained  by  the  simulation 
technique  and  data  obtained  in  an  actual  irradiation. 
The  results  show  that,  by  combining  proper  statisti¬ 
cal  methods  with  computer  simulation,  the  performance 
of  electronics  systems  under  radiation  may  be  suc¬ 
cessfully  predicted. 


INTRODUCTION 

Literally  thousands  of  components  are  available  today  to  the 
military  electronics  designer;  unfortunately,  only  a  few  will  ever 
be  adequately  tested  in  the  many  possible  combinations  of  elec¬ 
tronic  application,  temperature,  humidity,  and  pressure.  For 
this  reason,  the  functioning  of  a  mechanism  as  complex  as  many 
of  the  present-day  electronic  systems  is  always  a  question. 

When  the  system  is  exposed,  in  addition,  to  such  a  new  and  here¬ 
tofore  unstudied  environment  as  nuclear  radiation,  the  question 
is  compounded. 

Predictions  of  the  effects  of  radiation  on  electronic 
systems  must  take  into  direct  account  the  I'adiation-induced 
changes  in  the  operating  characteristics  of  the  components.  The 
tolerance  of  the  system  to  these  changes  and  the  possibility  of 
feedback  or  cumulative  effects  -  that  is,  effects  not  directly 
associated  with  the  radiation  damage  itself,  but  changes  caused 
by  a  damage  to  one  component  affecting  the  operation  of  another  - 
must  be  considered. 

One  approach  to  this  problem  has  been  through  the  use  of 
simulation  techniques  on  the  analog  computer.  The  method  con¬ 
sists  essentially  of  simulating  the  system  to  be  studied  and 
then  varying  the  characteristics  of  the  components  as  they  are 
found  to  vary  under  irradiation. 
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SIMULATION  TECHNIQUES 


Analog  computers  may  be  used  to  simulate  physical  systems 
in  several  ways,  depending  on  the  nature  of  the  problem.  In 
general,  as  shown  in  Figure  1,  the  principal  means  of  simula¬ 
tion  are  the  following: 

(1)  Solution  of  the  equations  governing  the  system; 

(2)  Direct  substitution  for  system  components,  i.e., 
system  components  are  replaced  with  analog  computer 
elements  directly, 

(3)  Simulation  of  system  output,  in  which  the  output 
from  the  computer  simulates  that  of  the  physical 
system.  The  method  of  obtaining  the  simulated 
output,  however,  bears  no  direct  resemblance  to 
the  physical  system. 

The  solution  of  system  equations  is  the  most  widely  used  method 
however,  it  usually  requires  the  largest  number  of  computing 
elements  and,  therefore,  is  the  most  expensive  method  of  simu¬ 
lation. 

On  the  other  hand,  the  method  of  direct  substitution  for 
the  system  components  is  simpler  and  easier  to  apply.  It  has 
the  added  advantage  that  the  equations  of  the  system  need  not 
be  written.  The  simulation  of  system  output  is  generally  used 
when  either  of  the  other  methods  is  impractical  or  impossible 
to  apply.  The  advantage  of  this  method  is  that  very  complex 
equipment  need  not  be  simulated  in  its  entirety. 

Of  the  three  methods  generally  used,  experience  has  in¬ 
dicated  that  the  direct  substitution  method  is  the  one  most 
suited  for  electronic  systems.  Obviously,  it  is  simpler  to 
simulate  one  LC  network  with  another  LC  network  than  it  would 
be  to  write  the  equations  and  then  study  the  behavior  of  the 
equations.  In  effect,  then,  in  the  case  of  passive  elements, 
the  simulation  consists  of  building  the  system  under  study  out 
of  the  precision  variable  components  of  the  computer. 

Operation  of  the  vacuum  tubes  and  other  active  elements 
in  an  electrical  circuit  may  be  accurately  simulated,  as  in 
Figure  2,  with  current  generators  or  voltage  amplifiers  con¬ 
nected  to  give  the  familiar  equivalent  circuits.  To  take  into 
account  the  nonlinearities  in  the  transfer  characteristics  of 
such  elements,  a  function  generator  may  be  used  in  conjunction 
with  the  other  elements. 
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METHODS  OF  SiMUL^mO^ 
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TIME — 1> 

SYSTEM  RESPONSE 
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NONUNEAR  SIMULATION 


Figure  3  llluBtratea  hew  a  biasing  arrangement  is  used  to 
simulate  the  cutoff  point  of  a  tube.  Notice  the  similarity  be¬ 
tween  the  actual  and  the  sxmulated  systems.  The  "feel"  for  the 
system  is  not  lost,  and  each  component  in  the  original  system 
has  a  counterpart  in  the  simulated  system.  The  value  of  any 
of  the  components,  or  the  characteristics  of  the  tubes,  may  be 
varied  as  needed. 


APPLICATION  TO  SYSTEMS 

The  reliability  of  the  predictions  is  only  as  good  as  the 
data  that  are  used.  In  this  respect,  the  analysis  of  a  system 
begins  with  an  analysis  of  the  available  irradiation  data  on 
the  components  used  in  the  system. 

Suppose,  for  pui’poses  of  i  ilusti  aticn,  tliat  a  group  of 
resistors  of  a  certain  type  have  been  irradiated.  Then  the 
mean  resistance  of  the  group  may  be  plctred  as  a  function  of 
irradiation.  Next,  by  standard  st'at,.otical  'rethods.  a  corj- 
fidence  envelope  may  be  drarvn  around  the  mean  value  of  resis¬ 
tance,  within  which  a  given  percenr  of  all  the  resistors  of 
that  class  would  lie,  as  in  Figure  L.  Tl»e  width  of  the  con¬ 
fidence  band  depends  upon  t.he  per  sen',  o;-'  the.  total  number-  of 
components  that  we  want  to  fail  within  the  band  and  upon  the 
confidence  desired  in  the  predictions. 

Now  suppose  that  we  .have  .olrralareJ  a  system  on  a  computer, 
the  output  of  which  is  a  function  of  the  resistors  within  it. 

By  varying  the  values  of  •:hese  resist  .arc,  there  may  be  either 
an  increase  or  a  decrease  in  the  uy.iality  of  the  system  perform¬ 
ance,  In  other  words,  the  partial  de’-ivative  of  the  system 
output  with  respect  to  the  i-esistance  is  either  positive  or 
negative.  If  it  is  positive,  the  -.'aiaes  of  the  resistance 
used  in  the  simulated  irradiated  system  wil 1  be  those  taken 
from  the  bottom  of  the  confidence  en’-elope.  If  tne  partial 
derivative  i...  negative,  the  top  of  the  ccnflderce  band  is 
used.  In  this  way,  by  taking  the  values  at  che  extreme  of 
the  envelope,  we  have  ensured  with,  say,  h'd  9h  percent 
confidence  that  the  resistor  iri  tne  actual  systeui  will  be  in 
better  condli-ion  than  in  the  simuiabed  system.  Incidentally, 
in  the  determination  of  the  pax’tial  deviatives  one  finds  im¬ 
mediately  the  most  critical  co.npcnec.ts  where  the  system  will 
require  close  parametric  control. 

Figure  5  gives  tne  results  of  an  aralys.Ls  of  a  video 
transistor  amplifier.  The  actual  amplifier  was  irradiated  in 
the  Penn  State  reactor  and  was  a  two -stage,  capacitance-coupled 
.amplifier  with  voltage  feedback  and  unbypassed  emitter  resis¬ 
tors.  The  deviation  between  the  curves  given  for  the  simulated 
and  actual  system  is  believed  to  be  caused  partly  by  the  lack 
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of  uniformity  normally  found  in  tranolotorc; .  In  addition, 
some  healing  of  the  amplifier,  when  the  irradiation  of  the 
real  system  was  interrupted  ovornight,  was  apparent.  The 
healing  effect  was  not  Included  in  the  dats  used  in  simula¬ 
ting  the  system  on  the  computer, 

OVER-ALL  PERFORMANCE 

Finally,  it  must  be  recognized  at  some  time  the  system 
may  refuse  to  perform  at  all.  That  Is,  there  may  be  a  "cata¬ 
strophic"  component  failure,  such  as  a  capacitor  short  which 
renders  the  system  completely  inoperative. 

Ordinarily,  standard  statistical  methods  are  used  to 
evaluate  a  system  subject  to  catastrophic  failures  and  the 
reliability  is  given  as  an  exponential  function  of  time  (Fig. 
6),  In  Irradiated  systems  the  picture  is  altered  somewhat, 
the  probability  of  failure  in  a  system  being  not  only  a 
function  of  time  but  of  the  radiation  administered  to  the 
system. 

If  it  is  assumed  that  the  catastrophic  failure  rate  of 
the  components  in  a  cystem  can  be  approximated  by 


f(t)  r  kf  , 


then,  when  the  system  is  oi^ei’ated  continuously  under  radia 
tloii,  the  reliability  may  be  given  by 

P(t)  = 


where  P(t)  represents  the  probabix'ty  of  the  system  operating 
for  a  time  t  without  failing,  and  and  k  are  parameters  which 
are  functions  of  the  ir .'rf^tad  doi^e  administered  to  the  sys¬ 
tem.  The  development  of  this  expression  is  discussed  in  the 
Appendix.  Portions  of  the  failure  rate  curve  may  be  approx¬ 
imated  by  this  function,  and  the  z*eliability  may  be  calculated 
if  the  failure  rates  for  the  various  components  are  known. 

Figure  7  gives  uhe  results  of  an  analysis  on  a  closed 
circuit  television  system,  wnose  can. era  is  subjected  uo  radia¬ 
tion,  Curve  2  shows  the  magnitude  of  the  optical  signal  into 
the  camera,  and  Curve  3  is  an  approximation  of  the  reliability 
of  the  system.  The  electronic  perfor'mance  (the  gain -bandwidth 
product)  was  obtained  by  the  computer  simulation  technique 
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and  the  optical  performance  from  data  on  irradiated  optical 
glasses . 

Such  an  analysis  of  a  system  is  based  on  answering  three 
questions , 

1.  As  the  system  is  irradiated^  how  will  it  operate? 

2.  What  is  the  chance  of  its  operating  at  all? 

3.  When  is  the  system  expected  to  fail? 

If  the  failure  of  the  system  is  caused  by  the  catastro¬ 
phic  failure  of  components,  ther3  its  mean  time  to  failure  is 
simply  the  mean  time  to  failure  of  the  most  susceptible  com¬ 
ponents.  If  the  performance  of  the  system  is  considered  un¬ 
satisfactory  at  some  arbitrary  level ,  then  its  mean  time  to 
failure  may  be  found  approximately  by  following  the  mean- 
value  curves  of  the  components  under  irradiation. 

The  reliability  curve  in  Figure  7  is  somewhat  arbitrary 
in  that  it  was  computed  from  the  probability  that  the  system 
is  to  operate  for  a  specified  length  of  time  without  a  failure. 
By  taking  shorter  or  longer  periods  of  time,  the  curve  may 
be  raised  or  lowered. 

Of  the  three  such  systems  irradiated  to  date,  the  re¬ 
sponse  seems  to  follow  the  optical  performance  curve  given 
in  the  figure.  There  are,  however,  other  considerations  to 
be  taken  into  account,  such  as  lighting  and  optics  adjust¬ 
ments,  that  were  not  included  In  the  analysis.  These  could 
shift  the  opti'a'l  curve  either  one  way  or  another.  Of 
the  three  cameras  irradiated  to  a  total  integrated  neutron 
flux  of  about  8  x  10^^  nvt,  no  change  in  electronic  perform¬ 
ance  was  discernible.  There  were,  however,  several  failures 
in  the  camera  control  unit. 


CONCLUSIONS 


It  is  believed  that,  by  combining  proper  statistical 
methods  with  computer  simulation  of  electronic  systems, 
electronic  performance  under  irradiation  may  be  successfully 
predicted.  The  accuracy  of  these  predictions  depends  of  course 
upon  the  quality  and  quantity  of  the  available  data  on  the 
components  and  upon  the  application  of  these  data  to  a  particu¬ 
lar  situation.  When  only  trends  can  be  found  in  the  component 
data,  the  simulation  can  still  give  valuable  insight  into  system 
stability  and  performance  in  a  nuclear  environment. 
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APPENDIX  A 


FAILURE  DATA  AND  RELIABILITY 

A  failure  is  here  defined  as  any  deviation  of  the  function 
Ing  of  a  part  or  subsystem  from  requirements  of  the  specifica¬ 
tions,  drawings,  purchase  descriptions,  contracts,  or  order.  A 
critical  failure  is  one  which  could  prevent  the  performance  of 
the  tactical  function  of  the  system. 

Generally,  it  is  assumed  that  the  probability  of  failure 
is  small,  and  that  failures  are  distributed  according  to 
Poiss6n's  formula 


p  (x)  = 

It 


where:  x  -  the  number  of  failures  and 
T  -  mean  time  to  failure 

Any  test  or  run  of  the  system  can  have  the  possible  out¬ 
come  of  "failure"  or  "no  failure".  Each  run  is  then  a  Poiss6n 
Trial,  and  the  probability  of  no  failures  in  a  time  t  of  the 
trial  or  run  is 


(1) 


a  continuous,  exponential  distribution. 

Hence  the  probability  of  one  or  more  failures  in  time  t  is: 

Q(0)  =  1  -  P(0)  =  1  - 

and  the  probability  of  failure  in  time  dt  is  the  rate  of 
occurrence  of  failures: 

1  .  t/T 

dQ  =  f  e  dt. 


Another  way  of  saying  this  is  that  the  instantaneous  proba¬ 
bility  rate  of  failure  at  time  t  conditional  on  no  failure 
before  time  t  is: 

Z(t)  =  PQ  t<  r  <  tf  dt/T>t3|=  f(t) 

1  -  P(t) 
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where : 


If 


then 


f(t)  >  probability  of  failure  In  time  t. 


Z(t)  =  k,  a  constant 


Z(t)dt  = 


kdt  ■  kt 


and 


kt 


o 


or 


hence 


and 


1  -  F(t)  r  e"^"^  =  reliability  of  the  system; 


P(t)  =  1  -  e”^^ 


dF(t)  =  f(t)  =  ke"^^ 


or 


f(t)  =  1 
T 


t 

T 
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Thus 


r 

2 

J  f(t)dt  =  1-e  r  — - 

0  T  ' 

'  t  (  1 

, "t  )  Ti 

I  ^ 

i  T 

«—  for 
T 


large 


T 


If  k  =1/Tis  not  constant  with  respect  to  time,  then  Figure  6A 
would  show  a  slope  0.  Suppose l/T=  kt^,  then  the  distribution 

of  failures  instead  of  following  the  simple  exponential  distribu¬ 
tion,  would  adhere  to  the  Welbull  distribution 

Z(t)  r  kt“ 

,.a 

l^n|l-F(T)j  s-kt**  and  1  -  F(t)  z  e  , 

which  may  be  taken  as  defining  the  reliability  of  the  system 
under  study. 
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COMBINED  TIME,  TEMPEli/iTUltE,  AND  RADIATION  EFFECTS 
ON  CBGANIC  MATERIALS 


by 


Co  Go  Collins 
General  Electric  Company 
Aircraft  Nuclear  Propulsion  Department 
Cincinnati  15,  Ohio 


Experimental  observations  of  combined  time,  temperature, 
and  radiation  effects  on  a  lubricating  oil  and  Teflon  aircraft 
engine  hoses  are  described  for  the  temperature  range  from  80° 
to  /j.00°F  and  for  radiation  rates  of  106  to  5  x  lO^  ergs  (gmC)-l 
hr-1. 


The  results  show  that  at  a  constant  radiation  rate  the 
failure  time  -  as  judged  by  the  oxidation-induction  period  of 
the  oil  and  by  leaks  in  the  Teflon  hoses  -  follows  an  Arrhe¬ 
nius-type  relationship  with  temperature.  At  a  constant  tem¬ 
perature,  failure  time  was  found  to  follow  a  logarithmic  re¬ 
lationship  with  radiation  dose  rate.  The  over-all  results 
can  be  described  in  a  summary  equation  in  terms  of  time,  ab¬ 
solute  temperature,  and  radiation  dose  rate. 

Several  implications  of  the  summary  equation  are  dis¬ 
cussed  in  reference  to  radiation  effects  in  organic  materials. 
One  implication  is  that  observations  of  equal  damage  at  con¬ 
stant  dosage  irrespective  of  dose  rate  may  be  valid  at  only 
one  temperature  or  at  high  dose  rates. 


INTRODUCTION 


One  of  the  major  problem  areas  confronting  radiation  effects  inves¬ 
tigators  at  present  is  that  of  defining  materials  performance  under  com¬ 
bined  radiation  and  temperature  environments.  Since  time-at-temperature 
is  an  important  factor  in  the  deterioration  of  organic  materials,  the 
importance  of  temperature  xn  radiation  effects  has  been  suspected  -  and 
demonstrated  -  in  a  number  of  instances.  Progress  in  this  area,  however, 
has  been  impeded  by  the  lack  of  mathematical  methods  of  treating  the 
data,  and  large  numbers  of  experiments  have  been  necessary  to  obtain  data 
for  each  parameter  over  a  range  of  values  of  practical  interest. 
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This  report  discusses  experiments  in  which  the  time  required  for  oc¬ 
curence  of  a  constant  amount  of  oxidation  of  an  organic  liquid  was  studied 
over  a  range  of  dose  rates  and  temperatures.  The  quantitj’'  -  time  required 
for  constant  amount  of  reaction  -  is  readily  analyzable  with  the  familiar 
logarithm  time  versus  absolute  temperature  relationship  of  accelerated 
life  tests.  It  proved  fairly  straight-forward  to  analyze  the  radiation 
experiments  in  a  similar  manner  and  the  results  of  the  experiments  can  be 
described  by  an  equation  involving  terms  for  time,  temperature,  and  ra¬ 
diation  rate.  Further  experimentation  is  necessary  to  clearly  establish 
the  validity  of  the  eouation;  however,  the  general  form  of  the  equation 
seems  clear  and  bears  sufficient  resemblance  to  fundamental  relationships 
that  it  may  be  applicable  to  other  materials.  In  this  respect,  the  im¬ 
plications  were  considered  of  sufficient  interest  to  warrant  discussion 
despite  limited  experimental  verification.  Some  of  the  implications  of 
the  eouation  regarding  the  relationships  between  radiation  rate,  radia¬ 
tion  dosage,  temperature,  and  time  in  the  degradation  of  organic  materials 
are  discussed  following  a  description  of  the  experimental  procedure  and 
results.  A  few  data  that  indicate  that  the  same  type  of  equation  applies 
to  the  life  of  Teflon  hoses  are  included. 


EXPERIMENT.AL  PROCEDURE 


Experiments  were  performed  cn  a  di-ester  oil  and  utilized  the  oxida¬ 
tion-induction  period  as  the  indicator  of  the  occurrence  of  a  given  amount 
of  chemical  reaction.  The  oxidation -induction  period  is  defined  as  the 
time  required  to  reach  the  "break"  in  the  curve  of  oxygen  absorption 
versus  time,  as  shown  by  curve  a  of  Figure  i.  At  high  temperatures,  the 
oxygen  absorption  curve  rises  without  a  definite  break,  as  in  curve  b  of 
Figure  1,  so  that  the  induction  period  is  usually  defined  as  the  time  to 
absorb  the  same  amount  of  oxygen  that  is  absorbed  in  a  true  induction 
period  (curve  a).  In  tbese  experiment.-,,  the  effect  of  increasing  radia¬ 
tion  rate  was  observed  to  affect  the  curves  in  a  manner  similar  to  that 
of  higii  temperature  and  the  usual  definition  of  the  induction  period  was 
followed.  This  corresponded  to  the  absorption  of  O.h  mole  of  oxygen  per 
426  grams  of  fluid. 

The  apparatus  employed  was  of  the  type  described  by  Dcrnte.  This 
consisted  of  circulating  oxygen  through  the  oil  in  a  closed  system  con¬ 
taining  a  bellows  pump,  adsorbent  tubes  for  oxidation  products  (water. 
C02)v  and  a  mercury  valve  for  admitting  oxygen.  The  oxygar.  flow  rate, 
intended  to  be  ^.0  liters  per  hoar,  was  found  to  be  3.5  liters  per  hour 
after  a  few  measurements  were  completed.  All  of  the  experiments  reported 
here  were  conducted  at  this  flow  rate. 

Samples  consisted  of  250  millt liters  of  di -2-ethylhexyl ■ sebacate 
containing  0.5  weight  percent  phenothiazine.  This  quantity  of  oil  filled 
the  2  inch  diameter  stainless  steel  irradiation  capsules  to  a  height  of 
approximately  5  inches.  Oxygen  bubbled  through  about  As  inches  of  this 
height  from  perforations  in  a  horizontal  portion  of  the  gas  inlet  tube. 
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Dosimetry  for  the  irradiations,  which  were  performed  in  a  cobalt-60 
source,  was  determined  by  means  of  the  cerous-ceric  oxidation. 

Measurements  were  made  inside  the  furnace-capsule  arran.f^e.nent  used  with 
the  oil  and  corrections  applied  for  the  thickness  of  the  walls  of  the 
glass  cerous  sulfate  container.  Energy  absorption  in  the  samples  was 
computed  by  graphically  integrating  the  curve  of  dose  rate  versus  dis¬ 
tance  (from  the  source)  for  the  sample  geometry.  An  average  energy  of 
1.25  raev  was  assumed  for  the  gamma  rays,  and  the  absorption  coefficients 
were  those  reported  by  Snyder  and  Powell. (l) 
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TYi'lG.iL  OXYGiils  .ibouilPTlOh  CUhV.Gd  li.  OilD  -.TlOh 
IKDlOTlb!  PFhIOD 


Curve  (,a)  ;  at  low  temperature  or  lo'w  radiation  dose  rate 
Curve  (b)  ,  at  high  temper iture  or  hig!.  radiation  dose  rate 
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RESULTS 


The  results  of  the  measurements  are  summarized  in  Table  I.  Since 
the  oxidation  reaction  is  significantly  affected  by  temperature,  it  is 
not  surprising  that  the  results  are  not  simply  correlated  in  terms  of 
radiation  dosage. 

Figure  2  is  a  plot  of  the  results  in  the  type  of  graph  employed  in 
accelerated  life  tests.  The  logarithm  of  tlie  time  reauired  for  the  con¬ 
stant  amount  of  reaction  indicated  by  the  induction-period  is  plotted 
versus  the  reciprocal  of  absolute  temperature  (^Kelvin),  results  at  any 
one  radiation  rate  being  connected  by  a  straight  line.  That  the  lines 
are  straight  is  demonstrated  for  only  a  limited  temperature  range  by 
these  data  but  the  idea  is  supported  by  the  results  of  K,  A.  Golub  on 
the  cis-trans  isomerization  of  poly-bucadiene, (2)  by  the  Teflon  data  of 
Figure  6,  and  follows  from  the  of  chemical  kinetics. 

Under  non-radiation  conditions,  iho  time,  t,  reouired  for  a  given 
constant  amount  of  reaction  is  related  to  temperature  by  the  equation 

In  t  ~  +  constant  (l) 

Ki 

wtiere  is  an  experimental  activati-.n  energy,  R  is  the  molar  gas  con¬ 
stant,  T  is  absolute  temperature,  and  the  constant  .’.nvolves,  in  part  at 
least,  the  ratio  of  the  concentration  of  a  reacting  species  initially  and 
after  time  t.  This  eouation  is  derived  from  Arrhenius  relationship  be¬ 
tween  reaction  velocity  and  temperature,  and  is  the  familiar  basis  of  all 
accelerated  life  tests  of  materials  and  components. 

An  eouation  of  this  same  form 

In  t  =.  I  +  b  (II) 

is  applicable  to  the  experimental  points  obcained  at  constant  radiation 
dose  rate.  Considering  that  a  giver,  dose  rate  produces  the  reactant 
species  at  a  constant  rate,  the  >'ea-'T,on  valccity  vn.ll  he  accelerated 
over  that  orxurring  with  temperat.ire  alone.  The  reaction  velocity  enters 
(l)  in  the  terms  for  activation  en8r5,y  and  in  the  cc.istant,  hsnce,  one 
would  expect  the  basic  form  of  1  tc  be  the  same,  ana  the  effect  of  ra¬ 
diation  rate  in  the  empirical  relat j.onships  tc  be  apparent  in  the  values 
of  and  tne  constant,  i,e„  in  a  and  b.  On  this  basis,  also,  it  would 
be  expected  tnat  since  the  reaction  rate  is  accelerated  more  by  a  higher 
radiation  rate,  an  increase  in  dose  rate  should  decrease  the  numerical 
value  of  a. 

In  for.mulating  an  over-all  eouation  for  the  experimental  results, 
the  relationship  between  dose  rate  and  the  acceleration  of  reaction 
velocity  is  necessary;  this  proved  to  be  fairly  simple.  Straight  lines 
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TABLE  I 


OXIDATION- 

INDUCTION  PERIOD 
DURING 

OF  AN  INHIBITED  DI- 
lURADIATION 

ESTER  OIL-*> 

OXIDATION  INDUCTION 
PERIOD 
(hrs. ) 

TEI'IPERA'fURE 

"'F 

RADIATION  DOSE  RATE 
(absorbed) 
ergs  gm"!  hr~i 

RADIATION  DOSAGE 
(absorbed) 
ergs  gm”! 

77 

300 

2.2  X 

1.69  X  109 

49 

M 

5.3  X  10'^ 

2.6  X  109 

106 

1.1  X  lo'^ 

1.17  X  10® 

79 

tl 

2.-7  X  10^ 

2.13  X  10® 

36 

tl 

2.2  X  lo'^ 

7.92  X  10® 

23.5 

It 

5.6  X  lO"^ 

1.32  X  109 

31.5 

400 

1.0  X  10^ 

3.15  X  lo'^' 

16 

tt 

2.1  X  lO'^ 

3.36  X  10® 

*Di-2-ethylhexyl-sebacate  plus  0.5  per  cent  (weight)  phenoghiazine.  Oxygen 
flow  rate  vbs  3.5  liters  per  hour  through  250  cc. 
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\irere  obtained  when  the  constants  a  and  b  of  enuation  II  were  plotted  versus 
thelogarithm  of  radiation  rate;  the  constants  are  therefore  of  the  form: 


a  =  In  r  +  C2 


b 


C3 


In  r  + 


C3 


vriiere  r  is  dose  rate  and  the  "c's"  are  empirical  constants^  Substituting 
these  values  for  a  and  b  in  eouation  II  •'/a.elds  the  summary  eouation: 


In  t  =  ..  ( 


1,95434)  In  r  -  4?. 40723 


III 


where  t  is  hourSj,  r  is  radiation  rate  {ergs  per  gram  absorbed  per  hour), 
and  T  is  absolute  temperature  (oK).  This  equation  fits  the  experimental 
data  shown  in  Figure  2  within  ± 


ul 


.  joU  S3I0K 


'.i/hile  a  more  fundamental  approach  to  reaction  rates  might  have  yielded 
a  somevjhat  less  complicated-appearing  relationship,  the  empirical  equation 
derived  here  does  not  conflict  vdth  the  basic  ideas  of  photochemistry  and 
bears  sufficient  resemblance  to  I’undamental  relationships  that  an  eouation 
of  the  same  general  tyoe  possibly  will  apply  to  many  organic  materials 
under  radiation  conditions.  In  this  respect,  it  is  of  interest  to  consider 
the  ’’view"  o:  radiation  damage  presented  by  the  eouation  and  its  implica¬ 
tions  regarding  the  relationships  between  time,  temperature,  dose  rate, 
and  dosage. 

Over-all  Indications  -  The  results  obtained  from  calculations  with 
the  eouation  provide  the  over-all  "picture”  of  radiation  effects  shown  in 
Figure  3«  The  heavy  lines  on  the  right  and  at  the  bottom  of  the  figure 
represent  the  apparent  limits  of  applicability  of  the  enuation  and  are  dis¬ 
cussed  below,  'iVithin  the  region  bounded  by  the  limits,  the  combined  tem¬ 
perature  and  radiation  effects  are  indicated  by  the  lines  shown  for  a  few 
selectea  radiation  rates.  At  constant  dose  rate,  straight  lines  are  ob¬ 
tained  when  In  t  is  plotted  versus  II,  The  slope  of  these  lines  is 
greatest  at  low  dose  rates,  thennni  activation  being  the  major  factor  in¬ 
fluencing  the  reaction.  At  successively  higher  dose  rates,  the  lines 
tend  to  a  slope  of  zero,  the  indication  being  that  radiation  effects  com¬ 
pletely  over-ride  thermal  effects. 

The  limits  of  the  region  over  wnich  the  aquation  is  applicable  are 
related  to  Voth  temperature  and  dose  rate.  The  slanting  line  on  the  right 
represents  results  obtained  wi:!;  temperature  alone, (s)  Since  the  lines 
representing  radiation  results  extrap'Olate  across  this  line,  they  lead 
to  the  obviously  erroneous  conclusion  that  the  time  reouired  for  a  given 
amount  of  reaction  under  certain  combined  radiation  and  temperature  con¬ 
ditions  is  less  than  that  reouired  under  temperature  alone.  The 


(a)  The  line  representing  results  without  radiation  was  not  measured  In 
these  experiments,:  the  results  shewn  are  auoted  from  reference  3  and 
are  for  sli.ghtly  different  experimental  conditions,  i,e.  for  an  oxygen 
flow  rate  of  5  liters  ner  hour. 
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Oxidation-Induction  Period  -  hours 


temperature  only  line,  therefore,  is  the  upper  temperature  limit  of  the 
equation's  applicability;  it  is  unfortunate  that  this  linrdt  cannot  be  de¬ 
termined  from  the  radiation  results  alone  or  from  the  equation. (b) 


The  second  limit  is  related  to  radiation  rate.  This  is  indicated  by 
the  heavy  horizontal  line  -^at  the  bottom  of  Fijmre  3  -  which  represents 
a  radiation  rate  (5.58  x  10°  ers^s  hr~^)  of  such  map;nitude  that  thermal 

effects  are  ner;li»ible.  whether  such  a  cotidition  actually  exists  is  inde¬ 
terminate  with  the  existing  data.  However,  matheraitically,  the  equation 
predicts  a  line  of  positive  slope,  i.e..  a  lesser  effect  at  a  high  tempera¬ 
ture  than  at  a  low  temrerature,  if  comnutations  are  made  for  radiation 
rates  '’ater  than  5.58  x  10°  ergs  gm'-  hr“^  for  the  oil  oxidation.  Hence, 
it  is  presumed  that  the  dose  rate  that  yields  -  by  computation  -  constant 
time  irrespective  of  temperature  is  the  dose  rate  limit. 


The  region  of  the  excessively  high  dose  rates  is  certain  to  be  of 
practical  interest  in  some  instances,  but  relationships  between  time  and 
dose  rate  in  this  region  are  strictly  a  matter  of  conjecture  at  present. 

It  might  be  presumed  that  one  would  obtain  a  series  of  horizontal  lines 
in  this  region  with  at  least  two  possibilities  for  the  time-dose-rate 
relationship.  First,  the  time  for  a  given  amount  of  reaction  may  be  in¬ 
versely  proportional  to  the  dose  rate,  i.e.  the  dosage  required  for  a 
’iven  amount  of  reaction  v^ould  be  the  same  irrespective  of  the  dose  rate, 
decond,  the  reaction  under  study  might  possess  a  finite  reaction  velocity 
occasioned  by  concentration,  mobility  of  the  reacting  species,  or  other 
factors,  such  that  much  of  the  energy  provided  by  radiation  is  not  utilized 
in  the  particular  reaction  beir.g  studied.  In  this  instance,  the  dosage 
for  a  constant  amount  of  damage  may  vary  for  different  radiation  rates. 

itadiation  Rate  and  liadiation  Dosage  -  Hadiation  dosage  rather  than 
radiation  rate,  has  been  the  basis  on  .;hich  the  limits  of  radiation  toler¬ 
ance  of  many  materials  liave  been  esta!  lished,  the  basis  being  supported 
by  many  experimental  obs.^rvations.,  The  indiratior'is  of  the  suiriaary  equa¬ 
tion  are  such,  however,  as  to  raise  the  possibility  that  the  experimental 
observations  may  be  applicable  only  under  specific  conditions.  The  rela¬ 
tionship  indicated  by  the  equation  between  time  and  radiation  rate  at 
constant  temperature  is 


In  t  -  a  In  r 


or 


t  -- 


1 


where  the  constant  "a"  is  a  function  of  temperature.  This  relationship 
is  shown  iri  Figure  4.  f.ult inlying  by  r  to  obtair^  dosage  (rt)  gives 

(b)  Presumably,  the  slope  of  Uie  ter.perature-only  line  would  be  given  by 
the  first  term  of  the  summary  equation  if  the  empirical  constants  were 
delerndned  for  the  proper  units  of  tiir.e  and  radiation  rate,  e.g.  seconds 
and  ergs  rnole"^  sec~^.  bince  temperature-ordy  experiments  were  not 
completei’.,  this  point  could  not  be  checked. 
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Therefore,  at  different  temperature's,  the  dost;:e  is  proportioral  to  a 
power  functior:  of  radiation  rate  and  straight  lines  are  obtained  when  In 
dosage  is  plotted  versus  In  ndsaticn  rate,  as  shown  in  Figure  5»  .A  re¬ 
lationship  of  this  type  has  been  reported  by  Goocir-iari  and  Coleman'^/  who 
noted  that  the  dosage-to-failure  of  Kel-F  as  a  dielectric  was  proportional 
to  the  square  root  of  radiation  rate  at  temperatures  of  30°  to  90°G. 

at  a  temperature  at  which  "a''  equals  one,  the  dosage  is  proportional 
to  a  constant,  i.e,  the  dosage  at  a  constant  amount  of  dainage  is  the  same 
irrespective  of  dose  rate,  l-ieasuremerit s  at  this  temperature  and  in  the 
region  of  "excessive"  raciiation  rate  (the  region  below  the  horizontal 
line  of  Figure  3)  appear  to  be  two  conditions  under  which  the  dosage 
correspondirig  to  a  given  amount  of  reaction  is  indeperident  of  dose 
rate. 


AfFLicHbiLirr  ivj  oTiuai  Mfriiti/LLo 


At  the  time  of  this  report,  studies  of  other  materials  in  the  frame 
of  a  reference  of  a  coristint  cunourl.  of  reaction  had  net  been  pursued  to  a 
point  that  an  over-all  equation  could  be  derived  with  confidence.  A  few 
results  obtained  with  Teflon  aircraft  engine  hoses  are  discussed  here  for 
the  comparisoi.  they  provide  wit.h  the  oil  oxidation  data. 


Test  results  are  shovai  in  Figure  6  in  which  the  time-to-failure  is 
plotted  versus  the  reciprocal  absolute  temperature.  Failure  of  the  hoses 
was  determined  by  leakage  or  rupture,  tlie  samples  being  held  under  a  static 
pressure  of  1100  psig  with  MIL-L-7f08  lubricating  oil  during  the  test. 

Host  of  the  measurements  were  made  at  a  radiation  exposure  rate  of  1,3  x 
10°  ergs  (gmC)~*  lir  )  because  of  poor  performance,  especially  in  a  subse¬ 
quent  test  in  which  the  pressure  was  cycled  (c),  only  two  additional  tests 
were  made  under  the  stauic  conditiot.s  for  the  purpose  of  checking  the  ap¬ 
plicability  of  the  equation.  The  experimental  me-surements  are  describable 
by  the  equation 


In  t 


2d79.25  ,156.43 

T  "  ~  '  T 


0.26:-'9l)  In  r  »•  7.00331 


/  \  ’  -1 

where  r  =  exposure  dose  rate,  ergs  vgmG/~‘  hr  .  Calculations  with  this 
equation  yield  the  results  shown  in  Figure  7.  ,ilthough  the  "picture"  pre¬ 
sented  is  self-cor:s:  stent ,  there  are  too  fev/  data  to  assure  its  validity. 
Certainly,  tfie  corin', tr^d  reisults  ouaiit.at:  vol y  disagree  witli  the  data  of 
iiarringtcri(  5 )  who  corcluded  from  ^^xperiirients  with  a  nimiber  of  elastomers 
and  plastics,  including  Teflon,  that,  damage  is  deuendert  only  on  total 
dosage  at  dose  rates  fron,  IcA  to  10*^  roentgen  hr"^.  There  may  be  no 


(c)  Teflor;  hoses  tested  at  J50^F  and  a  r  utiation  rate  of  1.3  x  IcA  ergs 
(gmG)~^  hi'  ‘  failed  ;fter  9.7  hours  when  the  pressure  v^as  cycled  from 
0  to  lOUU  pbig  au  13  minute  intervals. 
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Figure  6 

EXPERIMENTAL  MEASUREMENTS  OF  TIKE- TO-FAILURE  OF 
TEFLON  AIRCRAFT  ENGINE  HOSES 
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disagreement  at  hif^h  dose  rates.  As  discussed  previously,  the  damaye  may 
be  dependent  only  on  total  dosaye  if  the  degradation  reactions  are  not 
time  dependent  and  the  dose  rate  is  equal  to  or  greater  than  that  at  which 
thermal  effects  are  negligible.  In  Teflon,  this  dose  rate  is  2.74  ;  10'^ 

ergs  (gmC)“l  hr-l  ('~3.19  x  105  r  hr“l)  -  if  the  equation  is  correct.  The 
disagreement  at  the  lower  dose  rates  may  be  due  to  sample  geometry,  do¬ 
simetry,  or  other  factors;  but  the  writer  tends  to  the  belief  that  dif¬ 
ferences  in  the  mechanical  stress  levels  during  the  two  types  of  irradia¬ 
tion  experiments  are  responsible. 

Since  Teflon  has  excellent  thermal  and  oxidation  stability  under  non- 
radiation  conditions,  the  thermal  effects  apparent  in  the  experimental  re¬ 
sults  can  probably  be  ascribed  to  the  reactions  initiated  by  the  radiation. 
In  this  respect,  it  is  interesting  to  note  that  thermal  effects  in  the  Tef¬ 
lon  are  minor  in  comparison  to  those  in  the  oil  oxidation  over  the  tempera¬ 
ture  range  considered.  In  fact,  the  computed  results  for  a  dose  rate  of 
lO'^  ergs  (gmC)-l  hr~l  indicate  such  a  small  thermal  effect  that  good  experi¬ 
mental  statistics  would  be  required  for  its  detection. 


SUMMARY  AND  CONCLUSIONS 


An  analysis  of  experimental  results  has  been  described  which  indicates 
that,  where  the  time  renuired  for  a  constant  amount  of  reaction  is  the  vari¬ 
able  studied,  the  combined  time,  temperature,  and  radiation  effects  in  or¬ 
ganic  materials  can  be  described  by  an  eouation  of  the  form 

In  t  ,  +  ^3)  in  r  +  c, 

where  t  is  elapsed  time,  T  is  absolute  temperature,  r  is  radiation  dose 
rate,  and  the  c's  empirical  constants.  Several  relationships  between  time, 
temperature  and  radiation  dose  rate  ars  apparent  in  the  equation.  The  re¬ 
lationship  between  the  dosage  corresponding  to  a  given  amount  of  damage 
and  dose  rate  can  be  derived  and  is  indicated  to  be:  dosage  •'*'  rl-a,  where 
a  is  eoual  to  cp/T  +  03  and  is  temperature  dependent.  The  implications 

of  the  eouation  are  that  there  are  two  conditions  under  which  the  dosage 
for  a  given  amiount  of  damage  is  independent  of  aose  rate;  these  include  a 
specific  temperature  and  high  dcse  rates,  both  being  dependent  on  the  par¬ 
ticular  material  under  study. 

Although  the  validity  of  the  over -all  eouation  and  its  applicability 
to  a  variety  of  organic  materials  ai’e  yet  to  be  established,  it  appears 
reasonably  certain  that  combined  radiation  and  temperature  effects  can 
profitably  be  studied  in  terms  of  a  constant  amount  of  reaction.  The  con¬ 
stant  amount  of  reaction,  which  can  be  interpreted  as  the  "life"  of  materials, 
appears  to  obey,  under  constant  dose  rate  conditions,  the  logarithm  time 
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versus  absolute  temperature  relationship  familiar  in  accelerated  life  tests. 
Accelerated  life  tests  have  been  found  to  be  applicable  to  numerous  materials 
and  components  such  as  o-ring  seals,  electrical  insulation,  electric  motors, 
solenoids,  and  other  components  vdiose  failure  is  dependent  on  the  deterior¬ 
ation  of  a  "weak -link*'  material.  Such  an  approach  to  radiation  testing 
would  appear  to  offer  the  advantages  of  yielding  results  of  practical  value 
and  of  opening  up  the  possibility  of  extrapolating  experimental  results  by 
mathematical  analysis. 
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RADIATION  DAMAGE  OF  AIRPLANE 
TIRE  MATERIALS 
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SUMMARY 

An  airplane  tire  Is  a  relatively  conplex,  highly  stressed 
structure.  Radiation  damage  effects  are  unusually  conpllcated 
In  this  case  because  of  the  number  of  materials  Involved  and 
the  cooperative  effort  required  of  each.  Damage  to  tire  com¬ 
ponents  typical  of  present  production  was  studied  with  cobalt- 
60  gamma  radiation.  The  influence  of  various  environmental 
factors  was  also  Investigated.  Tire  fabric  and  fabric-rubber 
adhesion  are  far  more  susceptible  to  radiation  damage  than  are 
the  rubber  conpounds  in  general  use.  Damage  to  the  fabric  Is 
especially  critical,  in  terms  of  tire  performance,  since  It  bears 
the  brunt  of  the  stresses.  Cord  fatigue  and  dynamic  adhesion 
are  the  most  serious  aspects  of  radiation  damage  to  tires. 

Various  types  of  nylon  tire  cords  possess  slgnificcuit  differences 
in  their  ability  to  withstand  irradiation.  Presence  of  a  nitro¬ 
gen  atmosphere  was  found  to  be  very  efficacious  for  reducing 
radiolysis  in  nylon  cords.  It  was  thus  conceived  that  nitrogen 
inflation  would  provide  a  practical  means  of  raising  the  damage 
threshold  of  tubeless  airplane  tires.  Several  nitrogen  inflated 
tires  were  ejqiosed  to  a  dose  of  107  rad,  a  level  above  the 
damage  threshold,  for  nylon  cords  when  irradiated  in  air.  Their 
performance  in  laboratory  tests  confirmed  the  protective  benefits 
of  nitrogen  inflation. 


INTRODUCTION 

An  airplane  tire  is  the  end  result  of  the  skillful  blend¬ 
ing  together  of  rubber,  adhesive,  textile,  and  wire  into  the 
integrated  type  of  structxire  shown  in  Figure  1.  Its  relatively 
conplex  structure  is  subjected  to  high  stresses  and  repeated 
thermal  snocks  during  operation.  Its  performance  will  depend 
upon  the  degree  to  which  all  constituents  mutually  cooperate  in 
absorbing  these  mechanical  and  thermal  stresses. 
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Since  tires  have  provided  the  only  workable  means  to  date  for 
takeoff  and  landing  of  high  speed  land  based  airplanes,  we  quite 
naturally  assume  they  will  also  be  essential  to  a  nuclear  propelled 
plane.  In  anticipation  of  the  need  for  information  on  the  perform¬ 
ance  of  airplane  tires  operating  in  a  radiation  environment.  The 
Goodyear  Tire  &  Rubber  Company  initiated  in  its  Research  Division, 
late  in  19^6,  a  study  of  radiation  dairjige  to  typical  present  day 
tire  production  materials. 

This  work  has  been  carried  out  in  Akron  where  we  have  a  versatile 
facility  containing  3800  curies  of  cobalt-60.  Our  arrangement  com¬ 
bines  the  desirable  features  of  irradiations  in  air  with  the  ease  of 
handling  and  storage  of  the  source  in  a  water  shield.  Any  predetermin¬ 
ed  environmental  and  tenperature  conditions  can  be  maintained  during 
exposure. 


RADIATION  DAMAGE  -  GENERAL 

It  is  well  known  that  all  organic  materials  are  affected  adverse¬ 
ly  sooner  or  later  if  given  a  svifficient  exposure  to  radiation.  This 
Inpalrment  of  a  material's  ability  to  perform  its  specified  functions, 
l.e.,  a  foreshortening  of  its  useful  lifetime,  is  called  radiation 
damage.  The  minimum  cumulative  exposure  which  will  produce  a  detect¬ 
able  degree  of  damage  in  a  property  essential  to  an  article's  normal 
service  life  is  spoken  of  as  the  threshold  damage  dose.  Likewise,  a 
dose  which  will  cause  a  loss  in  such  a  property  is  termed  the 
threshold  level. 

Radiation  damage  to  tire  materials  takes  the  form  of  lowered  re¬ 
sistance  toward  abrasion,  tensile  or  conpressive  forces  and  flexure. 
Severe  damage  obviously  leads  to  premature  failure  in  service.  Thresh¬ 
old  levels  for  readily  detectable  damage  vary  for  these  constituents 
from  perhaps  10^  rad  for  textiles(^)  and  adhesives  to  2X  10°  rad  for 
rubber  conmounds.(^)  Similarly,  25^  damage  levels  might  be  10°  rad 
and  2.?xl0'  rad,  respectively. 

RADIATION  DAMAGE  -  RUBBER  COMPOUNDS 

Several  investigators  have  published  results  of  radiation  damage 
studies  on  rubber.  (2,-1^)  After  receiving  a  dose  of  about  2  meg¬ 
arads  typical  rubber  compounds  begin  to  show  decreased  tensile  strength 
and  elongation  and  increased  modulus  and  hardness.  Such  changes  in 
the  physical  properties  of  a  tread  compound  are  shown  in  Fig\u*e  2  for 
both  a  natural  rubber  and  a  (K-S  type.  This  behavior  is  quite  similar 
to  normal  deterioration  in  use,  though  obviously  accelerated,  as  shown 
by  Figure  3- 
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The  fundamental  reaction  responsible  for  ageing  of  rubber  is 
thought  to  be  oxidation  of  the  rubber  hydrocarbon  by  air. (17)  The 
rate  at  which  ageing  progresses  depends  on  the  presence  of  such 
accelerators  as  heat«  light,  moisture,  and  ozone. 

There  are  five  variables  to  be  recognized  in  any  ageing  study: 

(1)  type  of  rubber,  (2)  degree  of  cure,  (3)  presence  of  antioxidant, 

(U)  tenperature,  and  (5)  environment.  Each  of  these  can  exert  a  pro¬ 
nounced  effect  on  the  rate  of  ageing.  These  same  variables  must  be 
considered  in  any  radiation  damage  study. 

The  type  of  rubber  compound  is  an  inportant  factor  in  radiation 
damage.  Natural  rubber  and  GR-S  tend  to  harden  under  exposure,  where¬ 
as  butyl  becomes  progressively  softer.  The  precipitous  drop  in  tensile 
strength  of  butyl,  Figure  U,  admirably  illustrates  this  point.  For 
this  reason,  when  choosing  tubeless  airplane  tires  which  might  be  sub- 
;)ect  to  irradiation,  one  must  be  careful  to  select  those  whose  satis¬ 
factory  functioning  do  not  depend  upon  butyl  inibber.  Static  tests  on 
an  irradiated  tire  containing  a  butyl  gas  barrier  might  prove  quite  mis 
leading  because  the  sticky  damaged  butyl  might  still  act  as  a  gas 
barrier  under  this  particular  condition.  The  tire  would  certainly 
fall,  however,  when  tested  under  conditions  simulating  actual  operation 

The  degree  to  which  a  rubber  compound  is  vulcanized  prior  to 
Irradiation  also  governs  the  degree  of  damage  it  will  sustain  after 
a  given  dose.  Note,  in  Figure  how  the  same  natural  rubber  gum 
conpound  appears  more  radiation  resistant  when  undercured  and  less 
resistant  when  overcured  in  conparison  with  the  best  cure.  The  under¬ 
cured  sanple  actually  continues  essentially  to  cure  or  cross-link 
during  the  early  exposure  period,  reaching  its  peak  tensile  after  50 
megarad.  Because  of  this,  conparisons  between  conpounds  may  be  more 
difficult  than  one  would  normally  think. 

The  presence  of  an  antioxidant  or  antlozonant  in  the  conpound 
under  study  may  have  a  significant  bearing  on  its  ability  to  with¬ 
stand  irradiation.  Such  damage  inhibitors  have  been  the  subject  of 
a  continuing  investigation  by  Born  and  coworkers.  (8»9»10)  We  are 
inclined  to  think  that  chemicals  which  act  as  radiation  damage  in¬ 
hibitors  conform  to  the  normal  mechanism  of  protection  offered  by 
antioxidants  or  antiozonants  since  they  usually  belong  to  those 
classes  of  chemical  structures  known  to  possess  such  properties. 

Radiation  damage  to  most  rubbers  is  quite  tenperature  dependent. 
Note  in  Figure  6  how  damage  to  a  natural  rubber  tread  stock  increases 
rapidly  in  severity  as  the  tenperature  during  exposure  is  increased. 
This  response  lias  also  been  noted  with  nitrile  and  chloroprene  rubbers. 
A  new  tyne  of  rubber  developed  by  Goodyear,  called  Adduct  rubber. 

Figure  7,  however,  apnears  to  possess  an  improved  resistance  toward 
radiation  damage  at  room  tenperature  and  a  surprisingly  lack  of  any 
particular  tendency  toward  tenperature  dependence  throughout  the 
range  of  -120 *F  to  200 ’F. 
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The  envirotiment  in  which  a  rubber  is  e^osed  can  also  be  in¬ 
fluential  Iri  the  matter  of  radiation  damage.  Gamma  radiation  pro¬ 
duces  ozone  in  air  and  this  reactive  gas  is  perhaps  rubber's  worst 
enemy.  It  is  especially  dafnaging  when  the  compound  is  under  stress. 
We  have  found  that  this  type  of  surface  attack  can  be  prevented  by 
flushing  with  an  inert  gas  prior  to  irradiation  and  holding  the 
sample  in  such  a  protective  atmosphere  during  exposure.  Figure  8 
Indicates  the  degree  of  protection  one  can  obtain  with  this  method 
when  irradiating  a  natural  rubber  tread  stock. 

Rubber's  dynamic  properties  also  suffer  during  irradiation. 
Johnson,  et.  al.,(^5)  found  the  flex  life  of  a  natural  rubber 
vulcanizate  to  decrease  by  affer  a  50  megarep  dose  and  that  of 
a  GR-S  vulcanizate  to  drop  by  86;?. 

RADIATION  DAMAGE  -  FABRIC-RUBBER  ADHESITO 

Radiation  damage  to  adhesives  of  the  type  enployed  in  bonding 
textiles  to  rubber  offers  a  virgin  field  of  exploration.  At  least, 
one  must  so  conclude  from  the  paucity  of  literature  on  the  subject. 
The  adhesive  which  bonds  rubber  to  cord  performs  a  very  critical 
function  in  a  tire.  As  such,  Its  susceptibility  to  radiation  damage 
must  not  be  overlooked. 


We  have  a  few  preliminary  tests  on  two  proprietary  adhesives 
coded  A  and  B.  '^able  I  illustrates  how  the  dynamic  adhesions  between 
rubber*  and  various  tire  cords  behave  after  irradiation.  Data  in  the 
second  and  third  columns  were  obtained  from  test  pieces  which  had 
been  built  and  then  irradiated,  whereas  those  in  columns  four  and  five 
resulted  from  specimens  in  which  the  adhesive  coated  cord  was  erqjosed 
before  the  test  pieces  were  made  up.  The  greater  damage  shown  by  the 
latter  type  treatment  must  be  due  to  radiation  induced  oxygen  attack 
during  the  e:3q303ure  in  air.  After  a  lU  megarad  dose  nylon  shows  a 
high  loss  in  adhesion,  rayon  shows  a  small  loss,  and  Dacron,  with  poor 
initial  adhesion,  retains  virtually  none  after  exposure.  Table  II 
illustrates  what  happens  when  a  similar  approach  is  followed  in  study¬ 
ing  the  loss  of  adhesion  through  hot  static  tests.  Here  again  the 
samples  made  up  from  exposed  cords  showed  greater  damage.  The  degree 
of  damage  again  increases  in  the  order  Dacron,  rayon  and  nylon. 


In  summarizing  the  liata  in  these  two  tables  we  can  say  that  after 
a  dose  of  the  order  of  15  megarad,  (1)  rayon  retains  most  of  its  ad¬ 
hesion  strength, ^2)  nylon  shows  a  large  loss  in  adhesion  but  may  still 
be  in  the  "get  by"  range,  and  (3)  Dacron,  which  is  difficult  to  bond 
to  rubber,  retains  virtually  no  adhesion.  Also,  dynamic  properties  are 
more  sensitive  to  radiation  damage  than  are  static  properties. 
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TABI£  I 


Effect  of  lU  Msgarad  Qanma  Ray  Esqposxire 
on  Dynamic  Strip  Adhesions  of  Tire  Cords 

(Rate  of  Separation)"^ 


Type 

Cord 

Type 

Ad¬ 

hesive 

Con¬ 

trol 

Separation 
Ex-  Rate 

posed-aa  Increase  % 

Ex¬ 
po  sedaaa 

S^>aration 

Rate 

Increase  % 

Nylon  700  L 

A 

5.55 

2.12 

161 

0.72 

672 

Nylon  700 

A 

5.55 

2.12 

161 

0.55 

555 

Rayon  T258 

A 

5.55 

5.00 

11 

U.55 

22 

Rayon  T258 

B 

1.39 

1.16 

18 

1.39 

0 

Dacron  ^100 

A 

0.6U5 

virtually  no  adhesion 

*  subjected  to  3  lbs.  @  250*F 
-a*  test  pieces  built  before  Irradiation 

test  pieces  built  after  irradiating  dipped  cords 
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TABIE  II. 

Effect  of  lli  Megarad  Qamna  Ray  Exposure 
on  Hot  Static  Adhesion  of  Tire  Cords 

Hot  Static  Adhesion,  lbs  &  2^0 *F 

Type 


Type 

Ad¬ 

Con¬ 

Ex¬ 

Adhesion 

Ex¬ 

Adhesion 

Cord 

hesive 

trol 

posed* 

Loss  a 

posed** 

Loss  % 

Nylon  700L 

A 

11.0 

9.3 

15. U 

7.5 

31.8 

!^lon  700 

A 

13.1 

9.0 

31.5 

3.0 

77.0 

Rayon  T258 

A 

17.1 

15.1 

11.8 

IU.5 

15.2 

Rayon  I&$8 

B 

18.3 

15.5 

15.3 

13.5 

26.2 

Dacron  '5100 

A 

U.o 

U.l 

0 

5.5 

♦37.5 

*  test  pieces  built  before  Irradiation 
test  pieces  built  after  irradiating  dipped  cords 
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RADIATION  DAMAGE 


TIRE  FABRIC 


Radiation  damage  to  tire  fabric  is  especially  critical,  in  terms 
of  performance,  since  this  material  bears  the  brunt  of  all  stresses  to 
which  the  tire  is  subjected.  We  believe  this  to  be  the  factor  which 
limits  a  tire's  usefulness  when  operating  in  a  radiation  environment. 
This  is  the  weak  link  in  the  chain.  Currently  used  textiles  may  turn 
out  to  be  inadequate  because  of  poor  radiation  resistance.  Other  cord- 
type  materials  tetter  able  to  withstand  irradiation  have  never  teen 
used  successfully  in  tires  due  to  poor  tension- compression  fatigue  amd/ 
or  poor  adhesion.  Therefore  the  major  emphasis  at  this  time  should  be 
placed  on  iii;)roving  the  radiation  resistance  of  tire  fabric  rather  than 
tire  confounds. 

Quite  a  few  publications  have  dealt  with  effects  of  radiation  on 
textiles,  some  of  which  are  used  in  tire  fabrics(l*  Id*  19,  20).  When 
estimating  radiation  damage  to  a  synthetic  textile,  one  must  be 
cautioned  against  assuming  that  data  obtained  from  molded  or  pressed 
samples  of  the  polymer  will  be  the  equivalent  of  data  gathered  on 
fiberous  type  specimens.  The  large  surface  to  volume  ratio  of  fibers, 
the  highly  oriented  condition  in  which  they  are  drawn,  and  the  twisted 
state  in  which  cords  are  normally  used  all  contribute  to  an  increased 
sensitivity  to  radiation  damage.  Figure  9  demonstrates  the  differences 
one  can  expect  in  stress-strain  properties  of  nylon  cords  and  nylon 
molded  sheets.  Note  how  the  two  cases  in  many  respects  show  exactly 
apposite  trends. 

There  are  also  five  valuables  which  should  be  Investigated  in  a 
thorough  analysis  of  radiation  damage  to  tiro  fabrics:  (l)  type  of 
fabric,  (2)  environment,  (3)  stress,  (h)  temperature,  and  (5)  effect 
of  antioxidants  and  antlozonants.  We  have  a  study  underway  to 
evaluate  the  separate  and  combined  effects  which  these  variables  play 
during  irradiation  of  tire  fabrics. 

Different  fabrics  possess  varying  resistances  to  radlolysls. 

Even  different  types  of  the  same  fabric  show  this  phenomenon.  Note 
in  Figure  10  how  the  relative  resistance  varies  with  fabric  material 
and  with  type  of  nylon  as  well.  Since  an  airplane  tire  may  have  to 
withstand  a  cumulative  dose  somewhat  greater  than  107  rad,  we  have 
exposed  these  cords  to  doses  in  the  range  of  to  UOO  megarads. 

Cotton  and  nylon  are  virtually  destroyed  by  a  100  megarad  dose  where¬ 
as  Dacron  suffers  only  a  ll:^  loss  at  twice  this  dose.  Since  types 
66,  300,  and  700  nylon  are  the  same  chemically  and  vary  only  in  such 
respects  as  molecular  weight,  crystallinity,  surface  finish,  and 
antioxidant  type  and  concentration,  the  sigiilficant  differences  in 
their  radiation  sensitivity  is  most  Interesting.  This  subject  will 
be  considered  shortly. 
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The  most  serious  aspect  of  radiation  damage  to  tire  fabric  lies 
in  cord  fatigue.  Table  III  shows  that  nylon  and  rayon  suffer  virtual¬ 
ly  a  conplete  loss  in  ability  to  withstand  repetitjve  tension- 
conpresalon  stresses.  A  cord  which  will  stand  about  1  million  cycles 
is  considered  acceptable  in  this  respect.  Dote  that  Dacron  is  poorer 
initially  than  nylon  after  a  dose  of  16  megarads.  These  data  are  in 
substantial  agreement  with  those  published  by  HarmonC^). 

The  environment  in  which  a  cord  is  exposed  is  definitely  related 
to  the  damage  it  suffers.  Figure  11  compares  damage  to  type  3CX)  nylon 
in  air  to  that  when  esqjosed  in  nitrogen.  If  a  30^  loss  in  initial 
tensile  amd  elongation  can  be  tolerated,  a  tenfold  inci*ea8e  in  life¬ 
time,  dosagewise,  can  be  gained  by  bathing  nylon  in  nitrogen  during 
irradiation.  The  8i;perior  retention  of  properties  of  fabrics  exposed 
in  vysuum  over  those  treated  in  air  has  been  reported  in  the  litera¬ 
ture  (l).  Radiation  induced  oxidation  is  very  detrimental,  as  shown 
by  the  data  in  Table  IV.  A  1-25^  drop  in  stress-strain  properties 
occurs  in  nylon  type  300  exposed  to  ^  megarad  in  nitrogen  ,  and  a 
3-5^  drop  occurs  if  the  cords  are  under  stress  during  this  period. 
However,  cords  e:q)08ed  in  air,  at  relatively  low  dose  rates  in  order 
to  enphasize  oxidative  effects,  show  losses  of  when  unstressed 

and  69- 71^  if  stressed. 

The  degree  of  radiation  damage  to  tire  cords  depends  very  much  on 
the  magnitude  of  both  temperature  and  stress  during  irradiation.  We 
are  currently  investigating  the  combined  effects  of  dose,  teirperature, 
stress,  and  environment.  Some  preliminary  data  am  listed  in  ^able  V. 
Controls  held  in  air  showed  some  degradation,  especially  when  stressed; 
but,  samples  exposed  to  28  megarad  were  destroyed.  A  specimen  held  at 
200*F  under  a  stress  of  2  kg  per  cord  failed  statically  at  a  dose  less 
than  26  megarad.  These  tests  are  being  repeated  in  nitrogen  to  clarify 
the  overall  Irportance  of  oxygen  in  this  damage. 

If  the  previously  mentioned  ability  of  some  nylons  to  resist 
damage  is  due  to  varying  types  or  concentrations  of  antioxidants, 
either  incorporated  throughout  the  polymer  or  added  only  to  the  sur¬ 
face  finish,  it  is  conceivable  that  yet  greater  protection  by  this 
means  is  still  possible.  To  test  the  validity  of  this  hypothesis  we 
surface  coated  nylon  cords  •.d.th  several  antioxidants  prior  to  irradia¬ 
tion.  These  did  not  demonstrate  any  significant  degree  of  protection. 
It  may  be  necessary  to  actually  innoculate  nylon  with  an  appropriate 
antioxidant  during  manufacture  in  order  to  get  effective  protection 
by  this  means. 

Compounding  rubber  stock  in  which  fabric  is  to  be  cured  with  a 
damage  inhibitor  of  the  antioxldant-antiozonant  type  did  confer  some 
protection  to  the  fabric  in  some  of  our  experiments.  Figure  12  in¬ 
dicates  the  percent  of  initial  tensile  strength  retained  by  adhesive 
dipped  nylon  700P  cords  when  subjected  to  Irradiation  while  bare, 
while  covered  with  a  natural  rubber  vulcanlzate  lacking  an  antioxidant. 
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TABLE  ni. 

Damaging  Effect  of  16.  Magarad  Qaimna  Ray  Exposure 
on  Fatigue  Characteristics  of  Tire  Cords 


Number  Test  Cycles  Required 
for  Failure* 


Type 

Type  Ad- 

Cord  heslve 

Control 

Exposed 

Loss 

in 

Fatigue 
Resistance,  % 

Nylon  700  L  A 

2,205,500 

56,900 

97.2 

Nylon  700  A 

2,672,600 

58,867 

97.7 

Rayon  T258  A 

553,700 

2,125 

99.6 

Rayon  T258  B 

1,186,550 

1,975 

99.8 

Dacron  5100  A 

20,825 

22,1450 

+7.0 

*Mallory  Fatigue  Tester  operating  at  8^0  RPH  with 
So  lbs*  inflation* 
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RADIATION  DOSE 

Effect  of  environment  on  radiation  damage  to 
exposed  at  75 "F. 


TABLE  IV. 


Effect  of  Environment  on  Damage  to  Type  300 
Hylon  Tire  Cords  Induced  by  $  Negarad  e^qposure 
to  Co-60  Oairnna  Rays  at  7$*F 


Experimental  Conditions  Stress-Strain  Properties 


Time 

Required 
for  Ex¬ 
posure  days 

Environ¬ 

ment 

Stress, 

gev^cord 

Ultimate 

Tenalle, 

pal. 

%  of 
Initial 
Tensile 

Elonga¬ 
tion  at 
Break 
% 

%  of 
Initial 
Elonga¬ 
tion 

90  (control) 

air 

0 

28.0 

100 

17.5 

100 

90  (control) 

air 

1600 

28.7 

102 

17.9 

102 

90 

air 

0 

10.1 

36 

8.1 

I46 

90 

air 

1600 

8,2 

29 

5.5 

31 

5 

nitrogen 

0 

27.14 

98 

17.3 

99 

5 

nitrogen 

1600 

26.5 

95 

17.0 

97 

TXBI£  7. 


Iffeet  of  Tenperaturo  and  Stress  on  Badiation 
Damage  to  Nylon  700M  Exposed  In  Air 

Experimental  Conditions  Stress-Strain  Properties 


Dose> 

negarad 

Tern- 

per- 

atiire 

•F 

Stress 
kg  per 
Cord 

Ultimate 

Ten-  %  of 

slle  Initial 
psi.  Tensile 

Elonga¬ 
tion  at 
Break  % 

%  of 
Initial 
Elongation 

0 

80 

0 

28.6 

100 

21.3 

100 

0 

80 

2 

28.7 

100 

19.3 

91 

0 

200 

0 

27.0 

9U 

20.8 

98 

0 

200 

2 

27.0 

9k 

16.5 

78 

28 

80 

0 

10.1 

35 

10.1 

U7 

lU 

lliO 

1 

8,0 

30 

7.7 

36 

28 

80 

2 

5.U 

19 

U.3 

20 

28 

200 

0 

2.7 

9 

2.5 

12 

Less 
than  28 

200 

2 

broke 

during  exposure 

S  S  S  8 

S! 

HXONSUdLS  aUSN^I  IVIJINI  iO  % 
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RADIATION  DOSE  (MEGARAD) 

Effect  on  radiation  damage  to  nylon  cords  of  antioxidant  added  to 
rubber  coi!f>ound  used  to  coat  fabric. 


and  while  covered  with  similar  rubber  compounds  containing  differ¬ 
ent  antioxidants*  This  method  did  not  prove  as  effective  in 
preventing  damage  to  nylon  cords  as  did  the  use  of  a  nitrogen 
atmosphersi  since  with  these  antioxidants  present  a  ^0^  loss  In 
tensile  still  occurred  whereas  in  a  nitrogen  environment  only  a  20^ 
loss  was  detectable. 

RADIATION  DAMAGE  -  TIRES 

Having  observed  the  efficacy  of  a  nitrogen  environment  in  re¬ 
ducing  radiolysis  in  nylon  cords,  it  was  conceived  that  nitrogen 
inflation  would  provide  a  practical  means  of  raising  the  damage 
threshold  of  tubeless  airplane  tires.  The  gas  with  which  a  tube¬ 
less  tire  is  Inflated  normally  diffuses  to  some  extent  through  the 
gas  barrier  and  into  and  through  the  fabric  layers.  Cords  In  a  tire 
Inflated  with  nitrogen  should  thus  be  bathed  in  nitrogen  after  a 
sufficient  storage  period. 

On  the  assunptlon  that  doses  of  from  10  to  100  megarad  may  be 
realistic  service  requirements,  we  have  been  exposing  nitrogen  In¬ 
flated  tubeless  airplane  tires  to  levels  within  this  range.  Our 
experimental  setvqp  is  shown  in  Figure  13.  The  tires  are  rotated 
during  esqposure  in  order  to  improve  the  uniformity  of  the  dosage 
they  accumulate. 

Tires  e:]qposed  to  10^  rad,  a  level  above  the  damage  threshold 
for  nylon  cords  when  irradiated  in  air,  made  over  200  high  qpeed 
landings  smd  gave  a  normal  static  burst  test.  Their  performance 
thus  confirmed  the  protective  benefits  of  nltrogsn  inflation  during 
irradiation. 
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RADIATION  EFFECTS  ON  FLIGHT  CONTROL 
SUBSYSTEM  DESIGN 

by 

D  0.  Gunson 

Lockheed  Aircraft  Corporation 
Marietta,  Georgia 


Successful  application  of  nuclear  proimlsion  to  aircraft 
requires  the  development  of  a  high  performance  flight  control 
subsystem  that  is  not  only  more  reliable,  serviceable,  and 
maintainable  than  the  best  flying  today  but  is  at  least  as 
efficient  while  operating  in  a  radiation  environment <  The  at¬ 
tainment  of  these  aims  demand  close  and  continued  cooperation 
between  the  radiation  effects  specialist  and  controls  designer 
throughout  the  design,  develoiment,  and  testing  of  this  sub¬ 
system.  The  basic  problems  considered  ares  the  establishment 
of  the  important  characteristics  of  a  flight  control  subsystem, 
how  radiation  environment  affects  these  characteristics  and  the 
type  of  radiation  effects  data  required  by  the  controls  designer. 

A  practical  control  subsystem  is  as  fundamental  to  a  useful  flying 
machine  as  the  lifting  surfaces  and  power  plants,  Upon  its  development  has 
hinged  the  success  of  every  advancement  in  aircraft  performance.  Its  role 
will  not  diminish  in  the  future,  for  a  flying  machine  is  of  no  value  unless 
it  can  be  controlled  under  all  Msign  flight  conditions  v/ith  sufficient 
accuracy  to  perform  the  desired  i^ssion.  Thus,  each  advancement  in  aircraft 
impose  increased  requirements  on  l^e  control  system.  During  the  past  decade 
and  a  !ialf,  the  requirements  for  dytaanic  performance  have  advanced  at  an 
ever  increasing  rate.  Now  an  entir^y  new  environment  and  increased  flight 
times  eire  being  added  to  the  already 'stringent  requirements. 

The  successful  application  of  nuclear  propulsion  to  aircraft  and  its 
full  utilization  requires  the  development  of  a  completely  computable  fligr.t 
control  subsystem  that  is  also  tolerant  ^  the  radiation  environment  to 
which  it  will  be  subjected  under  operating  conditions  -  This  subsystem  must 
have  as  high  performance  capabilities  as  aAy  flying  today  and  be  more  re¬ 
liable,  serviceable,  and  maintainable-  This  is  a  big  order  for  the  controls 
designer. 

Since  the  radiation  environment  ai’fects  the  physical  configuration  of 
the  subsystem  as  well  as  the  chemical  and  physical  properties  of  the  com¬ 
ponent  materials,  close  and  continued  cooperation  between  the  radiation 
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Radiation  Effects  on  Flight  Control 
Subsystem  Design  (cont-) 


effects  specialist  and  the  controls  designer  is  necessary-  Each  must  have 
a  thorough  understanding  in  mutual  areas  and  be  conversant  v;ith  the  problems 
of  the  other.  The  controls  designer  must  specify  those  characteristics  and 
material  properties  that  are  important,  specify  the  test  programs,  and  eval¬ 
uate  the  test  results  as  to  subsystem  performance.  The  radiation  effects 
specialist  must,  in  turn,  evaluate  the  subsystem  under  dynamic  conditions 
in  the  radiation  environment  for  changes  in  these  characteristics  and  proper¬ 
ties.  He  must  also  aid  in  the  selection  of  mterials  and  components  and  per¬ 
form  pre  and  post  Irradiation  analyses. 

Even  though  a  considerable  amount  of  radiation  effects  data  on  applic¬ 
able  materials  and  components  is  presently  available,  it  is  incomplete  from 
an  engineering  view  point.  Most  of  these  data  were  measured  under  static 
conditions.  The  amount  of  data  taken  under  dynamic  conditions  is  small  and 
only  sufficient  to  show  that  there  is  a  great  difference  between  the  effect 
on  materials  and  components  alone  under  static  conditions  and  the  same  items 
in  a  subsystem  under  dynamic  conditions.  In  many  cases  the  effects  for  the 
later  case  are  much  less  severe.  If  a  practical  flight  control  subsystem  for 
a  nuclear  powered  airplane  is  to  be  developed,  the  actual  changes  in  material 
properties  and  physical  characteristics  must  be  evaluated  under  realistic 
operating  conditions  which  include  other  environments. 

At  the  present  time,  some  programs  are  being  undertaken  to  determine 
certain  effects  iinder  dynamic  conditions  on  small  representative  subsystem 
test  panels.  Other  and  larger  programs  are  being  planned.  If  these  programs 
are  to  produce  the  greatest  possible  gain  in  basic  system  knov/ledge,  both 
the  radiation  effects  specialist  and  controls  designer  must  actively  partici¬ 
pate.  Thus,  it  is  appropriate  to  point  out  the  more  important  characteristics 
of  control  subsystems,  how  a  radiation  envij.cnm®nt  affects  these  characteris¬ 
tics,  the  type  of  data  needed  by  the  controls  designer  and  some  possible 
methods  of  overcoming  the  problems. 

The  more  important  of  control  subsystem  characteristics  can  be  summa¬ 
rized  as  listed  in  Figure  1. 


FIGURE  1 

CIU'JIACTERISTICS 

1. 

Controllability 

2 

Stability 

3- 

Reliability 

4. 

Safety 

5. 

Versitility 

6 

Serviceability 

7. 

Maintainability 

Radiation  Effects  on  Flight  Contro] 
Subsystem  Design  (cont  ) 


This  list  could  be  expanded  greatly,  but  it  is  felt  that  the  other  pertinent 
characteristics  can  be  suitably  defined  as  an  attribute  of  one  or  more  of 
the  basic  seven. 

Controllability  is  the  measure  of  the  pilot  effort  required  for  controll¬ 
ing  the  aircraft  in  the  manner  and  with  the  accuracy  required  to  perform  the 
mission  To  free  the  pilot  from  giving  continual  attention  to  flying. 
stability  is  required.  That  is,  the  subsystem  mst  respond  only  to  the  pilot's 
command.  Reliability  is  an  absolute  necessity  and  requires  no  definition  in 
view  of  the  emphasis  presently  being  placed  on  it  The  primary  aim  must  be 
toward  substantially  failure  free  subsystems  and  components.  Safety  is  an 
Important  characteristic  which  mist  be  provided  through  "fail-safe”  design 
techniques  and  elimination  of  materials  that  retain  high  levels  of  activation. 
This  recognizes  that  due  to  damage,  hunmn  errors,  or  material  failures 
emergencies  do  occur,  but  they  must  be  so  controlled  that  safe  flight  can  be 
maintained  even  after  an  emergency  occurs  Safety  of  ground  personnel  during 
service  and  maintenance  periods  mist  also  be  provided ■  Versatility  is  a 
measure  of  the  ease  with  which  additional  features  can  be  incorporated  into 
a  subsystem  without  undue  duplication  of  functions  and  complexity  of  design 
anc*  installation  Serviceability  and  maintainabil ity  are  measures  of  the 
time  required  by  ground  personnel  for  servicing  and  maintaining  the  subsystem- 
These  last  factors  have  assumed  primary  importance  not  only  due  to  the  desire 
for  a  maxiimim  aircraft  utilization,  but  because  of  the  radiation  environment 
in  which  most  of  the  flight  controls  subsystem  components  are  located.  Host 
of  the  components  will  have  residual  activation  after  having  been  irradiated 
on  the  aircraft  All  of  these  characteristics  must  be  attained  by  design, 
development,  and  tests  of  the  subsystem  and  its  components. 

The  accomplishment  of  these  characteristics  involve  problems  that  can  be 
categorized  as  shown  in  Figure  2 


FIGURE  2 


PROBLEH  TYPES 


1  Human  Factors 

2  •  Mechanical 

3  Hydraulic 
4 ■  Electrical 

5  Aerodynamic 

6  Support  Equipment 

Human  factors  and  aerodynamics  establish  the  requirements  for  controllability 
and  stability.  Human  factors  and  support  equipment  problems  are  predominant 
in  obtaining  serviceability  and  maintainability.  Mechanical,  hydraulic,  and 
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Radiation  Effects  on  Flight  Control 
Subsystem  Design  (oont.) 


electrical  design  problems  are  involved  in  accomplishing  all  the  character¬ 
istics  . 

Existence  of  the  radiation  environment  affects  the  characteristics  and 
complicates  the  problems  in  the  manner  shewn  in  Figure  3- 

FIGURE  3 

RADIATION  EFFECTS 
1,  Changes  in  Material  Properties 
2t  Requirements  for  Remote  Handling 
3.  Increased  Danger  to  Personnel  and  Equipment 

Changes  in  certain  of  the  chemical  and  physical  properties  of  the 
mechanical,  hydraulic,  and  electrical  component  materials  affect  the  controll¬ 
ability,  stability,  and  reliability  characteristics  of  the  subsystem.  Changes 
in  friction,  hardness,  rigidity,  strength,  and  dimensions  of  the  mechanical 
components  are  of  primary  importance  to  the  dynamic  performance  of  the  sub¬ 
system.  For  the  same  reason,  changes  in  the  bulk  modulus,  visocity,  and 
lubricity  of  the  hydraulic  fluid  are  importeint.  Changes  in  such  properties 
as  flexibility,  compliance,  dimensions,  and  friction  coefficient  are  important 
in  hydraulic  seals.  Such  properties  of  electrical  components  as  resistance 
of  insulation,  conductivity  of  wires,  magnetic  characteristics  of  core  materials, 
and  certain  operating  characteristics  of  components  are  important.  These  lists 
can  be  greatly  expanded  by  being  more  specific-  The  radiation  effects  engineer 
must  determine  the  specific  materials  and  significant  properties  through  close 
cooperation  with  the  controls  designer.  The  important  point  to  be  made  is 
that  not  eiU  chemical  and  physical  properties  of  any  material  are  significant 
in  a  particular  application 

Reliability  is  of  prime  importance  to  the  controls  designer  It  csin 
be  achieved  only  through  the  assistance  of  the  radiation  specialist  in  select¬ 
ing  the  materials  whose  significant  properties  are  most  tolerant  to  radiation 
for  use  in  the  control  subsystem  design,  and  through  developmental  and  reli¬ 
ability  testing  under  dynamic  conditions  in  realistic  environments. 

Radiation  from  the  power  plants  and  activation  of  the  component  materials 
affect  the  safety,  serviceability,  and  maintainability  characteristics  of  the 
subsystem.  The  equipment  must  be  serviced  and  maintained  v;hile  on  the  air¬ 
craft  as  well  as  installed  and  removed  by  gi'ound  personnel  vfhose  safety  must 
be  assured  at  all  times.  Thus,  the  human  factor  establishes  the  requirements 
for  remote  checkout  provisions  and  test  equipment,  remote  servicing  provisions 
and  equipment,  designs  for  quick  installation  and  removal  of  the  subsystem 
components,  and  for  the  remote  handling  equipment  for  the  installation  and 
removal  of  those  components  most  susceptible  to  failures. 
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Radiation  Effects  on  Flight  Control 
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Sfifety  of  flight  requirements  are  affected  by  the  radiation  environ¬ 
ment  through  the  effects  on  equipment  reliability  since  they  are  net  primarily 
through  design.  This  considers  the  human  factors  as  related  to  the  flight 
crew.  The  aim  is  to  provide  safe  control  even  after  a  primary  failure  in 
the  subsystem. 

From  analyses  of  the  available  data,  the  apparent  problem  areas  are 
as  shown  in  Figure  4. 


FIGURE  L 

KNOWM  PROBLEM  AREAS 

1.  Piercing  of  Radiation  Shields 

2.  Lubrication 

3.  Hydraulic  Fluid 

4c  Hydraulic  Seals 

5.  Electrical  Insulation 

6.  Remote  Checkout  and  Malfunction  Detection 

7 .  Activation 

All  of  the  above  problem  areas  except  activation  and,  to  some  extent, 
remote  checkout  and  malfunction  detection,  are  similar  to  those  faced  in 
the  designs  for  very  high  temperature  applications..  In  most  cases  the  same 
solutions  are  applicable.  Many  of  the  techniques  being  employed  by  Lockheed 
were  developed  in  experimental  programs  on  high  temperature  hydraulic  and 
pneumatic  subsystems. 

Crew  shield  penetration  will  be  accomplished  by  the  use  of  torque  tubes 
and  linkages  as  shown  in  Figure  5  With  this  method,  a  minlTmim  of  free 
path  area  for  pen etT»ation  of  radiation  can  be  provided.  The  free  path  area 
consists  only  of  that  occupied  by  the  torque  tube  itself  plus  the  clearance 
around  the  outside  diameter  of  the  tube.  By  the  use  of  shield  rings  on 
each  end  of  the  tube,  the  clearance  area  can  be  effectively  reduced  to  help 
maintain  the  level  of  radiation.  This  method  also  adds  a  minimi  mi  of  friction 
to  the  controls  subsystem. 
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FIGURE  5 


Friction  always  is  an  important  problem  in  control  subsystems.  It 
increases  the  pilot  effort  required,  decreases  the  sensitivity  and  positional 
accuracy  obtainable,  and  has  a  destabilizing  effect,  Du®  to  the  long  flight 
times  allowable  in  nuclear  powered  aircraft,  the  pilot  effort  required  for 
satisfactory  control  must  be  a  minimum.  Hence,  adequate  lubrication  must  be 
provided.  To  lessen  the  reliance  on  lubrication,  the  features  shown  in 
Figure  6  must  be  employed  in  the  subsystem  design. 

FIGURE  6. 

FRICTION  REDUCTION  FEATURES 

1.  Rotary  Motion  through  Seals 

2.  Elimination  of  Cable  Pulleys 

3.  Large  Bellcranks  and  Quadrants 

4.  Straight  Control  Runs 

All  of  these  features  employ  the  principle  of  increasing  the  mechanical 
advantage  over  the  friction  forces  that  cannot  be  eliminated. 

Hydraulic  fluid  remains  the  best  proven  power  transmitting  media;  there¬ 
fore,  the  subsystems  must  be  designed  to  give  the  required  dynamic  performance 
for  conditions  where  the  fluid  properties  vary  under  Irradiation,  The  primary 
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effects  anticipated  are  changes  in  bulk  modulus,  visocity,  lubricity,  and 
sludging.  During  recent  tests  on  advanced  fluids  in  dynaMc  subsystem 
test  panels,  the  changes  in  these  properties  have  been  small,  but  the  test 
conditions  have  not  been  entirely  realistic*  Therefore,  some  difficulty 
must  be  anticipated.  The  filters  normally  in  the  subsystem  are  expected 
to  be  adequate  for  removing  any  sludge  that  forms.  Most  of  the  gas  generat¬ 
ed  in  the  fluid  by  radiation  will  remain  entrained  under  working  pressures, 
and  will  become  separated  in  the  return  lines  and  reservoirs.  This  gas  can 
be  utilized  in  maintaining  reservoir  pressurization  and  bo  bled  off  to  the 
atmosphere  throu^  pressure  relief  valves  when  an  excessive  (juantlty  has 
collected. 

Hydraulic  seals  remain  one  of  the  most  troublesome  problem  areas  be¬ 
cause  of  the  almost  vmivorsal  reliance  on  elastomeric  seals,  packings,  and 
back-up  rings  in  the  past.  Only  with  the  recent  interest  in  high  tempera¬ 
ture  applications  has  development  of  other  typos  of  seals  been  luidortakon 
on  an  intensive  scale.  Various  types  of  metal  and  mineral  seals  are  being 
developed  for  the  high  temperaixire  application.  These  seals  will  be  radia¬ 
tion  tolerant  by  nature.  It  is  felt  that  within  the  near  future,  complete 
subsystems  having  no  elastomeric  seals  will  be  developed.  This  must  be 
the  aim  for  flight  control  subsystems  of  nuclear  propelled  aircraft.  In 
the  interim,  some  elastomeric  seals  must  be  used.  At  present,  Buna-N  seals 
are  proposed  as  the  best  under  irradiation. 

The  problem  of  electrical  insulations  fall  in  generally  the  same  cata- 
gory  as  hydraulic  seals.  The  materials  found  good  for  application  in  very 
hl(^  ten^eratTxres  are  generally  inorganic;  hence,  they  are  also  radiation 
tolerant.  Some  of  the  standard  insulation  materials  are  adequate  in  many 
regions  of  the  aircraft  due  to  low  radiation  levels.  Therefore,  the  air¬ 
craft  can  be  divided  into  zones  and  the  types  of  insulation  can  vary  by 
the  zones  in  order  to  minimize  the  cost. 

Since  most  of  the  hl^  strength  metal  alloys,  ball  and  roller  beaiv 
ings,  and  magnetic  core  materials  contain  percentages  of  the  heavier 
dements,  residual  aiativation  will  be  a  hazard  to  the  ground  support 
persoxmel  who  must  service  and  maintain  these  subsystems  and  components. 

This  coupled  with  the  residual  activation  of  the  aircraft  requires  that  the 
oonqponent  Installation  designs  must  allow  quick  installation  and  removal. 

In  some  cases  remote  handling  techniques  must  be  employed.  Since  it  will 
not  be  practical  to  design  all  components  for  assembly  and  disassembly  by 
remote  handling  equipment,  the  Bu;^y  and  overhaul  logistics  problem  is  also 
complicated.  Provisions  must  be  made  at  the  using  facility  for  storing 
repairable  components  until  they  cool  off  enough  for  manual  overhaul,  and  for 
repairing  these  items.  The  non-repairable  components  must  be  disposed  of  in 
accordance  with  AEC  requirements.  The  personnel  utilization  problem  is 
complicated  by  the  safe  dose  rate  allowable. 

The  reqxiirements  for  remote  checkout  and  malfunction  detection  tech¬ 
niques  and  equipment  are  no  less  than  those  for  the  largest  missiles. 


Radiation  Effects  on  Flight  Control 
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Provisions  must  be  made  for  checking  the  complete  flight  controls  subsystem 
for  function  and  malfunction  from  a  central  point  located  in  the  nose  of  the 
aircraft  since  this  is  the  only  position  inside  the  shielded  hangar.  This 
equipment  must  be  able  to  Isolate  a  weak  or  malfunctioned  part  to  the  least 

removable  component  in  order  that  a  minimum  of  exposure  time  is  required  by 

the  ground  support  personnel  for  service  and  maintenance.  The  problem  is 
further  complicated  by  the  necessity  for  continued  reliability  of  the  test 

sensors.  Failure  of  a  sensor  will  cause  at  least  as  much  work  as  an  equip¬ 

ment  failure. 

As  has  been  pointed  out,  if  the  development  of  nuclear,  propelled  air¬ 
craft  is  to  be  successful,  the  radiation  effects  specialist  and  the  controls 
designer  are  faced  with  problems  requiring  close  cooperation  and  mutual 
support  throughout  the  design,  development,  and  testing  phases  of  the  flight 
controls  subsystem  development  program.  It  is  hoped  that  you  who  are  respon¬ 
sible  for  the  radiation  effects  now  have  a  somewhat  better  appreciation  of 
the  problems  and  requirements  for  this  particular  subsystem. 
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PNEUMATICS  -  A  TOOL  FOR  THE  DESIffllER 
OF  NUCLEAR  POl/ERED  AIRCRAFT 


by 

J.  A.  Osterman 

Lockheed  Aircraft  Corporation 
llarietta,  Georgia 

Pneumatics  has  been  demonstrated  to  have  promise  in  the 
field  of  high  temperature  aircraft  system  operation.  Air  has 
greater  resistance  to  nuclear  radiation  and  is  capable  of 
operating  over  a  wider  range  and  higher  temperatures  than  any 
available  hydraulic  fluid.  This  paper  gives  the  essentiail 
results  of  a  program  of  study  of  a  high  temperature  pneumatic 
system  which  offers  a  promising  tool  to  the  nuclear  aircraft 
system  designer. 

The  susceptibility  of  hydraulic  fluids  to  radiation  damage  has  been 
one  of  the  most  important  problems  facing  the  aircraft  engineer  who  is 
charged  with  designing  control  systems  for  nuclear  powered  aircraft.  It  is 
only  natural  therefore,  that  considerable  attention  has  been  directed  toward 
the  possible  use  of  pneiimatics  as  providing  a  service  fluid  undamaged  by 
nuclear  radiation. 

The  sane  approach  has  been  used  to  solve  the  high  temperature  problem 
in  aircraft  system  design.  The  aircraft  engineer  is  faced  more  and  more 
by  higher  and  higher  temperatures  which  are  beyond  the  operating  limits  of 
hydraulic  fluids.  Here  again,  pneumatics  has  been  given  an  appredsal  with 
the  view  of  providing  a  working  fluid  compatible  with  high  temperatiu^e , 

For  the  past  two  years,  Lockheed  Aircraft  Corporation,  imder  contract 
with  Wright  Air  Development  Center,  has  been  engaged  in  the  field  of  high 
temperatiire  pneumatics.  I^eclfically  this  effort  has  been  aimed  at  the 
design,  development  and  test  of  a  5000  pound  per  square  inch  actuator  for 
lOOO^F  operation.  The  purpose  of  this  jjaper  is  to  present  the  results  of 
this  program  v;hich  adds  a  promising  tool  not  only  in  the  field  of  high 
temperature  but  in  nuclear  engineering  design. 

The  initial  effort  was  embodied  in  a  study  program  to  determine  the 
feasibility  of  pneumatics  and  to  design  an  aircraft  speed  brake  servo 
system  to  operate  in  ambient  temperatures  up  to  lOOOOF.  Two  basic  points 
of  philosophy  were  germane  to  the  program  right  from  the  beginning.  Since 
it  was  felt  that  much  of  the  failure  in  other  pneumatic  programs  arose 
because  the  "hardware”  consisted  of  modified  hydraulic  components,  it  was 
decided  that  we  should,  in  effect,  "design  for  pneumatics".  Secondly,  with 
the  high  operating  temperatures  particularly  in  mind,  it  was  decided 


to  use  the  best  naterials,  finishes,  and  surface  treatments  available. 

This  approach  was  taken  because  the  number  of  critical  design  problems 
demanded  the  utmost  care  in  dealing  with  each  one  so  that  all  possible 
"weak  links"  could  be  eliminated. 

Once  the  program  philosophy  had  been  established,  the  study  began  with 
a  comparison  of  various  gases.  While  the  term  pneumatics  generally  implies 
air  in  its  gaseous  state,  it  was  defined  here  to  include  all  likely  gases. 

The  properties  and  characteristics  of  thirteen  condidates  v/ere  studied  and 
compared.  The  list  included  air,  ammonia,  argon,  carbon  dioxide,  carbon 
monoxide.  Freon,  helium,  hydrogen,  neon,  nitrogen,  oxygen,  steam,  and  sul¬ 
phur  dioxide.  In  addition  to  the  physical  properties  of  gases  such  as 
bulk  modulus,  boiling  point,  specific  heat,  and  density,  various  other 
characteristics  were  studied  and  compared.  Some  of  these  were  toxicity, 
corrosiveness,  chemical  stability,  effect  of  nuclear  radiation,  and  avail¬ 
ability.  As  a  result  of  this  comparison,  no  one  gas  was  found  to  be  ideal; 
however,  the  greatest  favor  fell  to  air  and  to  steam,  with  air  being  given 
first  choice  and  steam  appearing  to  merit  study  should  air  prove  otherwise 
unattractive . 

Concurrently  with  the  study  of  gases,  a  preliminary  design  was  .lade  of 
an  airplane  of  sufficient  scope  to  support  the  study  and  be  appropriate  to 
the  operational  conditions  spelled  out  in  the  contract  with  Wright  Air 
Development  Center,  While  details  of  the  airframe  configuration  and  the 
mission  profiles  are  classified,  it  may  be  stated  that  a  high  speed  air¬ 
craft  was  assumed.  This  airframe  design  was  necessary  in  order  to  logically 
locate  hardware,  to  determine  the  thermal  history  of  the  air  as  it  passed 
through  the  system,  and  to  establish  appropriate  tubing  lengths  for  purposes 
of  weight  study. 

Before  going  on  with  the  history  of  the  study  progreun,  it  nay  be  well 
to  dwell  a  moment  on  one  of  the  characteristics  of  gases  which  has  been 
mainly  responsible  for  the  poor  light  in  v/hich  pneumatics  is  so  often  con¬ 
sidered.  That  characteristic  is  compressibility.  The  direct  measure  of 
compressibility  of  any  fluid  is  bulk  modulus,  which  is  to  fluids  what 
modulus  of  elasticity  is  to  solids.  On  theoretical  grounds,  the  bulk 
modulus  of  a  gas  is  the  product  of  its  specific  heat  ratio  and  the  applied 
pressure. 

If  a  graph  is  made  of  the  bulk  moduli  of  several  fluids  under  2000  and 
4000  pounds  per  square  inch  pressure  at  various  temperatures,  a  number  of 
interesting  circumstances  will  be  noted.  First,  at  ordinary  temperatures 
a  wide  disparity  exists  between  the  bulk  modulus  of  hydraulic  fluid  and 
air  Secondly,  and  more  to  the  jxjint,  it  v;ill  be  seen  that  as  temr^erature 
goes  up,  the  bulk  moduli  of  the  hydraulic  fluids  go  do\m.  In  fact,  at  1000®F, 
if  these  hydraulic  fluids  vrere  sufficiently  stable  chemically  to  operate  at 
that  temperature,  their  bulk  moduli  v/ould  have  fallen  better  than  90  percent 
and  into  the  region  of  air  at  the  same  tempera  !iure. 

On  the  other  hand,  air  over  this  temperature  range  has  a  virtually 
constant  bulk  modulus.  This  then  explains  why  many  hydraulic  servo  systems, 
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othen\rise  stable  at  ordinary  teinoerature  levels,  nay  '-o  unstable  when  the 
temperature  is  increased  In  other  v/ords,  bemperature  change  can  drop  the 
bulk  modulus  out  of  the  stability  lijiiits  inlierent  in  a  given  system  design. 

If  the  system  had  been  designed  for  pneumatics  and  for  the  appropriate  bulk 
modulus  however,  no  change  in  system  stability  would  be  expected.  As  will 
be  pointed  out  later,  the  speed  brake  system  designed  for  this  prograra  ran 
as  stably  at  high  temperature  as  it  did  at  room  temperature  conditions. 

Several  types  of  pneumatic  systems  were  studied  in  relation  to  the 
nature  of  the  power  source  and  power  actuation  components. 

Classified  by  power  source,  the  investigation  covered  engine-mounted 
compressors,  ground  charged  reservoirs,  and  liquefied  gas  converters.  The 
other  pneumatic  power  sources  (using  the  term  in  the  broad  sense)  were  a 
regenerative  steam  system  and  a  solid  propellant  system-  With  regard  to 
power  actuation  components,  the  pneumatic  devices  studied  were  air  turbines 
and  air  motors  as  well  as  pi si on- cylinder  actuators.  For  purposes  of  weight 
comparison,  hydraulic  an^.  oj.ecorical  systems  were  studied,  assuming  materials 
to  exist  which  were  not  actually  available. 

On  the  basis  of  weight  alone,  the  systems  which  proved  the  most  favorable 
were  the  engine-mounted  compressor  and  the  regenerative  steam  system;  however, 
those  weights  were  so  close  to  each  other  that  a  strong  position  in  favor  of 
one  over  the  other  in  this  regard  was  not  warranted. 

System  stability  investigations  indicated  strong  likelihood  that  the 
system  could  be  made  stable.  This,  together  with  high  temperature  material 
studies,  led  to  the  conclusion  that  pneumatics  was  feasible  for  aircraft 
system  operation  in  1000°F  ambients-  To  give  the  greatest  margin  of  stability 
it  was  decided  that  the  pressure  of  the  system  should  go  up  from  the  3000 
pounds  per  square  inch  of  conventional  hydraulic  systems  to  5000  pounds  per 
square  inch  supply  pressure.  Since  no  5000  poxmd  per  square  inch  compressor 
was  available  which  could  operate  at  high  temperatxire,  a  two-pronged  attack 
was  necessary. 

Therefore,  in  early  1957,  Lockheed  proceeded  with  the  Phase  II  program 
for  the  design,  development,  and  test  of  a  5000  pound  i^er  square  inch  system 
assuming  the  successful  developaaent  of  a  compressor  was  forthcoming.  The 
compressor  problem  was  not  only  one  of  designing  for  operation  at  high 
temperature  in  general,  but  included  the  problem  of  precooling,  interstage 
cooling,  and  aftercooling  at  altitudes  v/here  the  heat  transfer  correlations 
were  scant  or  conflicting  or  both.  Subsequently  a  contract  for  the  develop¬ 
ment  of  a  high  temperature  5000  pound  per  square  inch  compressor  was  awarded 
to  the  Walter  Kidde  Company  of  Belleville,  New  Jersey; 

The  pneumatic  speed  brake  system  dovmstream  of  the  compressor  consisted 
generally  of  two  parts.  In  the  so-called  accessory  compartment  are  to  be 
found  the  air  reservoirs,  filters,  relief  valves,  mechanical  moisture  sepa¬ 
rator,  chemical  drier,  pressure  regulator,  and  shut-off  valve  for  the 
necessary  storage  and  handling  of  the  compressor  discharge  air  prior  to 
speed  brake  operation.  The  definition  of  the  system  requirement  established 
the  maximum  ambient  temperattire  in  this  compartment  at  400°F.  The  actuator 
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compartment,  on  the  other  hand,  is  where  the  temperature  rises  to  1000°F. 

In  this  compartment  is  located  the  actuator,  servo  valve,  and  coiled 
tubinj  connecting  the  two,  in  addition  to  rod-end  bearings  and  associated 
linkage  and  fittings  r  Since  this  is  the  area  of  primary  interest  in  this 
paper,  we  will  continue  with  a  discussion  of  these  components . 

As  regards  materials  for  1000®F  usages  a  wide  range  of  materials, 
surface  treatments,  surface  coatings  and  hardnesses  were  considered. 
Approximately  seventy-five  combinations  were  run  off  in  a  wear  testing 
machine  in  an  ambient  temperature  of  1000°F  and  at  bearing  loads  up  to 
6150  pounds  per  square  inch. 

The  materials  tested  included  Inconel-X,  Haynes  Alloy  No.  6B  (both 
wrought  and  cast),  alundum,  and  carbon  graphite  compacts.  Protective 
coatings  included  flame  platings  of  aluminum  oxide  and  tungsten  carbldej 
Electrolizing,  electroless  nickel,  hard  chrome  plate,  nitriding,  chrome 
diffusion,  and  nickel  oxide.  Dry  film  lubricants  included  Electrofilm 
No.  1000  and  No.  2006,  Everlube,  Surf-Kote,  and  Drilube.  Where  a  particular 
combinatioi  showed  promise,  the  test  was  repeated  by  inverting  the  materials 
of  the  wear  block  and  running  race. 

As  borne  out  by  these  tests  and  by  the  subsequent  system  testing,  the 
so-called  flame  platings,  particularly  those  of  tungsten  carbide,  have 
exceptionally  good  wear  characteristics  at  these  temperatures  These 
platings  are  very  hard  coatings  applied  with  high  velocity  to  the  surface 
of  the  part  to  be  plated.  The  base  material  is  in  effect  immersed  in  the 
flame  resulting  from  the  detonation  which  imparts  the  high  velocity  to  the 
plating  material.  The  resulting  physical  and  chemical  bond  effects  a  coat¬ 
ing  which  in  the  "as  used"  condition  is  machined  to  about  .002  of  an  inch 
thickness.  While  the  alviminura  oxide  has  a  lower  coefficient  of  friction  and 
can  stand  higher  temperatures  than  the  timgsten  carbide,  it  is  less  dense 
and  more  brittle.  In  the  various  actuator  assemblies  tested,  failures  of 
the  aluminum  oxide  coating  occurred  from  time  to  time,  wliereas  the  tungsten 
carbide  showed  virtually  no  signs  of  visible  wear  or  scoring  . 

Host  parts  of  the  actuator  are  InconneJ.-X.  The  actuator  cylinders 
were  of  two  types.  One  was  Inconnel-X  Electrolized  on  the  inside  diameter, 
the  other  Inconel-X  with  a  Haynes  Alloy  No.  6B  liner  The  pistons  and 
piston  rods  were  flame  plated  with  either  aluminum  oxide  or  tungsten  carbide 
on  a  base  material  of  Inconel-X.  The  cylinder  wall,  which  is  approximately 
four  inches  on  the  inside  diameter  and  accomodates  an  eleven  inch  piston 
stroke,  has  a  wall  thickness  of  one-half  inch.  Actually,  from  a  pressure 
and  structural  point  of  view,  the  wall  is  twice  as  thick  as  necessary;  but 
to  prevent  breathing  \uider  pressure  the  cylinder  wall  is  one-half  instead 
of  one-quarter  of  an  inch  The  bearing  around  the  piston  rod  is  Haynes 
Alloy  No,  6. 

Piston  seals  tested  have  been  of  two  types,  the  conventional  piston 
rings  and  a  pack  of  segmented  carbon  base  rings ,  The  rod  seals  have  been 
of  four  types,  contracting  piston  rings,  a  pack  of  the  segmented  carbon 
base  rings,  stacked  chevron  packings  made  of  asbestos  contained  in  an 


Inconel  wire  mesh,  and  the  Boeing  ring-spring  seal.  Of  all  of  these,  the 
Boeing  ring-spring  seal  shows  the  most  promise.  \Jhile  it  seems  to  have 
more  friction  than  the  segmented  carbon  ring,  for  instance,  it  is  by  far 
the  best,  having  shown  virtually  no  leakage  on  bench  test. 

The  servo  valves  tested  were  of  two  types,  the  rotary  plate  valve  and 
the  linear  spool  valve-  The  rotary  plate  valve  utilized  a  plate  made  of 
Haynes  Alloy  No.  25  mounted  on  a  flame  plated  base.  Very  high  friction  was 
experienced  with  the  plate  valvej  therefore,  tlie  majority  of  the  system 
testing  was  accomplished  with  the  spool  type- 

To  minimize  the  chance  of  differential  expansion  both  sleeve  and  spool 
were  made  of  the  same  material.  With  respect  to  materials,  both  tool  steel 
and  Haynes  Alloy  No.  25  were  used.  As  it  turned  out,  the  bulk  of  the  cycling 
was  accomplished  with  tool  steel  valves  which  became  blued  with  temperature, 
but  showed  no  excessive  corrosion  and  no  evidence  of  erosion.  The  Haynes 
Alloy  Wo.  25  valve,  while  softer,  is  basically  of  a  high  temperature  material, 
which  would  be  expected  to  exlaibit  no  corrosion. 

Since  the  servo  valve  was  not  mounted  on  the  actuator,  the  two  were 
connected  by  coiled  tubing .  Three-eighth  inch  Inconel-X  tubes  were  woiind 
into  coils  of  four  turns  each  with  a  mean  diameter  of  four  inches  and  allr-"'d 
a  relative  rotation  of  twenty-five  degrees  of  the  actuator  at  its  fixed  pivot. 
These  coils  were  trouble  free  tliroughout  the  program.  Part  of  this  is  due 
to  the  fact  that  the  tubes  underwent  virtually  no  flattening,  when  coiled, 
which  is  directly  attributable  to  the  manufacturing  process.  Before  bending 
the  coils  the  tubes  were  filled  with  Cerrobend,  the  zero  coefficient  of 
expansion  alloy  used  for  special  inspection  techniques  and  other  purposes. 

This  technique  held  flattening  of  the  tubing  to  less  than  one  percent  of  the 
tubing  diameter,  being  on  the  order  of  one  thousandth  of  an  inch  on  an  out¬ 
side  diameter  of  three-eighths  of  an  inch. 

System  testing  was  accomplished  in  two  phases,  stability  and  endurance. 
Stability  testing  was  run  initially  at  room  temperature.  Various  stabiliz¬ 
ing  techniques  were  employed,  such  as  adding  velocity  feedback  to  the 
mechanical  displacement  feedback  signal.  Provisions  were  also  made  for 
varying  the  gain  of  the  feedback  signal  Initally,  the  system  was  run  under 
a  torsion  bar  load  to  simulate  the  aerod-^uiamic  forces.  Each  torsion  bar 
was  four  inches  square  in  section  and  had  an  effective  length  of  about  43 
feet.  At  full  speed  brake  actuation  this  represents  approximately  60,000 
foot-pounds  of  torque,  with  an  actuator  force  of  33,000  pounds.  Later  the 
load  was  removed.  In  all  cases  the  system  operation  was  stable,  not  only 
at  room  temperature,  but  also  at  the  maximum  ambient  temperature  v;hich  was 
in  excess  of  1050®F. 

The  system  endurance  testing  was  largely  accomplished  at  elevated 
temperature.  Except  for  stability  testing,  system  check-out  luns,  and  some 
cycling  at  intermediate  points,  endurajice  running  was  accomplished  in  a  1000°F 
ambient.  In  general,  the  endurance  testing  verified  the  results  of  the  wear 
tests  run  off  early  in  the  program-  The  tungsten  carbide  flame  plating  stood 
up  very  well,  as  did  the  Haynes  Alloy  No  6B  and  the  Electrolized  Inconel-X, 


-5- 


It  is  believed  that  the  aluminum  oxide  flame  plating  would  have  also  made 
a  good  showing  if  failure  of  it  had  not  been  aggravated  by  a  seizure  of  the 
rod-end  bearing.  ’Then  this  occurred,  the  ailuminum  oxide  failed  locally. 
Then,  by  virture  of  its  highly  abrasive  nature,  mating  surfaces  became 
heavily  scored,  causing  further  failure  of  the  aluminum  oxide.  On  the 
other  hand,  a  rod-end  bearing  failure  had  little  effect  on  the  tungsten 
carbide  flame  plate. 

A  case  in  point  to  the  credit  of  the  tungsten  carbide  flame  plate 
was  demonstrated  when  a  rod  packing  backup  ring  shrunk  more  than  the  design 
allowance  in  operation.  The  ring  was  bound  so  tightly  on  the  piston  rod 
that  it  could  only  be  removed  by  cutting  it  off.  As  would  be  expected,  the 
ring  was  severely  scored;  however,  the  flame  plate  on  the  rod  showed  practi¬ 
cally  no  effect. 

In  conclusion,  it  may  be  stated  that  pneumatic  systems  can  be  made  to 
operate  stably.  Each  system  must  be  carefully  designed  recognizing  the 
limitations  which  air  imposes  when  used  as  the  service  fluid.  With  some 
allowance  for  friction,  the  way  to  successful  sealing  of  the  piston  rod  is 
becoming  manifest.  Sealing  of  the  piston  remains  as  one  of  the  problem 
areas,  but  application  of  the  rod  seal  technique  to  the  piston  may  have 
considerable  merit.  Until  the  sealing  problem  has  been  ironed  out,  broad 
conclusions  are  yet  forthcoming,  but  it  has  been  demonstrated  that  pneumatic 
systems  can  be  made  stable. 

liaterials  and  surface  finishes  are  available  which  indicate  that 
operation  in  ambients  as  high  as  1000°F  are  possible  for  aircraft  usage. 
Generally,  the  metals  tested  are  fairly  high  on  the  stategic  materials 
list  or  often  involve  processes  which  border  on  the  exotic.  A  logical  next 
step,  of  course,  is  to  go  to  more  commonly  used  materials.  In  any  event, 
the  program  shows  high  temperatvire  operation  to  be  entirely  possible. 

In  summary,  it  is  felt  that  this  program  has  not  only  demonstrated  the 
feasibility  of  pneumatics  but  has  pointed  the  way  to  a  successful  solution 
of  the  radiation  problem  in  aircraft  actuation  systems. 
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AIRCRAFT  RADOME  DESIGN  PROBLEMS 
ASSOCIATED  WITH  A  NUCLEAR  ENVIRONMENT 


by 

Frank  W.  Thomas 

Lockheed  Nuclear  Products 
Lockheed  Aircraft  Corporation 
Georgia  Division,  Marietta,  Georgia 

The  effects  of  nuclear  radiations  In  relation  to  microwave  transmission 
through  a  dielectric  lens,  or  radome,  presents  a  number  of  design  prob¬ 
lems  associated  with  nuclear-powered  aircraft.  A  mathematical  pre¬ 
diction  technique  was  used  to  estimate  the  effect  of  radiation  on  multi¬ 
layer  dielectric  flat  panels. 

The  results  show  how  radiation  effects  data  may  be  interpreted  in  the 
design  of  a  radome  test  panel. 

RADOMES  IN  THE  ANP  PROGRAM 

The  advent  of  nuclear  power  has  introduced  many  special  technical  problems  into 
aircraft  radome  design.  Some  of  these  problems  concern  only  effects  of  radiation  on 
matter;  but  others  concern  also  additional  aircraft  environments,  such  as  temperature. 
The  design  of  radome  dielectric  lens  for  a  nuclear  aircraft  Is  considerably  more  complex 
than  would  be  the  case  with  chemically  powered  aircraft.  This  paper  Is  concerned 
with  the  analysis  and  the  effects  of  nuclear  radiation  on  radome  wall  materials  and 
especially  with  the  transmissivity  of  a  radome  wall  composed  of  these  materials. 

Electrical  parameters  of  a  radome  wall  are  changed  as  a  result  of  nuclear  radlotlon. 
Extensive  changes  can  limit  radome  performance;  therefore  the  scope  of  such  effects 
must  be  established.  Simple  changes  In  electrical  properties  can  best  be  Interpreted 
os  total  effects  on  transmission  through  a  particular  radome  wall,  so  some  means  must  be 
employed  to  Interpret  radiation  Induced  electrical  changes  In  materials  In  terms  of 
actual  radome  wall  performance.  The  primary  object  of  this  paper  Is  to  show  the 
application  of  a  predictability  technique  In  analyzing  typical  radome  wall  flat  panel 
sections  that  are  exposed  to  radiation.  The  proper  use  of  the  predictability  technique 
can  reduce  drastically  the  environmental  testing  requirements  that  would  otherwise 
be  needed  to  test  full-scale  microwave  lenses  In  a  nuclear  radi  itlon  environment. 
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Following  are  listed  radome  properties  that  are  affected  by  radiation  and  related 
aircraft  environmental  factors: 

Temperature  coefficient  of  dielectric  constant  and  loss  tangent 

Electrical  thickness  equivalent  skin  or  core 

Insertion  phase  difference 

Boresight  error 

Total  transmissivity 

Radloactivatlon  of  materials 

Phase  effects  of  Incident  and  reflected  waves 

Changes  In  electrical  constants  of  dielectric  arrays 

Antenna  near  field  patterns 

Far  field  energy  distribution 

Dielectric  anisotrophy  of  radome  wall  materials 

Beam  shift 

Polarization  of  transmitted  waves 
incidence  angle  of  beam 
Relaxation  effects 
Anisotrophy 
Reflection  coefficient 
Insertion  phase  shift 

RADIATION  AND  RADOME  DESIGN 

It  Is  first  necessary  to  describe  briefly  some  of  the  more  fundamental  aspects  of  the 
problem  of  correlating  radiation  effects  with  the  design  of  a  dielectric  lens.  Radomes 
are  conventionally  used  for  gun  laying,  Are  control,  bombing,  and  navigation.  They 
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njust  meet  stringent  mechanical  and  structural  specifications,  which  may  also  be  affected 
by  nuclear  radiations.  High  strength,  electrical  stability,  and  physical  serviceability 
over  a  wide  range  of  temperatures  and  microwave  frequencies  are  necessary.  Furthermore, 
the  shape  of  a  radome  must  conform  to  aerodynamic  and  microwave  optical  design  criteria. 
For  these  reasons,  a  tool  for  the  Interpretation  of  radiation  effects  Is  needed  In  the  pre¬ 
liminary  design  of  any  nuclear  aircraft  microwave  dielectric  lens  device. 

The  effects  of  the  various  kinds  of  nuclear  radiation  are  briefly  considered  In  selecting 
materials  for  a  radome  wall.  Numerous  tests  have  been  conducted  on  unirradlated 
materials,  and  Wright  Air  Development  Center  has  conducted  an  extensive  study  of  effects 
of  gamma  radiation  on  radome  materials.  However,  relatively  little  work  has  been  done 
on  pertinent  electrical  parameters  of  radome  material  during  and  after  neutron  Irradiation. 
Wave  guide,  or  cavity,  measurements  made  on  materials  being  Irradiated  In  a  nuclear 
environment  at  microwave  frequencies  are  not  folly  adequate  to  assess  actual  radome  wall 
performance. 

Tests  made  before  and  after  Irradiation  do  not  show  the  transient  effects  on  transmission 
that  may  occur  in  a  radome  wall.  Little,  therefore,  has  been  determined  on  what  happens 
during  nuclear  Irradiation  of  an  entire  radome  module.  Actually,  transient  effects  may  be 
of  particular  Importance,  since  Irradiation  and  transmission  must  take  place  simultaneously. 

RADOME  RADIATION  ENVIRONMENT 

Various  processes  for  the  interaction  of  particle  and  electromagnetic  radiation  In 
radome  wall  materials  are  covered  only  briefly  in  this  paper. 

Interactions  produce  two  main  effects:  absorption  of  the  photon  and  scattering  or 
reflection  from  either  the  atomic  nucleus  or  atomic  electron  cloud.  Of  the  atoms  within 
the  dielectric,  at  least  a  dozen  combinations  of  Interactions  and  effects  are  possible; 
but  for  dielectric  lens  design  purposes  only  the  following  need  be  considered: 

Compton  Effect:  The  electrical  properties  of  a  radome  dielectric  are  changed 
by  the  Compton  process.  This  is  caused  by  ionization  resulting  from  Compton 
scatter  electrons  and  their  interaction  with  other  atoms. 

photoelectric  Effect:  The  photoelectric  and  Compton  processes  are  competitive, 
the  latter  is  perhaps  the  more  important  In  Its  effect  on  dielectrics  for  the 
lower  photon  energy  ranges  encountered. 

Pair  Production:  The  effect  of  the  pair  production  process  Is  to  augment  the 
Compton  process  damage  in  radome  dielectric  materials,  since  the  0.51-mev 
photons  from  positron  annihilation  are  largely  absorbed  In  radome  materials. 

The  pair  process  Is  not  likely  to  affect  radome  materials  unless  photon  energies 
are  adequate. 
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Neutron  Radiation:  Many  radome  materials  contain  hydrogen  and  nitrogen.  A 
reaction  with  hydrogen  produces  o  gamma  photon.  Reaction  with  nitrogen  can 
yield  a  proton  and  carbon-14.  Neutron  alpha  reactions  cause  short  but  Intense 
Ionization  tracks  In  the  dielectric.  Scattering  of  fast  neutrons  can  produce 
nucleus  and  proton  recoil. 

RADIATION  ON  RADOME  MATERIALS 

Greatest  concern  In  the  selection  of  materials  suitable  for  radome  walls  Is  with  the 
photoelectric,  Compton,  and  neutron  Interactions. 

Organic  radome  materials  such  as  resin  fiberglas  laminate  skins;  foam  plastic  cores; 
and  resin  fabric  or  fiberglas  honeycombs  composed  of  carbon,  hydrogen,  oxygen,  and 
nitrogen  may  be  particularly  affected. 

Recoil  protons  derived  from  fast  neutron  collisions  In  organic  materials  lose  their 
energies  by  Interaction  with  orbital  electrons,  resulting  In  Ionized  atoms  or  molecular 
fragments.  This  can  effect  dielectric  changes  and  result  In  molecular  fragment 
recombinations,  yielding  new  molecular  species  with  resultant  change  In  electrical 
properties  of  the  radome  irraterials. 

Final  degradation  of  thermal  energy  of  a  neutron  and  Its  resultant  capture  by  a 
hydrogen  atom  will  produce  a  gamma  ray  of  2.11  mev,  thus  Increasing  the  original 
gamma  ray  flux  reactions.  The  use  of  lithium,  boron,  or  excessive  nitrogen  In  radome 
materials  may  be  undesirable  because  of  thermal  neutron  capture  and  resultant 
Intense  ionization  caused  by  the  secondary  radiations.  Elements  that  yl  eld  high 
energy  beta  radiation  may  also  be  undesirable  In  radome  design.  Complex  dielectric 
properties  may  be  especially  sensitive  to  flux  rate. 

DATA  ANALYSIS  BY  MEANS  OF  THE  PREDICTABILITY  TECHNIQUE 

Transient  and  post-Irradlatlon  test  data  on  basic  radome  materials  alone  do  not 
show  the  ultimate  effect  of  nuclear  radiation  on  radome  performance.  The  net 
effect  of  changes  In  loss  tangent  and  dielectric  constant  must  be  stated  in  terms  of 
power  transmission  for  a  specific  wavelength  through  a  specific  wall.  Such  a  wall 
may  be  symmetrical  or  unsymmetrtcal  and  composed  of  a  mono  layer  or  multi -dielectric 
layers.  Three-layer  unsymmetrical  walls  were  selected  as  possible  examples  for 
analysis  of  radiation  effects. 

The  complex  parameter  changes  resulting  from  radiation  were  related  and 
presented  in  final  form  as  power  transmission  through  the  dielectric  lens. 

For  the  analysis  of  typical  radome  wall  flat  panel  sections,  the  use  of  a  predicta¬ 
bility  computer  program,  developed  by  William  Schroder,  Lockheed  Aircraft  Corporation, 
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Calffornia  Division  ,  Is  employed. 


An  IBM  704  compufer,  which  prints  transmission  coefficients  In  tabular  form,  can  be 
readily  programmed  to  accept  as  Input  data  the  loss  tangent,  dielectric  constant  of 
each  sheet,  and  thickness  and  number  of  dielectric  layers. 

Output  data  will  then  be  the  power  transmission  coefficient  as  a  function  of  the 
Indicent  angle  and  polarization  of  the  beam.  Details  of  the  actual  IBM  routine  and 
the  equations  used  to  establish  the  IBM  program  employ  Maxwell's  relations  and 
related  expressions  for  energy  absorption,  reflectance,  and  transmission  In  terms  of 
complex  electrical  parameters.  Manual  computation  Is  not  practical  because  of  the 
time  required.  Theoretical  estimates  of  the  results  of  nuclear  radiation  on  radome 
dielectric  sandwich  panels  were  made  by  the  use  of  this  program.  Changes  of  +10% 
and  -10%  In  dielectric  constant  and  loss  tangent  of  the  skin  and  core  of  the  panel 
were  considered  equivalent  to  those  Induced  by  radiation,  and  power  transmission 
was  calculated  for  three  unsymmetrical  radome  flat  well  panels. 

The  theoretical  analysis  of  three  typical  panels  Is  presented  In  Figures  1  through 
8.  Each  flat  panel  consists  of  a  core  and  two  skins  or  faces.  The  electrical  and 
physical  properties  of  the  core  and  each  skin  are  defined  In  these  figures. 

Figures  1,  2,  and  3  show  that  power  transmission  for  the  test  panels  Is  changed 
about  plus  or  minus  3%  when  a  change  of  Initial  parameters  Is  +  10%.  This  may  not  be 
too  significant  for  systems  where  normal  efficiencies  are  high,  but  It  could  be  critical 
for  those  with  marginal  design  level  as  regards  range  and  sensitivity  and  energy  reflection. 

Figures  4,  5,  and  6  show  that  transmission  Is  less  affected  by  induced  changes  for 
parallel  polarization. 

Figures  3  and  6  show  that  a  tapered  panel  optimized  as  to  core  thickness  and 
Incidence  angle  provides  Increased  transmissivity  and  that  changes  that  might 
result  from  radiation  would  be  less  significant. 

Figures  7  and  8  show  that  polarization  influences  transmission  more  at  high 
Incidence  angles. 

Although  emphasis  was  placed  on  assumed  changes  In  electrical  properties  and 
the  span  of  these  changes  was  arbitrarily  chosen,  the  predictions  Indicate  the 
ultimate  results  of  radiation  on  a  dielectric  lens  assembly  when  actual  changes  are 
used  Instead  of  simulated  ones. 

Advanced  mathematical  routines  for  the  IBM  704  are  being  devised  to  analyze 
radomes  of  complex  shape  or  contour.  Such  programs  will  permit  rapid  analysis  of 
proposed  radome  designs,  allowing  faster  correlation  of  aerodynamic,  structural, 
nuclear,  and  microwave  design  requirements. 
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Ultimate  prediction  capability  of  the  program  for  any  shape  radome  as  a  function 

of  radiation  would  Include  the  following: 

Beam  shift 
Phase  front  shift 
Transmission 
Bore  sight  error 
Far  field  pattern 

CONCLUSIONS 

I*  Radome  parametric  studies  Involve  such  a  number  of  variables  that  direct  measure¬ 
ment  of  all  parameters  simultaneously  In  a  nuclear  environment  for  an  actual  radome 
Is  difficult. 

2.  Mathematical  routines  using  known  perturbation  of  dielectric  wall  properties  as  a 
function  of  nuclear  radiation  Integrated  dose  can  predict  radome  performance. 

3.  Estimates  using  these  techniques  Indicate  that  the  overall  effect  of  radiation  may 
not  be  too  great  on  actual  radome  performance  for  a  ±  10%  variation  of  electrical 
properties  within  the  limits  of  the  design  cases  studied. 

4.  Correlation  of  radiation  effects  In  dielectric  materials  with  actual  performance  In  a 
dielectric  lens  can  best  be  made  using  a  computer  program. 

5.  Correlation  of  dose  and  change  In  transmission  to  any  design  panel  can  be  performed. 
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FIG.  1  -  CONVENTIONAL  UNSYMMETRICAL  WALL  ELECTRIC  PARAMETERS 
(PERPENDICULAR  POLARIZATION) 


RG.  2  -  CONVENTIONAL  UNSYMMETRICAL  WALL  CORE  OPTIMIZED  FOR  MAXIMUM  LOW 
INCIDENCE  ANGLE  TRANSMISSION  (PERPENDICULAR  POLARIZATION) 
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RAY  INCIDENCE  ANGLE  IN  DEGREES 

FIG.  4  -  CONVENTIONAL  UNSYMMETRICAL  WALL  ELECTRIC  PARAMETERS  (PARALLEL  POLARIZATIONS) 


RAY  INCIDENCE  ANGLE  (9)  IN  DEGREES 
RG.  5  CONVENTIONAL  UNSYMMETRICAL  WALL  CORE  OPTIMIZED  FOR  MAXIMUM 
LOW  INCIDENCE  ANGLE  TRANSMISSION  (PARALLEL  POLARIZATION) 


RAY  INCIDENCE  ANGLE  (e)  IN  DEGREES 
RG.  6  UNSYMMETRICAL  TAPERED  PANEL  WALL  OPTIMIZED  FOR  EACH 
INCIDENCE  ANGLE  (PARALLEL  POLARIZATION) 
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FIG.  7  CONVENTIONAL  UNSYMMETRICAL  WALL  ELECTRIC  PARAMETERS 
(PERPENDICULAR  &  PARALLEL  POLARIZATION) 


a. 


RAY  INCIDENCE  ANGLE  0  IN  DEGREES 

FIG.  8  -  CONVENTIONAL  UNSYMMETRICAL  WALL  CORE  OPTIMIZED  FOR  MAXIMUM  LOW 

INCIDENCE  ANGLE  TRANSMISSION  (PERPENDICULAR  &  PARALLEL  POLARIZATION) 


FIG.  9  -  UNSYMMETRICAL  TAPERED  PANEL  WALL  OPTIMIZED  FOR  EACH  INCIDENCE  ANGLE 
(PERPENDICULAR  AND  PARALLEL  POLARIZATION) 
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ON  THE  ENERGY  LEVELS 
IN 

NEUTRON -IRRADIATED  P-TYPE  SILICON 


by 

C.  A.  Klein  and  W.  D.  Straub 
Research  Division,  Raytheon  Manufacturing  Company 
Waltham,  Massachusetts 


High  resistivity  p-type  silicon  samples  were  irradiated 
in  the  Brookhaven  pile  for  short  periods  of  time  resulting  in 
integrated  fast  neutron  fluxes  of  the  order  of  a  few  10  n/cm. 
Hall  coefficient  and  conductivity  measurements  were  then 
performed  over  extended  temperature  ranges.  A  detailed 
analytical  investigation  of  these  data  provides  evidence  for 
at  least  two  bombardment -induced  hole  trap  levels  in  the 
lower  half  of  the  energy  gap.  The  deep  trap,  at  0.29  ev  from 
the  valence  band,  is  introduced  at  a  rate  which  appears  to  be 
proportional  to  the  integrated  flux.  Moreover,  it  has  been 
established  that  for  this  level  the  product  of  the  statistical 
weight  factor  by  the  temperature  shift  factor,  as  defined 
below,  is  practically  equal  to  one.  The  shallow  trap,  at 
0.16  ev  from  the  valence  band,  is  not  yet  as  fully  describable. 
The  work  is  still  in  progress,  and  only  the  established  results 
at  the  time  of  this  writing  (August  1958)  are  given  here.  The 
whole  approach  is  of  a  "phenomenological"  nature,  m  the  sense 
that  no  attempts  were  made  to  correlate  the  presently  available 
data  with  possible  types  of  defects  or  defect  models. 
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I.  INTRODUCTION 


The  electrical  and  optical  properties  of  semiconducting  materials 
are  known  to  be  extremely  sensitive  to  any  high-energy  nuclear  radia¬ 
tion.  It  has  been  recognized  that  this  sensitivity  is  primarily  due  to 
the  introduction  of  new  energy  states  in  the  forbidden  band,  as  a  result 
of  the  disorder  induced  in  the  lattice  structure  during  the  bombardment. 
Considerable  effort  has  been  spent  on  the  determination  of  the  position 

and  characteristics  of  the  "defect  states"  produced  in  germanium  by  a 
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wide  variety  of  energetic  particles  ’  ’  .  Much  less  extensive  work 
has  been  done  on  the  radiation  effects  in  high-purity  silicon,  mainly 
because  single -crystal  specimens  have  only  recently  become  generally 
available,  and  also  because  it  was  felt  that  the  position  of  the  defect 
states  in  the  forbidden  band  of  silicon  is  such  that  it  is  more  difficult 

4 

to  locate  them  . 

The  present  paper  is  concerned  with  the  nature  of  defect  states 

introduced  in  the  lower  half  of  the  energy  gap  of  silicon  by  pile -neutron 

irradiation.  The  samples  studied  were  cut  from  a  high -resistivity 

undoped  single  crystal  (K-121)  and  of  the  standard  "bridge"  form. 

Hall  coefficient  measurements  performed  before  irradiation  showed 

that  they  were  all  p-type  with  a  net  impurity  concentration  of  about 
14  -3 

2  X  10  cm  with  boron  as  majority  impurity.  These  samples  were 
packed  in  dry  ice  and  irradiated  in  the  BNL  graphite  reactor  for  short 
periods  of  time  (a  few  minutes  or  less),  using  one  of  the  pneumatic 
tube  facilities^.  Hall  coefficient  and  conductivity  have  then  been  measured 
as  a  function  of  temperature  after  room -temperature  storage  periods 
extending  from  several  days  to  more  than  a  year.  The  observed 

^The  exposure  times  were  determined  by  L.  Davis,  Jr  .  of  Raytheon's 
Research  Division,  on  the  basis  of  results  obtained  at  the  Oak  Ridge 
National  Laboratory  indicating  that  under  similar  conditions  the  number 
of  defects  per  fast  incident  neutron  is  approximately  1. 6  in  germanium, 
and  chosen  so  that  the  concentration  of  bombardment -induced  defects 
should  be  comparable  to  the  net  chemical  impurity  density. 
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effects,  therefore,  should  be  practically  identical  with  those  produced 
by  a  fast -neutron  bombardment  at  or  near  room  temperature. 

The  purpose  of  this  paper  is  to  present  information  on  the  defect 
states  in  irradiated  silicon.  This  information  could  be  gained  by  a 
detailed  quantitative  analysis  of  some  of  our  data.  It  was  not  attempted 
to  correlate  our  results  with  existing  ideas  or  models  concerning  the 
properties  of  the  energy  states  associated  with  Frenkel  defects  or  any 
defect  expected  to  result  from  irradiation.  Theoretical  fits  of  the  charge 
carrier  concentration  in  bombarded  samples,  from  about  100®K  up  to 
450*K,  based  on  well-established  laws  and  obvious  assumptions,  per¬ 
mitted  us  to  locate  and  partly  describe  two  of  the  defect  states  in  the 
lower  half  of  the  energy  gap  of  neutron -bombarded  p-type  silicon. 

II.  EXPERIMENTAL  RESULTS 


For  two  of  the  irradiated  samples  {6a  and  8b),  the  charge  carrier 
concentration  as  a  function  of  the  reciprocal  temperature  is  given  in 
Fig.  1.  The  same  figure  also  exhibits  typical  carrier  concentration 
data  for  an  unbombarded  sample  cut  from  crystal  K-121.  These  concen 
trations  were  obtained  from  the  Hall  coefficient  measurements  and  the 
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ratio  of  Hall  mobility  to  drift  mobility  according  to  the  formula 


P  = 


1 

Re 


(1) 


The  value  of  the  mobility  ratio  at  a  given  temperature  above  100“ K  was 

estimated  on  the  basis  of  the  results  of  Morin  and  Maita  for  high -purity 
6  b 

p-type  silicon  ’  It  can  be  seen  from  the  carrier  concentration  data  of 


Recent  experiments  by  Long  seem  to  indicate  that  the  linear  negative 
temperature  dependence  of  Hall  to  conductivity  mobility  and  the  result 
that  than  unity  at  300* K  are  truly  intrinsic  properties 

of  p-t^e  silicon,  which  may  be  due  to  the  warping  of  the  constant - 
energy  surfaces.  ^ 
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the  unbombarded  sample  that  such  a  mobility  ratio  correction  leads  to 

the  expected  absence  of  any  temperature  dependence  in  the  hole  density 

above  IBO'K,  where  all  the  chemical  impurities  are  ionized.  On  the 

other  hand,  the  ratio  is  believed  to  change  very  little  with  bombardment, 
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as  long  as  impurity  scattering  is  not  too  important  ’  .  It  was  felt, 
therefore,  that  above  100  or  150“K,  the  hole  concentration  in  bombarded 
samples  should  be  as  given  by  Eq.  (1)  with  the  pre -bombardment  mobility 
ratio  values.  The  marked  temperature  dependence  of  the  charge  carrier 
concentration  in  irradiated  samples,  especially  the  very  large  removal 
observed  at  low  temperatures,  implies  the  existence  of  bombardment - 
induced  hole  trapping  centers  which  manifest  themselves  by  a  rapid 
’’freezing  out”  of  free  carriers,  as  the  temperature  is  lowered.  A 
quantitative  analysis  of  this  trapping  process  was  the  main  object  of 
our  investigation,  and  will  be  presented  in  part  III  of  this  report. 

As  mentioned  in  the  introduction,  the  irradiations  were  performed 
in  a  pneumatic  tube  facility  (PN-6)  of  the  Brookhaven  reactor.  Sample 
6a  was  bombarded  for  two  minutes  in  May  1957,  and  the  "damaging 
dosage”  was  estimated  of  the  order  of  1. 2  x  10  n/  cm  .  In  view  of 
the  difficulties  in  determining  the  effective  dose  and  the  general  con¬ 
fusion  which  exists  in  this  matter,  it  might  be  of  interest  to  report 
how  this  number  was  obtained.  According  to  Primak,  the  fast -neutron 
damage  produced  in  a  nuclear  reactor  has  to  be  related  to  the  uppermost 
two  decades  of  energy  below  the  energy  at  which  the  pile  neutron  spectrum 
begins  to  fall  greatly  from  a  "resonance"  distribution^*^.  The  damaging 
dosage  (nvt)^  is  therefore  defined  as  the  time  integral  of  the  resonance 
flux  over  a  two -decade  range  of  energies  at  the  upper  end  of  the  abundant 
part  of  the  pile  neutron  spectrum.  Assuming  a  cadmium-ratio  of  33, 
the  resonance  flux  in  any  energy  decade  is  known  to  be  just  one  seventh 

of  the  thermal  neutron  flux^^.  In  May  1957,  the  thermal  flux  in  PN-6 

12  -2  -1 

was  given  as  3.  5  x  10  cm  sec  .  The  resonance  flux  in  this  facility 

This  value  was  generally  quoted  for  the  BNL -reactor  lattice  before  the 

recent  fuel  element  conversion^  1 . 
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was  then  5x10^^  cm  sec  per  energy  decade,  so  that  a  two -minute 

14 

exposure  probably  resulted  in  a  damaging  dosage  of  about  1.2x10 
n/ cm^.  Sample  8b  was  irradiated  in  November  1957  under  similar 
conditions,  but  for  four  minutes.  At  that  time  the  fuel  element 
conversion  process  was  already  started,  and  accurate  neutron -flux 
data  were  not  available.  Nevertheless,  the  Brookhaven  staff  felt 
that  the  resonance  flux  in  PN-6  should  not  have  significantly  changed, 
so  that  the  damaging  dosage  received  by  sample  8b  was  tentatively 
assumed  to  be  twice  the  damaging  dosage  of  sample  6a,  i.  e.  , 

2. 4x  10^'^  n/ cm^. 

In  connection  with  sample  8b,  it  will  be  noted  (Fig.  1)  that  two 

temperature  runs  were  performed,  the  first  in  December  1957,  about 

a  month  after  irradiation,  and  the  second  in  April  1958.  At  first 

glance  it  might  appear  that  the  post -bombardment  carrier  concentration 

was  practically  not  affected  by  five  months  of  room -temperature  aging, 

and  this  seems  to  confirm  the  exceptional  stability  of  radiation -induced 

12 

defects  in  silicon  already  noticed  in  electron -bombarded  specimens 
Closer  examination,  however,  reveals  that  at  low  temperatures  the 
hole  concentration  slightly  recovered,  possibly  as  a  result  of  an  energy 
level  shift  in  the  bombardment -induced  traps.  This  point  will  be  examined 
in  detail  later  on. 

Attempts  to  fit  the  charge  carrier  concentration  data  made  it  clear 
that  a  reliable  analysis  had  to  be  based  on  a  thorough  understanding  of 
the  Fermi -level  behavior  in  the  studied  specimens,  before  and  after 
irradiation.  This  behavior  is  illustrated  in  Fig.  2,  which  gives  the 
Fermi -energy  as  a  function  of  temperature  for  the  three  investigated 
samples.  The  Fermi-energies  were  computed  from  the  carrier 
concentration  data  using  the  standard  expression  for  the  hole  density 
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in  a  non -degenerate  semiconductor  , 

Zirm  kT  -  E. 

p  =  — )  exp  ( — 

h 

where  E^  -  E^  is  the  energy  interval  between  the  Fermi -level  and  the 
edge  of  the  valence  band,  and  m^  the  density-of-states  mass  for  the 
valence  band.  This  density-of-states  mass  is  given  by  the  following 
relationship: 

where  the  subscript  L  refers  to  the  light  holes  and  the  subscript  H 
13 

to  the  heavy  holes  .  With  m  .  =  0.14m  and  m  „  =  0.  53  m  , 

pL  o  pH  o 

as  calculated  for  silicon  by  Lax  and  Mavroides  from  cyclotron  resonance 
data  at  4.  2*K,  one  obtains 

3/2  3/2 

)  =  5.50x10^® 

an  approximation  we  use  throughout  this  paper.  The  resulting  Fermi- 
energies,  if  plotted  versus  temperature  on  a  linear  scale,  as  in  Fig.  2, 
strongly  indicate  that  the  main  effect  of  radiation  damage  in  silicon  is 
the  introduction  of  new  energy  levels  in  the  forbidden  gap.  These  new 
energy  levels  manifest  themselves  either  by  anchoring  the  Fermi -level 
high  above  the  valence  band,  as  in  sample  8b,  or  by  an  appreciable 
slowing  down  of  the  level's  movement,  as  in  sample  6a. 

II.  THEORETICAL  ANALYSIS 
A.  SAMPLE  8  B 

The  carefully  measured  carrier  concentrations  in  sample  8b,  as 
yielded  by  the  second  temperature  run  (Fig.  3),  were  the  first  to  be 
submitted  to  an  analytical  treatment.  More  specifically,  it  was  attempted 


2irm  kT 
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to  interpret  these  data  in  terms  of  a  single  discrete  bombardment - 
induced  hole  trap  level  possibly  affected  by  a  temperature  dependence. 
According  to  this  assumption,  and  in  the  investigated  temperature 
range  (150 -400 ®K),  the  charge  carrier  concentration  should  be  given  by 

p  =  Nc  -  N^.  (T)  .  (5) 

14 

where  N  is  the  net  chemical  impurity  density  (N  =  N  -  N,  =  2.02x  10 

^  ^  Si  Cl 

cm  ,  as  indicated  by  pre -irradiation  measurements),  and  N^,  the  density 
of  the  filled  hole  traps.  The  law  of  mass  action  shows  that 

^  “T  +  K(T)/p  ’ 

where  is  the  total  trap  concentration,  and  K(T)  the  equilibrium  constant 
for  this  trapping  process.  If  the  free  carriers  are  non -degenerate,  one 
can  write 


K(T) 


(7) 


the  Q's  representing  the  respective  partition  functions  of  free  holes,  empty 
and  filled  traps.  Therefore, 


2  0.,  2T7m  kT  E. 

K(T)  =  (  --I  )  exp  i--^)  .  (8) 

if  is  the  ratio  of  the  statistical  weight  of  the  empty  trap  to  that 

of  the  full  one,  and  E^  the  energy  level  of  the  trap  measured  from  the  top 
of  the  valence  band.  It  is  usual  to  assume  a  linear  temperature  variation 
for  such  an  energy  level: 


E^  (T)  =  E^(0)  +  aT 


(9) 
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If  this  is  correct  in  the  case  of  a  bombardment -induced  level,  we 
find  that  the  equilibrium  constant  can  be  written  as 


K(T)  =  rC(T)exp(-^)  ,  (10) 

where  F  is  the  product  of  the  statistical  weight  factor  y  =  2 

by  the  temperature -shift  factor  t  =  exp  (q/k)  of  the  level,  a  level  whose 

zero -temperature  position  is  at  an  energy  from  the  valence  band. 

C(T)  =  5.50  X  10^®  (T/300)^/^  cm  ^  as  previously  calculated  for  holes 
in  silicon  (Eq.  4). 


After  bombardment,  and  at  temperatures  below  250“K,  the  charge 
carrier  concentration  in  sample  8b  is  seen  (Fig.  3)  to  be  exceedingly 
small,  if  compared  to  the  initial  concentraion  N  .  We  may,  therefore, 
write 


^t 

^t*  ■  1  +  K(T)/p 


(11) 


for  1000/T>  4.5,  or 


N 


N.  -  N 
t  c 


K(T) 


Substitution  of  Eq.  (10)  in  Eq.  (12)  yields 


p(T/300) 


■3/2 


5.50x10 

77 


18  .3 
-  exp  ( -  )  cm 


(12) 


(13) 


where 


7?  = 


VNq  -  1 

r 


(14) 


is  a  dimensionless  quantity  critically  dependent  on  the  ratio  of  the 
radiation -induced  trap  density  to  that  of  the  net  chemical  impurities. 
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The  experimental  values  of  log  [  p  {T/300)"^^^]  at  low  temperatures  are 
plotted  versus  the  reciprocal  temperature  in  Fig.  3;  they  exhibit  a 
characteristic  slope  corresponding  to  an  activation  energy  of  0.287  ev. 

On  the  basis  of  Eq.  (13),  this  activation  energy  may  be  considered  as 
the  zero -temperature  energy  level  of  a  bombardment -induced  hole 
trap  in  our  p-type  silicon  sample  8b  at  the  time  the  measurements 
were  performed  (about  six  months  after  irradiation) .  The  associated 
value  of  the  parameter  7?  is  then  determined  as  equal  to  0. 350  ±0.005. 

A  single  level  approach  has  to  be  substantiated  by  proving  that, 
with  the  parameters  obtained  from  the  low -temperature  slope,  the  model 
provides  correct  carrier  concentrations  in  the  higher  temperature  range. 
According  to  Eqs.  (5)  and  (6),  the  hole  density  in  an  irradiated  sample  is 
given  by 


N. 


P  =  N, 


T  +  K(T)/p 


a  quadratic  equation  in  p,  the  solution  of  which  is 


(15) 


p  =  -  i  IK  (T)  +  -  N^]  +  J^  [K{T)  +  +  N^K(T) 

(16a) 

or 

p  =  -  ^  LK(T)  +  r  17NJ  +  y^LK:(T)  +  r77N^]^  +N^K(T) 

(16b) 

if  we  introduce  the  parameter  77  .  Let  us  define  the  function  K  (T) 
as  follows 


IC(T) 


C(T)  exp  (- -j^) 


(17) 
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Then  p  is  given  by 


p  =  { -  L  K  (T)  +  77  ]  +  /ik:{t>477]^  +  }  , 

(18) 

and  this  expression  may  be  used  to  test  the  single -level  model  at  high 
temperatures,  since  the  only  adjustable  parameter  is  F  .  Results  of 
some  of  the  computations  are  illustrated  in  Fig.  3.  They  strongly  suggest 
that  our  8b  data  (run  2)  are  correctly  described  by  a  single  discrete  hole 
trapping  level  at  0.287  ev  above  the  valence  band.  Moreover,  this  bom¬ 
bardment-induced  level  is  characterized  by  a  F  practically  equal  to  one; 
in  other  words,  the  trap  is  such  that 

y  .  exp  (-^)  “  1  ,  (19) 

where  y  refers  to  the  spin  situation  of  its  quantum  state  and  exp  (o/k)  to 
the  temperature  dependence  of  its  energy  level. 

As  recently  emphasized  by  Wertheim  in  connection  with  electron- 

bombardment  damage  studies,  carrier  concentration  measurements  can 

only  serve  to  determine  the  combined  effect  of  the  statistical  weight  ratio 

g 

and  the  energy-level  shift  .  This  is  strikingly  demonstrated  by  the  present 

analysis, but  the  knowledge  of  F  allows  an  accurate  evaluation  of  N^,  the 

total  density  of  those  defect  centers  leading  to  localized  states  at  an 

energy  above  the  valence  band.  This  solves  one  of  the  main  problems 

the  present  investigation  was  considering.  F  =  1  means  =  1.35N^  on 

the  basis  of  Eq.  (14),  if  77  =  0.  35.  In  other  words,  this  means  that  the 
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four-minute  exposure  of  sample  8b  (N  -  2.02  x  10  cm  )  in  the  BNL 

^  14 

reactor  has  resulted  in  the  introduction  of  2. 73  x  10  defect  centers 

3 

per  cm  ,  which  all  act  as  hole  traps  with  the  same  zero -temperature 
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depth  of  0.287  ev  . 

The  carrier  concentrations  measured  in  sample  8b  about  one 

month  after  irradiation  are  shown  in  Fig.  4.  The  previous  treatment 

again  yields  an  excellent  fit,  if  F  is  taken  as  equal  to  one.  It  should, 

nevertheless,  be  noted  that  the  low -temperature  data  now  points  towards 

an  activation  energy  of  0.291  ev  instead  of  0.287  ev,  with  7?  =0.30  ±0.01. 

The  implications  of  this  result  are  not  clear  at  the  present  time,  since 

the  observed  discrepancies  may  possibly  be  due  to  a  small  systematic 

error  during  one  of  the  runs.  Yet  the  overall  quality  of  the  data  leads 

us  to  believe  that  either  the  first  temperature  rurt  (though  limited  to  350®K) 

or  the  room -temperature  aging  process  actually  resulted  in  a  small  increase 

in  the  density  of  the  detected  traps  coupled  with  an  apparent  decrease 

in  the  trap  depth.  Similar  phenomena,  but  much  more  pronounced,  have 

2  3 

been  reported  by  Cleland  et  al.  for  neutron -bombarded  germanium  ’  . 

At  this  point  it  might  be  of  interest  to  consider  the  behavior  of  the 
Fermi-level  in  sample  8b,  as  exhibited  in  Fig.  5.  The  Fermi -energies 
determined  from  the  two  sets  of  carrier  concentration  data,  if  plotted 
versus  the  reciprocal  temperature,  seem  to  indicate  that,  after  crossing 
the  expected  trap  level  at  about  0.29  ev  from  the  valence  band,  the  Fermi- 
level  moves  up  again  with  decreasing  temperature.  It  will  be  shown  that 
this  remarkable  behavior  can  be  very  well  explained  on  the  basis  of  a  single  - 
level  model  with  the  previously  determined  characteristics . 

By  taking  Eqs.(2)  and  (10)  into  account,  Eq.  (15)  becomes 


2C(T)  exp  (-^)  = 


E^  -  E 

exp  ( 


(20) 


where  E^  is  the  Fermi -energy  and  E^  the  zero -temperature  level  of  the 

"^Evidently,  this  statement  does  not  imply  that  the  totality  of  the  radiation 
induced  defect  states  are  located  on  this  unique  level.  It  only  refers  to 
that  part  of  the  forbidden  gap  that  could  be  reasonably  explored  by  the 
carrier  concentration  analysis. 
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hole  trap,  both  measured  from  the  valence  band  edge.  Eq.  (20)  yields 
a  quadratic  equation  in  exp  (E^/kT),  the  solution  of  which  provides  the 
following  expression  for  the  Fermi-energy  as  a  function  of  temperature; 

E^  =  E^  +  kT  In  {  K  (T)  +  r?  +  ^[\C{T)+V]^  +  >  ,  (21) 

if  K  (T)  and  V  are  defined  according  to  Eqs.  (14)  and  (17).  At  low  tem¬ 
peratures  the  function  1C  (T)  becomes  exceedingly  small.  The  behavior 
of  the  Fermi-energy  is,  therefore,  essentially  fixed  by  the  value  of  the 
parameter  7? .  It  is  then  easy  to  realize  that  a  rise  of  the  Fermi-level 
by  lowering  the  temperature  can  only  be  brought  about,  if  0<  77  <  0.  5. 

(or  In (2  77  ) <  0)  .  Curves  obtained  with  77  =  0.2,  0.3,  and  0.4,  assuming 
r  =  1,  are  drawn  in  Fig.  5.  For  each  of  these,  the  trap  depth  was 
chosen  in  such  a  manner  as  to  fit  the  Fermi-energy  at  300 ®K  (E^  =  0, 308  ev), 
in  order  to  get  a  kind  of  "normalization"  making  a  selection  of  the  correct 
parameter  possible.  Clearly,  an  acceptable  overall  fit  is  provided  by 
V'VtO.'i,  a  value  which  corresponds  to  an  energy  level  at  0.289  ev  above 
the  valence  band.  This  appears  to  be  in  excellent  agreement  with  the 
results  we  gained  from  the  straightforward  analysis  of  the  carrier  concen¬ 
tration  data,  taking  into  consideration  the  fact  that  the  Fermi-level  analysis 
was  performed  simultaneously  on  both  sets  of  data.  It  should  also  be 
mentioned  that  fitting  attempts  with  F  I  were  unsuccessful,  and  in  this 
respect  the  Fermi-level  analysis  reveals  itself  as  an  especially  sensitive 
technique. 


B.  SAMPLE  6  A 

The  two -minute  exposure  of  sample  6a  resulted  in  an  estimated 

14/2 

damaging  dosage  of  1.2  x  10  n/cm  ,  or  half  the  dosage  given  to 
sample  8b.  We  may,  therefore,  expect  to  find  in  sample  8b  about  half 
as  many  defect  states  per  unit  volume,  especially  half  as  many  hole 
traps  of  the  type  detected  in  8b.  More  specifically,  if  we  assume  that 
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the  density  of  each  bombardment -induced  defect  goes  linearly  with 

the  integrated  flux,  we  expect  the  irradiation  of  sample  6a  to  result 

14  -3 

in  the  introduction  of  about  1.4  x  10  cm  of  those  defect  centers, 
which  the  foregoing  analysis  has  shown  to  give  rise  to  a  hole  trapping 
process  characterized  by  a  zero-temperature  energy  level  at  0.29  ev 
from  the  valence  band  and  a  F  -factor  of  one.  Before  exposure, 
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the  net  impurity  content  N  of  sample  6a  was  found  equal  to  1 . 88  x  10 
-3  ^ 

cm  .  Hence  we  expect  a  parameter  7?  of  about  -0.25,  and  on  the 
basis  of  Eq.  (21)  the  Fermi  level  should  cross  the  trap  level  and  move 
monotonically  towards  the  valence  band,  as  the  temperature  is  lowered. 
This  is  precisely  the  kind  of  behavior  exhibited  by  sample  6a  in  Fig.  2. 

The  successful  Fermi-level  analysis  of  sample  8b  (Fig.  5)  strongly 
suggests  that,  at  least  in  that  part  of  the  forbidden  gap  between  0.28  and 
0.36  ev  above  the  valence  band,  the  situation  is  correctly  described 
by  Eq.  (21)  with  F  =  1,  i.e. 

E^  =  E^  +  kT  In  {  K  (T)  +  77  +  /[K  (T)  +  77  +  4K(T)}  , 

(22) 

where  E^  is  the  zero -temperature  hole -trap  level,  a  level  which  is 
approximately  at  0.29  ev  but  whose  exact  position  might  depend  upon 
the  density  of  defect  centers,  density  which  is  characterized  by  the 
parameter  77  •  K  (T)  is,  of  course,  related  to  the  equilibrium  constant 
and  the  initial  carrier  concentration  according  to  Eq.  (17). 

A  close  examination  of  the  measured  Fermi -energies  in  sample 
6a  (Fig.  6)  reveals  that,  in  the  region  around  300 'K,  the  Fermi  level 
moves  linearly  with  temperature.  This  circumstance  allows  a  rather 
accurate  determination  of  the  Fermi-level  slope  at  300  ®K,  a  temperature 
at  which  E^  =  0.300  ev  and  should,  therefore,  be  determined  by  Eq.  (22). 
The  slope  yields  dE^/d(kT)  ]  =  9.90.  On  the  other  hand,  if 
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we  write 


X(T,  77)  =  K  (T)  +  77 


+  J[K{T)  +  77  ]^+  4K(T)  , 


(23) 


it  will  be  seen  that 


E 

-3^  =  in  U(T,n))  +  K  (T)  +  |] 


X  (T,  77  )  +  2 


X(T,77)  /lK(T)  +  77f  +4K:  (T) 

(24) 

Thus,  the  knowledge  of  and  dE^d  (kT)  at  300  “K  can  supply  the  para¬ 
meters  E^  and  77  of  sample  6a,  independently  of  any  previous  determin¬ 
ation,  as  a  result  of  the  simultaneous  consideration  of  Eqs.  (22)  and  (24). 
A  numerical  treatment  of  the  system  gave  E^  =  0.293  ev  and  77  =  -0.23, 
in  satisfactory  agreement  with  what  is  expected;  an  energy  level  at  about 
0.29  ev  from  the  valence  band  and  a  defect  density  apparently  proportional 
to  the  damaging  dosage  .  If  we  accept  the  quoted  values  and  put  them  into 
Eq.  (22),  the  resulting  curve  should 


a)  fit  the  measured  Fermi  energies  between  0.28  and  0. 36  ev,  if  the 
conclusions  of  the  8b  analysis  are  of  general  significance,  and 

b)  fit  all  the  measured  Fermi-energies,  if  our  6a  data  admit  of  an 
interpretation  in  terms  of  a  single -level  model. 


®It  should  be  realized,  however,  that  the  foregoing  procedure  cannot 
be  considered  as  a  reliable  one.  It  implies  exceedingly  accurate 
measurements  in  a  temperature  range  where  they  are  known  to  be 
difficult.  Moreover,  it  implies  that  with  a  Fermi  level  at  about 
0.30  ev  and  temperatures  near  300 'K  the  trapping  process  is  strictly 
due  to  a  single  discrete  level.  In  any  case,  an  accurate  determination 
of  E  can  only  be  made  on  the  basis  of  a  "true"  carrier  concentration 
slope.  This  requires  the  Fermi  level  to  be  pinned  down  on  the  level, 
in  other  words,  requires  a  parameter  77  which  is  positive(see  Eq.  13), 
and  this  condition  is  not  fulfilled  in  sample  6a. 
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In  fact,  Fig.  6  clearly  answers  yes  to  test  a  and  no  to  test  b.  A  dis¬ 
crete  energy-level  at  0.293  ev  provides  an  excellent  description  of 
the  situation  up  to  0.42  ev  frona  the  valence  band,  but  fails  as  soon 
as  the  Fermi  level  falls  below  0.28  ev.  Repeated  calculations  per¬ 
formed  with  somewhat  different,  but  still  acceptable  (E^,  7?  )  com¬ 
binations  only  helped  to  corroborate  this  evidence. 

The  carrier  concentrations  which  result  from  a  single -level 

approach  with  E,  =  0.293  ev,  71  =  -0.23,  and  F  =  1  are  plotted 
^  14  -3 

as  curve  1  in  Fig.  7.  1.45x10  cm  hole  traps  are  available  and 

they  will  all  be  filled  at  200  *K.  At  temperatures  below  200  ®K  we 
expect,  therefore,  the  bombardment -induced  deep  traps  to  manifest 
themselves  simply  by  increasing  the  compensation.  In  other  words, 
at  temperatures  below  200®K  sample  6a  should  behave  like  a  p-type 
specimen  with  a  net  impurity  content 

Nc  =  “c  -  =  -  ”  “c  ■ 

if  there  were  no  other  active  defect  states  in  the  lower  part  of  the 
energy  gap.  We  already  know  that  this  is  not  the  case.  At  low  tem¬ 
peratures  we  are  faced  with  a  new  trapping  process,  and  this  fact 
is  strikingly  borne  out  by  plotting  the  experimental  log  [p(T/300) 
values  versus  the  reciprocal  temperature  for  1000/T>7  (Fig.  7).  The 
straight  line  yields  an  activation  energy  of  0.159  ev.  In  other  words, 
our  6a  data  imply  the  existence  of  a  discrete  hole  trapping  level  at 
0.  159  ev  (zero-temperature  position)  from  the  valence  band.  In  the 
temperature  range,  where  all  the  deep  traps  are  filled,  and  N'  can 
be  considered  as  an  "initial  carrier  concentration",  we  are  then 
obviously  tempted  to  describe  the  situation  according  to  the  single - 
level  model,  that  is,  to  write 
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The  primed  symbols  now  refer  to  the  shallow  level  and  the  fraction  of  the 
initial  charge  carrier  concentration  which  was  not  frozen  out  by  the  deep 
traps.  For  p  <  <  ,  that  is  1000/T>7,  Eq.  (15'  )  leads  to 

_/m/,n«x-3/2  5.50x1018  ,  -3  z,,,, 

p(T/300)  =  -  exp(-i-7^)cm  ,  (13') 

77 

from  which  the  parameter  77  '  for  the  second  level  is  then  determined 
as  equal  to  10.  5  ±  0.  1,  if  E^'  =  0.159  ev. 


Repeating  the  procedure  adopted  for  sample  8b,  we  may  then  try 
to  determine  the  P  “  factor  of  the  shallow  level,  or  F  '  ,  by  calculating 


P 


N'F' 

2 


{- LK'  (T)  +  77']  +  /lk'(T)  +77]^  + 


(18' ) 


for  different  values  of  F  '  and  confronting  with  the  experimental  points 
in  the  temperature  range  where  Eq.  (15' )  might  be  meaningful.  Compu¬ 
tations  performed  with  F '  =  1  and  F'  =  4  resulted  both  in  virtually  the 

f 

same  curve,  the  curve  2  of  Fig.  7  .  Thus  it  became  apparent  that  our 
"Ansatz"  should  be  correct,  but  any  hope  to  determine  the  F  -factor 
for  the  0.16  level  on  the  basis  of  our  6a  data  alone  had  to  be  abandoned. 
This  unfortunate  state  of  affairs  is  due  to  the  large  value  of  77  '  for 
this  sample,  since  with  77'  ~  10  we  have  always 


LK '  (T)  +77' ]2  >>  (26) 

WhyF  '  =  1  or  4?  If  we  do  not  take  into  consideration  the  problematic 

temperature -shift  of  the  bombardment -induced  energy  levels,  the  F 
-factor  is  given  by  the  statistical  weight  factor  y  defined  as  the  ratio 
of  the  statistical  weights  of  the  reactants  to  that  of  the  product  of  the 
trapping  process.^'* It  is  probably  reasonable,  but  by  no  means  obvious, 
to  assume  that  y  can  take  a  value  of  either  unity  of  four  depending  on 
whether  the  state  occupied  by  the  hole  has  a  residual  electronic  spin 
or  not.  This  question  and  its  implications  in  respect  of  aF=  1,  as  found 
for  the  deep  trap,  will  be  discussed  in  a  forthcoming  paper. 
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for  any  reasonable  value  of  P  '  ,  and,  therefore, 


N’  K’  (T) 
c 

p  - 

K'  (T)  +  7?  ' 


(27) 


i.  e. ,  the  charge  carrier  concentration  is  practically  independent  of 
the  r  -factor  if  the  parameter  7?  becomes  much  larger  than  unity^. 


Under  those  circumstances,  it  is  easy  to  understand  why  even 
very  careful  fitting  attempts  performed  on  the  basis  of  a  simultaneous 
consideration  of  the  trapping  by  both  levels  (two -level  model'.),  and 
with  the  help  of  an  IBM-650  digital  computer,  failed  to  yield  the  true 
r  -factor  of  the  shallow  trap.  Nevertheless,  such  attempts  have 
led  to  a  very  accurate  determination  of  E^,  the  zero -temperature 
level  of  the  deep  trap  in  6a.  The  question  will  be  discussed  in  detail 
elsewhere,  but  it  can  be  seen  in  Fig.  8  that,  with  E^  =  0.289  ev  (a 
surprising  agreement  with  what  was  inferred  from  the  Fermi -level 
analysis  of  sample  8b'.),  the  theoretical  fit  provides  an  exact  repro¬ 
duction  of  the  experimental  data  at  all  temperatures .  A  general 
discussion, an  evaluation  of  the  results,  and  a  discussion  of  their 
implications  for  future  work  will  also  be  postponed  to  a  later  paper. 


^By  the  same  token  a  large  7?  '  entails,  even  at  relatively  high 
temperatures,  K  '  (T)  <  <  7?  '  and  p  ~  K'  (T)/77  '  ,  as  in  Eq.  (13'  ) 
which  provides  us  with  the  activation  energy.  There  is  no  need 
to  emphasize  that  the  detection  of  the  energy  level  is  appreciably 
facilitated  by  a  straight  slope  covering  an  extended  temperature 
range. 
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VI.  FIGURE  CAPTIONS 


Fig.  1: 

Fig.  2: 

Fig.  3: 

Fig.  4: 

Fig.  5; 

Fig.  6 

Fig.  7 


Charge  carrier  concentration  vs.  reciprocal  temperature 
before  and  after  pile  exposure  for  three  samples  cut  from 
a  60  Ohm-cm  p-type  silicon  crystal 

Fermi-energy  measured  from  the  top  of  the  valence  band 
vs.  temperature  before  and  after  pile  exposure  for  three 
samples  cut  from  a  60  Ohm-cm  p-type  silicon  crystal. 

Charge  carrier  concentration  vs.  reciprocal  temperature 
for  a  p-type  silicon  sample  with  an  initial  hole  density  of 
2.0  X  10^4  cm~^  after  a  2.4  x  10^4  n/ cm^  exposure  and 
six  months  room -temperature  aging.  The  dotted  straight 
line  fits  the  p(T/300)"V2  values  at  low  temperatures  and 
yields  the  activation  energy  The  theoretical  curves  in 
the  high  temperature  range  indicate  that  the  level  has 
a  r  -factor  practically  equal  to  one. 

Charge  carrier  concentration  vs .  reciprocal  temperature 
for  a  p-type  silicon  sample  with  an  initial  hole  density  of 
2.0  X  10^4  cm  after  a2.4xl0l4n/  cm^  exposure  and 
one  month  room -temperature  aging.  The  dotted  straight 
line  fits  the  p(T/ 300) "3/2  values  and  yields  the  activation 
energy.  The  curve  was  calculated  assuming  that  the  level 
has  a  r  -factor  equal  to  one. 

Fermi-energy  measured  from  the  top  of  the  valence  band 
vs.  reciprocal  temperature  for  a  p-type  silicon  sample  with 
an  initial  hole  density  of  2  0  x  10^4  cm"3  after  a  2.4x10^4 
n/cm^  exposure.  For  each  theoretical  curve  the  trap  depth 
and  trap  density  were  adjusted  so  as  to  fit  the  Fermi-energy 
at  300 ®K  (matching  point).  The  F  -factor  was  always  taken 
equal  to  one. 


Fermi-energy  measured  from  the  top  of  the  valence  band  vs. 
temperature  for  a  p-type  silicon  sample  with  an  initial  hole 
density  of  1 . 9  x  10i4  cm"^  after  a  1.2  x  10^4  n/cm^  exposure. 
The  theoretical  curve  was  obtained  by  fitting  the  Fermi-energy 
and  its  temperature  derivative  at  300 "K. 


Charge  carrier  concentration  vs.  reciprocal  temperature  for 
a  p-type  silicon  sample  with  an  initial  hole  density  of  1.9x10^'^ 
cm~^  after  1.2  x  10^4  ly  cm'^  exposure.  The  dotted  straight 
line  fits  the  p(T/ 300)"V2  values  at  low  temperatures  and 
yields  the  activation  energy  of  the  shallow  trap.  Curve  1 
corresponds  to  the  theoretical  Fermi  energies  of  Fig.  6..  and 
curve  2  was  calculated  according  to  Eq.  (27)  with  7?  '  =  10.  5. 
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Fig.  8:  Charge  carrier  concentration  vs.  reciprocal  temperature 
for  a  p'type  silicon  sample  with  an  initial  hole  density 
of  1.9  X  10^4  cm~3  after  a  1.2  x  lO^^  n/cm2  exposure. 

The  theoretical  curve  is  an  EDP  solution  of  a  cubic  equa¬ 
tion  in  p  resulting  from  two  uncorrelated  hole  trapping 
centers.  The  deep  trap  is  characterized  by  E+'  =  0. 289ev 
7?  =  -  0.21  and  r=  1;  the  shallow  one  by  E.'  =  0. 169ev 
and  77  '  =  10. 5.  ^ 
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EFFECT  OF  RADIATION  ON  THE  CRITICAL  SHEAR  STRESS 
OF  A  METAL  SINGLE  CRYSTAL 


by 

C.  E.  Morgan 
Conva Ir 

A  Division  of  General  Dynamics  Corporation 
Fort  Worth,  Texas 


A  mechanism  is  postulated  to  account  for  the 
radiation  hardening  of  a  metal  single  crystal.  The 
mechanism  is  based  on  a  lattice  defect  consisting 
of  interlocking  dislocation  rings.  According  to 
this  model,  the  critical  shear  stress  of  a  metal 
single  crystal  varies  as  the  cube  root  of  the  inte¬ 
grated  fast  neutron  flux.  This  agrees  well  with 
the  results  of  experiments  by  Blewltt  et  al -on 
Irradiated  high-purity  copper  crystals. 

INTRODUCTION 


The  critical  shear  stress  of  a  metal  single  crystal  is  great¬ 
ly  Increased  under  certain  conditions  by  Irradiation  with  neutrons . 
Such  an  effect  arises  from  structural  defects  Introduced  into  the 
lattice  of  the  metal  by  the  irradiation.  The  structural  defects 
that  supposedly  result  from  bombardment  of  a  metal  with  neutrons 
and  that  might  conceivably  harden  the  metal  are:  (l)  Frenkel  de¬ 
fects,  (2)  clumps  of  vacancies  produced  directly  by  the  irradia¬ 
tion,  (3)  dislocation  rings  produced  by  the  high  stresses  around 
thermal  and  displacement  spikes,  (4)  small  mlsoriented  regions  left 
behind  by  displacement  spikes,  and  (5)  clustered  vacancies  or  inter¬ 
stitials  formed  by  diffusion  of  the  Frenkel  defects.  It  is  of  in¬ 
terest  to  determine,  if  possible,  which  of  these  defects  cause  the 
increase  in  critical  shear  stress  of  an  Irradiated  metal. 
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DISCUSSION  OF  PREVIOUS  WORK 


Frenkel  defects  would  have  to  be  extremely  effective  in 
interacting  with  the  glide  dislocations  to  produce  the  large 
amount  of  hardening  that  is  observed  for  irradiated  metals. 

For  example,  it  is  apparent  that  the  Frenkel  defects  begin  to 
anneal  in  copper  at  temperatures  very  low  compared  to  the  tem¬ 
perature  at  which  the  radiation-induced  increase  in  critical 
shear  stress  anneals.  In  copper  the  interstitials  begin  to 
anneal  rapidly  at  30°K  and  the  vacancies  anneal  at  about  290°K. 
The  increase  in  critical  shear  stress  does  not  anneal  until  a 
temperature  of  about  550°K  is  reached.  Quenching  experiments^ 
show  that  randomly  dispersed  vacancies  do  not  appreciably  harden 
a  metal]  hardening  occurs  only  when  the  specimen  has  been  warm¬ 
ed  up  to  a  temperature  at  which  the  vacancies  can  migrate  to 
dislocations  or  else  cluster  to  form  voids.  Therefore,  it 
appears  that  randomly  dispersed  Frenkel  defects  cannot  be  the 
chief  cause  of  the  radiation-hardening. 

If  the  Frenkel  defects  form  Jogs  on  the  glide  dislocations, 
hardening  of  the  metal  would  be  expected,  since  the  jogs  effect¬ 
ively  increase  the  Peierls  force  for  the  dislocation.  There  are 
a  number  of  objections  to  this  view,  however.  Since  the  acti¬ 
vation  energy  due  to  the  increased  Peierls  force  is  not  likely 
to  be  very  large,  even  when  the  applied  stress  is  zero,  the  jog¬ 
ged  slip  dislocations  should  be  practically  unrestrained  at  any 
but  the  lowest  temperatures.  Thompson  and  Holmes^  have  observ¬ 
ed  that  the  internal  friction  of  copper  single  crystals  decreases 
after  reactor  irradiation  and  that  this  decrease  saturates  at  the 
very  low  dose  of  4(lol2)nvt.  Barnes  and  Hancock^  find  that  the 
internal  friction  of  copper  single  crystals  remains  unaffected  by 
a  short  neutron  bombardment  at  yS^K,  if  no  warming  is  allov;ed  be¬ 
fore  measurement.  Successive  pulse  annealings  produce  no  marked 
change  until  300°K  is  reached,  when  the  internal  friction  is 
greatly  reduced.  Since  randomly  dispersed  Frenkel  defects  should 
have  little  effect  on  the  damping,  the  decrease  of  internal 
friction  is  undoubtedly  caused  by  pinning  of  the  dislocations  by 
the  Frenkel  defects.  In  view  of  tne  fact  that  the  decrease  of 
internal  friction  and  the  change  in  elastic  constants  are  com¬ 
pletely  saturated  at  10^5  nvt,  and  that  Blewitt^  observes  no  ra¬ 
diation  annealing  or  saturation  of  the  increase  in  critical  shear 
stress  for  integrated  fluxes  up  to  10^®  nvt,  the  radiation-harden¬ 
ing  cannot  be  due  to  pinning  of  dislocations  by  point  defects. 
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That  the  Frenkel  defects  cluster  by  diffusion  to  form  ag¬ 
gregates  of  interstitials  or  to  form  voids,  tiowever,  remains  a 
possibility.  Such  clusters  could  be  formed  on  dislocations  or 
at  random  throughout  the  lattice.  In  either  case,  they  would 
be  very  effective  barriers  to  the  movement  of  slip  dislocations, 
provided  that  they  are  larger  than  a  few  atoms.  There  are  a 
number  of  important  difficulties,  however.  For  example,  if 
copper  is  irradiated  at  a  temperature  at  which  the  point  defects 
are  mobile,  it  hardens  at  the  maximum  rate  consistent  with  the 
bombarding  flux,  provided  no  self -diffusion  takes  place.  This 
indicates  that,  since  the  point  defects  can  easily  reach  the 
dislocations  and  disappear  at  jogs  on  the  dislocation  line,  the 
metal  is  hardened  even  when  a  high  steady-state  concentration 
of  point  defects  necessary  to  nucleate  aggregates  cannot  be 
built  up.  Quenching  experiments5  on  nickel  and  aluminum  show 
that  during  annealing  of  the  quenched-in  vacancies  the  volume 
of  the  metal  decreases  simultaneously  with  its  resistivity. 
Indicating  that  the  vacancies  are  actually  disappearing  and  not 
merely  aggregating  into  voids.  Finally,  Blewittl  has  found  that 
a  single  crystal  copper  specimen  bombarded  at  2Q°K  with  neutrons 
and  immediately  tested  at  20OK  has  the  same  critical  shear  stress 
as  specimens  bombarded  at  20°K,  then  pulse  annealed  to  tempera¬ 
tures  as  high  as  250°K,  and  at  last  tested  at  20°K.  It  will  be 
recalled  that  no  diffusion  of  structural  defects  takes  place  in 
copper  until  a  temperature  somewhat  above  20°K  Is  reached;  in 
fact,  the  interstitials  anneal  rapidly  between  30°K  and  40°K. 
Thus,  the  point  defects  cannot  aggregate  at  20°K,  and  yet  the 
copper  specimen  is  fully  hardened  by  the  irradiation. 

Hence,  one  is  led  to  the  conclusion  that  radiation-harden¬ 
ing  can  only  be  caused  by  a  defect  formed  directly  by  the  ir¬ 
radiation.  A  hardening  mechanism  based  on  the  Frenkel  defects 
is  unsatisfactory  for  the  reasons  which  have  been  mentioned. 

Brinkraan^has  suggested  that  small  misoriented  regions 
similar  to  small  grains  may  be  left  behind  when  a  displacement 
spike  resolidifies.  Since  the  freezing  of  the  molten  core  of 
the  displacement  spike  proceeds  from  the  outside  in  with  the 
surrounding  lattice  as  a  site  for  nucleation,  the  resolidified 
region  should  have  the  same  orientation  as  the  remainder  of  the 
lattice.  Thus,  one  would  not  expect  very  many  such  misoriented 
regions  to  be  formed  by  irradiation. 
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It  is  difficult  to  see  how  neutron  bombardment  could  direct¬ 
ly  produce  clusters  of  Interstitial  atoms  merely  by  displacing 
lattice  atoms.  On  the  other  hand,  it  is  quite  likely  that  small, 
stable,  multiple  vacancies  might  be  produced  directly  by  neutron 
bombardment.  According  to  Brlnkman°,  a  primary  knock-on  will  pro¬ 
duce  s+^able  multiple  vacancies  in  copper  when  its  kinetic  energy 
is  40C  ev  or  less.  Multiple  vacancies  produced  by  knock-ons 
with  higher  energy  than  this  will  be  unstable  and  collapse  to  form 
a  displacement  spike.  Snyder  and  Neufeld’T^  have  shown  that  the 
average  total  number  of  displaced  atoms  including  the  primary 
knock-on  is  approximately  where  E  is  the  energy  of  the  knock- 

on  and  Eq  Is  the  average  threshold  displacement  energy.  The 
threshold  displacement  energy  is  about  25  ev.  Therefore,  the 
largest  stable  multiple  vacancy  that  can  be  produced  directly  by 
bombardment  of  copper  is  about  eight  atomic  volumes.  There  is 
some  question  whether  or  not  such  small  cavities  can  cause  much 
hardening,  because  they  are  too  small  to  exert  a  sizeable  image 
force  on  a  glide  dislocation.  If  they  did,  the  resulting  hardening 
should  be  practically  Insensitive  to  temperature,  in  contrast  to 
experimental  results.  The  strain  field  around  the  multiple 
vacancy  is  probably  not  large.  That  means,  there  should  be  practi¬ 
cally  no  interaction  between  glide  dislocations  and  multiple  va¬ 
cancies  except  when  they  are  separated  by  only  a  few  interatomic 
distances.  If  an  edge  dislocation  passes  through  a  multiple  va¬ 
cancy,  a  Jog  will  be  formed  on  it.  These  Jogs  may  not  seriously 
Impede  edge  dislocations.  Jogs  are  probably  not  formed  on  screw 
dislocations  by  the  multiple  vacancies.  Moreover,  in  this  case 
it  may  not  be  necessary  for  the  partial  dislocations  of  an  extend¬ 
ed  edge  dislocation  to  come  together  to  cut  through  the  multiple 
vacancy,  since  complex  stacking  faults  may  not  necessarily  be  form¬ 
ed  by  individual  climbing  of  the  partlals.  If  they  are  formed, 
they  will  be  only  about  4  atoms  in  size.  If  multiple  vacancies 
are  effective  barriers  for  the  glide  dislocations,  they  might  well 
be  the  cause  of  radiation-hardening.  Holmes^  at  Oak  Ridge  has 
advanced  a  theory  that  radiation  hardening  is  caused  by  clumps  of 
point  defects.  It  appears  to  be  difficult  to  explain,  on  the 
basis  of  clustered  point  defects,  the  fact  that  work-hardened 
metals  are  not  hardened  nearly  as  much  as  annealed  metals  by 
irradiation. 


THE  THEORY 

The  possibility  that  small  shear  or  prismatic  dislocation 
rings  produced  by  the  high  stresses  around  thermal  and  displacement 
spikes  might  be  the  cause  of  hardening  seems  to  have  been  overlook¬ 
ed  so  far.  According  to  Seitz  and  Koehler°,  the  atoms  in  the  heat¬ 
ed  zone  of  a  spike  exert  a  pressure  on  the  rest  of  the  lattice  of 
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the  order  of  the  shear  modulus.  These  Investigators  show  that 
the  probability  of  a  dislocation  lying  in  the  region  about  the 
spike,  where  the  stress  is  high  enough  to  move  dislocations,  is 
negligible  for  well-annealed  metals.  Therefore,  this  stress 
will  probably  not  be  relieved  by  gliding  of  the  dislocations 
already  present  in  the  metal.  Consequently,  small  dislocation 
rings  should  be  formed  spontaneously.  If  shear  dislocation 
rings  are  formed  by  irreversible  slip  on  intersecting  glide 
planes  so  that  the  rings  interlock,  the  ring  may  not  collapse 
completely  after  dissipation  of  the  spike.  In  fact,  it  seems 
probable  that  multiple  slip  occurs  on  nonparallel  glide  planes 
even  for  h.c.p.  metals  because  of  the  exceedingly  large  stresses 
around  the  spike.  As  a  result  of  intersection  of  the  dislocation 
rings.  Jogs  will  be  on  opposite  sides  of  the  dislocation  rings. 

When  the  temperature  of  the  metal  la  below  that  required  for 
rapid  self -diffusion,  the  Jogs  prevent  the  collapse  of  the  dis¬ 
location  rings.  Since  Jogs  on  the  dislocation  rings  are  edge 
dislocations  and  since  there  la  a  force  of  attraction  between 
opposite  sides  of  the  dislocation  rings,  diffusion  of  vacancies  to 
or  from  the  Jogs  on  the  dislocation  rings  will  cause  these  Jogs 
to  climb.  The  dislocation  rings  will  then  disappear  from  the 
lattice  at  temperatures  in  the  self -diffusion  range.  An  alternative 
is  that  prismatic  dislocations  are  produced  by  the  spike  instead  of 
shear  dislocation  rings.  Each  of  these  will  probably  consist  of 
an  extra  half -sheet  of  material  which  is  roughly  circular  in 
shape  and  have  an  edge  dislocation  all  the  way  around  the  edge  of 
the  extra  half -sheet.  This  prismatic  dislocation  ring  will  also 
be  removed  from  the  lattice  by  climb.  Thus,  the  recovery  of  the 
critical  shear  stress  of  an  Irradiated  metal  will  take  place  by 
self -diffusion. 

Prom  the  phenomenon  of  work-hardening,  we  know  that  the  pre¬ 
sence  of  other  dislocations  in  the  lattice  can  seriously  Impede 
the  glide  dislocations  and  thereby  cause  a  large  amount  of  harden¬ 
ing.  ^’This  is  not  to  suggest  that  radiation-hardening  is  like 
work- hardening. )  The  Jogs  on  the  shear  dislocation  rings  should 
preveiit  collapse  or  expansion  of  the  dislocation  rings.  At  any 
rate,  in  order  to  expand  a  small  dislocation  ring  with  no  Jogs  on 
it,  the  applied  stress  will  have  to  oppose  a  very  large  stress  act¬ 
ing  to  collapse  the  ring.  For  example,  in  the  case  of  a  shear  dis¬ 
location  ring  in  copper  with  a  radius  of  20  interatomic  distances, 
a  stress  of  about  200  kg/mm^  is  needed  to  expand  it.  Prismatic 
dislocation  rings  of  this  size  should  be  practically  immobile. 
Therefore,  it  seems  likely  that  the  dislocation  rlnfes  will  be 
immobile  during  deformation  of  the  metal. 
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Now,  consider  an  irradiated  metal  in  which  the  density  of 
dislocation  ring  groups  is  N.  Two  ways  by  means  of  which  the 
principal  glide  dislocations  interact  with  these  groups  may  be 
envisioned.  First,  the  strain  field  around  the  dislocation 
ring  groups  may  hold  up  the  glide  dislocations  so  that  they  do 
not  cut  through  the  dislocation  rings  at  the  stress  level  dur¬ 
ing  deformation.  Since  the  dislocation  rings  are  small  and 
their  strain  field  does  not  extend  beyond  a  distance  of  the 
order  of  magnitude  of  the  radius  of  the  dislocation  ring,  the 
glide  dislocations  can  bow  out  between  them  In  the  manner  sug¬ 
gested  by  0rowan9  for  overaged  alloys.  In  this  case  the  criti¬ 
cal  shear  stress  will  be  inversely  proportional  to  the  average 
distance  between  dislocation  ring  groups,  or  directly  propor¬ 
tional  to  nV3,  The  critical  shear  stress  for  this  case  Is 
given  by: 


or 

z  GbN^/3^  (l) 

where  G  is  the  shear  modulus  and  b  is  the  Burgers  vector.  This 
type  of  hardening,  however.  Is  not  considered  to  be  very  sen¬ 
sitive  to  temperature. 


The  second  kind  of  interaction  would  occur  if  the  glide  dis¬ 
locations  cut  through  the  dislocation  rings.  When  the  glide  dis¬ 
locations  cut  through  the  dislocation  rings.  Jogs  will  be  formed 
on  them.  The  energy  required  to  form  a  Jog  will  beo<^Gb3,  where 
o<^l.  The  average  distance  between  dislocation  ring  groups  will 
be  N-V3,  and,  if  we  add  the  Jog  energy  and  the  contribution  from 
the  applied  stress,  we  obtain  the  activation  energy 

,  Obi 

aQb3  .  ;;i73 


For  plastic  flow  to  be  rapid,  the  flow  stress  must  be 

Cf  =  GbNl/3 


(3) 


which  is  the  same  as  Equation  1.  In  this  case,  the  flow  stress 
should  be  temperature -dependent .  In  order  to  cut  through  the  dis¬ 
location  ring,  however,  the  partial  dislocations  in  h.c.p.  and 
f.c.c.  metals  will  have  to  coalesce;  otherwise,  a  complex  stacking 
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fault  at  the  place  of  Intersection  will  be  formed.  The  stress 
necessary  to  force  the  two  partial  dislocations  of  an  extended 
dislocation  together  is  quite  high  even  for  metals  with  high 
stacking  fault  energies,  such  as  aluminum.  The  dependence  of 
the  width  of  an  extended  dislocation  on  stress  during  inter¬ 
section  of  dislocations  has  not  been  considered,  but  this  will 
affect  the  activation  energy  U. 

Equation  2  indicates  that  the  critical  shear  stress  at  con¬ 
stant  strain  rate  will  be  directly  proportional  to  the  absolute 
temperature  of  the  specimen.  In  the  case  of  metals  with  low 
stacking  fault  energy  such  as  copper,  the  temperature -dependence 
of  the  critical  shear  stress  will  be  more  complicated,  since 
different  rate -determining  processes  occur  in  different  tempera¬ 
ture  intervals.  Indeed,  Blewitt^O  observes  that  the  critical 
shear  stress  of  copper  single  crystals  Irradiated  with  neutrons 
is  approximately  proportional  to  the  square  root  of  the  absolute 
temperature  of  the  specimen. 

DISCUSSION  OP  THE  THEORY 


The  number  of  dislocation  ring  groups  should  be  proportional 
to  the  number  of  temperature  spikes  and,  therefore,  proportional 
to  the  integrated  fast  neutron  flux.  The  critical  shear  stress 
from  Equation  1  or  Equation  3,  then,  is  proportional  to  the  cube 
root  of  the  integrated  fast  neutron  flux,  in  agreement  with  the 
results  that  Blewlttl^  has  obtained  for  99.999  percent  pure 
copper  single  crystals  irradiated  up  to  an  Integrated  fast  neu¬ 
tron  flux  of  2(l0l9)nvt.  Calculation  of  the  number  of  displace¬ 
ment  spikes  produced  during  bombardment  shows  that  there  should 
be  enough  of  the  dislocation  ring  groups  to  account  for  the 
amount  of  hardening  observed.  Thus,  it  is  possible  to  explain  the 
general  features  of  radiation-hardening  on  the  basis  of  the  dis¬ 
location  rings  caused  by  temperature  spikes.  The  value  of  a  new 
theory  is  determined  by  how  well  it  can  account  for  observed  facts 
and  how  well  it  can  predict  experimentally  verifiable  phenomena. 
Some  of  the  effects  to  be  expected  on  the  basis  of  the  theory  are 
as  follows; 

1,  The  generation  of  dislocation  rings  around  temperature 
spikes  depends  on  the  existence  of  the  high  stresses  in 
the  vicinity  of  the  spike.  Furthermore,  a  glide  disloca¬ 
tion  in  this  high  stress  region  will  relieve  the  stress 
by  moving.  Therefore,  no  dislocation  rings  will  be  form¬ 
ed  if  a  dislocation  happens  to  be  in  the  high  stress 
region.  The  more  a  metal  has  been  work -hardened,  the  more 
chance  there  is  for  a  dislocation  to  be  in  the  high  stress 
region.  For  a  fully  work -hardened  metal,  exceedingly  few 


dislocation  rings  will  be  produced  by  the  tempera¬ 
ture  spikes.  Thus,  it  is  possible  to  explain  on  the 
basis  of  the  dislocation  ring  mechanism  the  fact  that 
irradiation  of  a  work -hardened  metal  raises  the  flow 
stress  by  a  small  amount  compared  with  the  change  due 
to  irradiation  of  an  annealed  metal. 

2.  Since  the  dislocation  rings  are  produced  directly  by 
the  Irradiation,  radiation-hardening  of  a  metal  should 
take  place  at  the  lowest  temperatures. 

3.  Because  the  dislocation  rings  will  disappear  from  the 
lattice  by  climbing  of  the  Jogs  on  shear  dislocation 
rings  or  by  climbing  of  the  small  prismatic  dislocations 
recovery  of  the  critical  shear  stress  to  its  pre-irra¬ 
diation  value  should  be  by  self -diffusion. 

4.  Only  heavy,  energetic  particles  will  produce  thermal 
spikes.  Therefore,  bombardment  of  a  metal  with  par¬ 
ticles  such  as  fast  neutrons  should  result  in  con¬ 
siderable  hardening  of  the  metal.  Bombardment  with 
electrons  that  produce  single  atomic  displacements 
only  or  with  gamma  rays  should  produce  at  most  a  slight 
hardening  of  the  metal. 

5.  For  annealed  metals,  there  should  be  no  saturation  of 
the  increase  in  critical  shear  stress  with  Integrated 
flux,  up  to  very  high  levels  of  irradiation.  Likewise, 
there  should  be  practically  no  radiation  annealing. 

6.  There  should  be  little  radiation-hardening  until  either 
the  separation  distance  of  the  dislocation  ring  groups 
becomes  less  than  the  length  of  the  Frank-Read  sources 
in  the  metal,  or  the  effect  of  the  dislocation  rings 
becomes  large  enough  to  exceed  the  normal  slip  resistance 
of  the  lattice. 

7.  The  radiation-induced  increase  in  the  critical  shear 
stress  will  be  sensitive  to  temperature  if  the  glide 
dislocations  cut  through  the  dislocation  rings. 

8.  According  to  the  dislocation  ring  model,  growth  of  slip 
lines,  under  certain  conditions,  should  be  more  diffi¬ 
cult  in  an  irradiated  metal;  i.e.,  there  should  be  a 
lattice -hardening.  Some  metals,  for  example  nickel, 
seem  to  show  this  effect. 
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The  coarse  slip  lines  and  the  yield  point  effect  sometimes 
observed  in  irradiated  metals  indicate  that  nucleation  of  slip 
lines  is  more  difficult  than  their  growth.  Ordinarily j  this 
would  not  be  expected  from  the  dislocation  ring  model.  However, 
if  a  glide  dislocation  cuts  through  a  dislocation  ring  group  in 
a  manner  such  that  point  defects  are  left  behind  by  the  Jogs  so 
formed  on  the  glide  dislocation,  the  next  glide  dislocation  that 
cuts  through  the  dislocation  ring  group  will  lay  down  a  row  of 
point  defects  next  to  the  row  left  by  the  first  dislocation. 

Since  the  formation  of  rows  of  point  defects  by  succeeding  glide 
dislocations  under  these  conditions  takes  less  and  less  energy, 
there  would  be  a  drop  In  stress  corresponding  to  the  yield  point 
effect.  This  should  happen  only  under  conditions  In  which  the 
glide  dislocation  can  produce  point  defects  as  a  result  of  being 
Jogged  and  In  which  the  strain  rate  Is  high  enough.  Other  pro¬ 
cesses  such  as  cross-slip  or  pinning  by  point  defects  may  be 
Involved  and,  since  detailed  observations  of  the  yield  point 
effect  have  not  been  made,  the  above  explanation  must  be  regard¬ 
ed  as  tentative. 

The  strength  of  metal  single  crystals  after  extensive  plastic 
flow  Is  nearly  the  same  for  both  Irradiated  and  nonlrradiated 
specimens.  Moreover,  the  work -hardening  rate  at  this  point  be¬ 
comes  nearly  the  same.  This  is  due  to  the  mutual  interference 
of  glide  dislocations.  As  deformation  proceeds,  more  glide  dis¬ 
locations  are  produced  and  the  mobility  of  the  dislocations  be¬ 
comes  more  sensitive  to  the  presence  of  other  glide  dislocations , 
than  to  dislocations  Introduced  by  Irradiation. 

The  purpose  of  this  paper  is  to  point  out  that  most  of  the 
observed  effects  of  Irradiation  on  the  critical  shear  stress  of 
a  metal  single  crystal  can  be  explained  by  a  model  of  radiation¬ 
hardening  caused  by  small  dislocation  rings  Introduced  into  the 
lattice  of  a  metal  by  the  irradiation.  It  is  difficult  to  see 
how  any  other  defect  that  might  be  produced  by  Irradiation  could 
be  responsible  for  both  effect  1  and  effect  2  mentioned  above. 

The  small  amount  of  radiation-hardening  displayed  by  a  work-hard¬ 
ened  metal  might  be  explained  by  a  theory  based  on  some  other 
radiation-induced  structural  defect,  if  it  is  assumed  that  the 
number  of  dislocations  due  to  the  work -hardening  is  considerably 
larger  than  the  number  of  such  defects  caused  by  irradiation. 

This  implies  that  the  strength  of  a  highly  work -hardened  metal 
could  be  considerably  Increased  by  Irradiating  it  to  a  sufficient¬ 
ly  high  level,  while  the  dislocation-ring  model  indicates  that 
the  strength  of  a  fully  work -hardened  metal  will  be  Increased 
only  a  little  for  any  amount  of  irradiation. 
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MEASUREMENT  OF  THE  RANGE  OF  RECOIL  ATOMS 

Ly 

R.  A.  Schmitt  and  R.  A.  Sharp 

John  Jay  Hopkins  Laboratory  for  Pure  and  Applied  Science, 

General  Atomic  Division  of  General  Dynamics  Corporation, 

San  Diego,  California 

An  important  problem  in  the  interpretation  of  radiation-damage  and 
sputtering  phenomena  is  the  evaluation  of  the  range  of  an  atom  which  moves 
through  a  lattice  after  having  received  an  initial  energy  of  10  to  100  kev. 
An  initial  energy  of  about  this  amormt  is  acquired  by  an  atom  which  has  been 
struck  by  a  fast  neutron  (^^1  Mev).  The  kinetic  energy  of  such  atoms  will 
be  transferred  to  the  lattice  largely  by  thermal  vibrations  and  the  creation 
of  lattice  vacancies  and  displacements.  If  the  collision  cross  section  is 
very  large,  the  displaced  atoms  and/or  vacancies  will  be  separated  by  only 
an  atomic  lattice  spacing.  On  the  other  hand,  separations  of  many  atomic 
diameters  are  expected  for  low  collision  cross  sections.  Annealing  of 
such  a  perturbed  lattice  is  achieved  by  the  diffusion  of  lattice  vacancies 
and  interstitial  atoms . 


A  theoretical  evaluation  of  the  collision  cross  section  is  difficvilt 
becaiise  (l)  such  a  calculation  essentially  involves  a  many-body  problem  of 
the  interpenetration  of  orbital  electron  clouds  during  collision,  and  (2) 
the  atomic  charges  of  moving  atoms  through  solid  lattices  are  unlaiown. 
Everhart,  et  al.,^~3  have  measured  the  differential  cross  sections  for  the 
scattering  of  singly  ionized  noble  gases  (25  to  I38  kev)  in  noble  gases 
and  also  have  determined  the  charge  states  of  the  scattered  atoms  after 
single  collisions.  Since  a  crystal  lattice  is  a  condensed  medium  and 
atomic  stopping  power  is  larger  in  crystalline  lattices  than  in  gaseous 
media, ^  it  is  doubtfxal  whether  the  scattering  data  of  gaseous  atomic 
interactions  may  be  used  to  theoretically  calculate  scattering  cross 
sections  in  solids. 


Experimentally,  a  direct  observation  of  recoil  atoms  in  solids  is 
difficult  because  of  the  uncertain  charge  states.  At  present,  tb.;  motion 
of  an  uncharged  atom  in  a  lattice  does  not  produce  effects  that  are  amenable 
to  detection.  Only  by  measxrrements  of  density,  electrical  conductivity, 
elastic  moduli,  etc.,  are  we  able  to  study  the  effects  of  radiation  damage. 5 
A  study  of  the  range  of  recoil  atoms  is  one  vmknown  in  the  general  problem 
of  radiation  damage. 

Investigations  of  the  ranges  in  metals  and  gases  of  the  light  atom, 
helium,  in  its  two  charge  states  have  been  reported. 8  Also,  considerable 
work  has  been  devoted  to  measurements  in  solids  of  highly  energetic  (in 
the  tens -of -Mev  region)  fission  products. T-9 
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By  using  an  electromagnetic  isotope  separator,  Thulin  '  has  measured 
the  penetration  depth  of  50-lcev  Xe^-35  atoms  in  Formvar  absorbers.  Following; 
the  discussion  of  Bohr^^  in  his  classic  paper  on  atomic-particle  penetration, 
Nielsen^^  shown  that  the  calculated  ranges  agreed  fairly  well  with  the 
observed  ranges  found  by  Thulin.  Furthermore,  Nielsen  successfully  correlates 
the  depth  of  ion  penetration  of  a  nijmber  of  ions  into  various  foils  with  the 
width  of  the  energy  distribution  of  light  charged  particles  scattered  from 
the  ions  deposited  in  the  foils. 13 

A  novel  experimental  technique  has  been  applied  to  the  problem  of  de¬ 
termining  the  ranges  of  atoms  in  the  kev  region.  Atoms  with  initial  energies 
in  this  region  are  produced  by  irradiating  suitable  targets  with  high-energy 
bremsstrahlung  gamma  rays,  and  the  products  of  photonuclear  reactions,  such 
as  photoneutron  (7,  n)  transmutations,  are  observed.  The  bremsstrahlung 
gamma  rays  used  in  this  work  were  generated  by  a  2i4-Mev  Allis -Chalmers 
betatron.®  The  photon  spectrm  of  gamma  rays  is  continuous  from  0  to  2h-  Mev; 
the  energy  spectra  of  evaporated  neutrons  from  copper!^  and  leadl5  are  cen¬ 
tered  at  about  I.5  Mev,  with  a  full  width  at  half  maximum  of  about; 

2  Mev.  Since  (l)  the  energy  thresholds  for  the  photoneutron  reactions  of 
the  ff.otjority  of  the  elements  reported  in  Table  I,  below,  are  well  below 
2k  Mev  and  near  the  copper  and  lead  photoneutron  thresholds  and  (2)  all  the 
photonuclear  absorption  cross  sections  of  elements  in  Table  I  nearly  resemble 
the  familiar  bell-shaped  or  resonance  cross-section  curve  of  copper  and  lead, 
it  is  anticipated  that  the  spectra  of  evaporated  neutrons  for  all  the  targets 
studied  here  are  centered  at  about  1.5  Mev,  with  corresponding  half -widths 
of  2  Mev. 


'The  nucleus  receives  a  negligible  amount  of  momentum  from  the  incoming 
gamma  ray  (^2  kev  and  -^O.T  kev  for  a  l6-Mev  photon  absorbed  by  copper  and 
gold,  respectively).  However,  the  nucleus  (atom)  recoils  from  neutron 
emission  with  an  energy  Ep,  given  by  the  formula 


E  -E 

M. 

A 


where  Ejj  is  the  neutron  energy  and  Mjj  and  M;^  are  the  masses  of  the  neutron 
and  atom,  respectively.  For  representative  atoms,  such  as  fluorine  (A  =  19), 
copper  (a  =  63  and  65),  and  goxd  (A  =  197),  the  most  probable  recoil  energies 
are  about  75^  25,  and  10  kev,  respectively.  As  a  result  of  photoneutron 
emission,  the  nuclei  of  the  recoiling  atoms  are  usually  unstable  to  beta 
(positron)  decay  and,  therefore,  are  amenable  to  detection  by  standard  scin¬ 
tillation  and  ceta-counting  methods. 

The  ranges  were  determined  for  carbon  in  polystyrene,  CH),  fluorine 

(plo  Teflon,  CF2),  chlorine  (Cl3^*-  in  saran,  CHCl),  and  the  metals  titanium, 
iron,  zinc,  copper,  molybdenum,  silver,  and  gold  in  their  respective  metallic 

We  wish  to  acknowledge  the  generous  cooperation  of  Dr.  Waldo  K.  I^^on 
and  Mr.  R.  B.  Doherty  of  the  U.S.  Navy  Electronics  Laboratory,  San  Diego, 
California,  in  the  performance  of  these  irradiations. 
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lattices .  The  method  used  was  somewhat  similar  to  the  stacked-foil  technique 
commonly  used  in  charged-particle  excitation-function  studies.  Thin  ele¬ 
mental  or  plastic  target  foils  about  0.001  in.  thick  were  interleaved  between 
0.0008-in. -thick  aluminum  catchers  (Alcoa  1199-0^  99-986^  pure).  Before  use, 
the  metallic  foils  were  thoroughly  cleaned  with  organic  solvents  and  etched 
with  appropriate  acids;  plastic  foils  were  detergent -cleaned.  In  order  to 
evaluate  the  background  activities  present  in  the  aluminum  catchers,  an  eq\ial 
nimiber  of  alminum  foils  without  targets  were  simultaneously  Irradiated. 

About  2000  roentgens  per  minute  of  23-Mev  bremsstrahlung  were  incident  on 
the  sandwiched  targets . 

Immediately  after  Irradiation,  the  sandwich  was  disassembled  and  the 
radioactivities  of  the  target  foils,  aluminum  catcher  foils,  and  aluminum 
background-monitoring  foils  were  measured  by  scintillation-coimting  the 
annihilation  radiation  (0.5I  Mev)  of  the  positrons.  The  recoiling  atoms 
which  are  close  to  the  target-foil  surface  are  cavight  in  the  aluminum  foils. 
If  the  range  of  the  recoil  atoms  is  small  compared  with  the  foil  thickness, 
the  range,  R,  is  related  to  the  fraction  N^/Nf  of  atoms  which  recoil  out  of 
the  target  foil  and  to  the  foil  thickness,  x,  by  R  =  2(Nc/Nf)x,  where 
and  Nf  are  the  numbers  of  atoms  that  recoil  into  the  catcher  and  that  are 
retained  in  the  target  foil,  respectively.  By  counting  nearly  equal  thick¬ 
nesses  of  target  foils  and  aluminxm  catchers,  the  errors  in  the  ratio 
Ng/Nf  are  less  than  %.  Subsidiary  experiments  showed  that  scattering  of 
the  recoil  atoms  back  into  the  target  foils  by  air  between  the  target  foils 
and  aluminimi  catchers  is  negligible. 


In  Table  I  are  listed  the  experimental  data  and  other  information  rele¬ 
vant  to  this  paper.  The  first  column  lists  the  recoiling  atoms  and  their 
mass  nvimbers.  The  most  probable  recoil  energy  (in  kev)  is  given  in  the 
second  column,  with  the  average  recoil  energies  for  all  atoms  assumed  to  be 
about  1.5  times  the  most  probable  recoil  energies,  as  was  observed  for  the 
copper  and  lead  photoneutron  spectra.  The  third  coliimn  gives  the  fraction 
of  activated  atoms  that  recoil  out  of  a  0.001 -in.  foil  with  23-Mev  bremsstrah¬ 
lung  incident  on  the  foils,  (A  graphic  representation  of  Ng/Nf  is  displayed 
in  Fig.  1,  below.)  In  the  fo\irth  and  fifth  columns  the  correspon^ng  ranges 
are  shown,  calculated  by  the  formulae  R  =  2(N^/Nf)x  (A  units,  10“°  cm)  and 
(pg/cni^)^  respectively. 


The  ranges,  calculated  by  the  formvila  2(N(,/Nf)x,  are  the  vector  sums  of 
the  tortuoiis  paths  traveled  by  the  recoiling  atoms.  Theoretical  ranges,  cal- 
cvilated  by  assuming  billiard-ball  collisions,  were  made  according  to  Seitz 
and  Koehler, 5  and,  in  general,  the  theoretical  ranges  were  many  times  higher 
than  those  given  in  Table  I.  Nielsen,12  following  the  lines  of  Bohr,^^  ob¬ 
tains  the  following  theoretical  expression  for  the  range; 


R  =  0.6  • 


2/3  ^  ,2/3 


1/2 


)  A  +  A^  2 

-  •  - —  •  A^  •  (pg/cm  )  , 


where  Z]^,  A^^  and  Z2,  A2  are  the  atomic  numbers  and  masses  for  the  incoming 
particles  and  the  target  atoms,  respectively,  and  Ei  is  the  energy  of  the 
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Tatle  I 


FRACTION  (Nc/Nf)  OF  ACTIVATED  AT0^6  RECOILING  OUT  OF  0.001-IN.  FOILS 
AND  CALCULATED  RANGES  OF  THE  RECOIL  ATOMS 
(23-Mev  BREMSSTRAHLUNG  INCIDENT  ON  FOILS) 


Recoiling 

Atom 

Most 

Probable 

Recoil 

Energy 

(kev) 

1 

/ 

Nc/Nf  X  10 

R,  Calc. 
Range  in 
%  (10“®  cm) 

R,  Calc . 
Range  in 

Hg/cm2^ 

Theoretical 
Range  in 
|ig/cm2c 

c^i 

130^ 

216  +20 

11,000®'^ 

117 

216® 

F^® 

85“^ 

182  ii8 

9,200® 

200 

118 

Cl^^ 

45*^ 

104  +10 

5,300® 

90 

31® 

Ti^5 

33^ 

6.1  io.8 

310 

14 

24 

Fe^^ 

30*^ 

11.0  +1.7 

560 

44 

19 

Zn^^ 

25^ 

2.7  I0.3 

137 

9.8 

15 

62 

Cu 

25  I16 

3.1  I0.2 

163 

14.5 

16 

Cu^"^ 

25  116 

3.2  +0.2 

157 

14.0 

16 

Mo^^ 

16*' 

1.4  I0.4 

71 

7.5 

8.1 

A  106 

Ag 

ll.^' 

1.4  io.3 

71 

7.5 

6.5 

A  196 

Au 

9  i6 

0.55  i0.06 

28 

5.4 

3.2 

a 


Errors  are  two  standard  deviations  of  counting  statistics  only. 


^Calculated  from  the  previous  column  by  the  formvila  Rp,  where  p  is  the 
density  of  the  target  foil. 

^Calculated  theoretical  expected  range  given  by  K.  0.  Nielsen  (Reference 
12,  p.  73,  Eq.  (13)),  where 


R  =  0.6 


(zf3  .  z|/3) 
^1^2 


„  1/2 


+  A„ 


•  AgE^  (pg/cm  )  , 


1  1 

a  formula  based  on  a  paper  by  N.  Bohr.  The  energy  Inserted  in  the  formvala 
is  the  average  recoil  energy. 


'Estimated  by  assuming  that  the  respective  evaporated-neutron  spectra 
were  similar  in  shape  to  those  observed  for  copper  and  lead.  The  average 
recoil  energy  is  assumed  to  be  about  1.5  times  the  most  probable  recoil 
energy. 


'"Range  for  C  is  in  polystyrene  (CH)j  for  Cl^  ,  in  saran  (CHCl);  and 
for  f18,  in  Teflon  (CF2).  Other  ranges  are  in  their  respective  metallic 
lattices . 


Table  I—continued 


f*  XI 

^Range  of  C  (for  same  recoil  energy)  in  polyethylene  (CH2)  and  in 
cellulose  acetate  (Ci.0%.300  7^  10,600  tl,100  A  and  11,900  21,800  A, 
respectively. 

®With  the  exception  of  the  hydrogen  atoms,  it  was  assumed  that  the 
target  atoms  shared  eqvially  in  stopping  the  recoiling  atoms.  Rydrogen 
atoms  were  neglected  in  the  calcvtlations. 


incoming  particle  in  kev.  In  the  application  of  Nielsen's  equation  to  the 
data  of  Table  I,  we  have  let  Z-^,  Aq^,  represent  the  recoiling  atom  and 
Z2>  A2  represent  the  target  atoms.  For  the  data  pertaining  to  the  metals, 
equals  Z2  and  A2  eqraals  Ai  +  1.  In  calculating  the  theoretical  ranges 
for  the  nonmetals  (cH,  f18,  and  Cl3^)  in  their  respective  plastic  matrices, 
Z2  and  A2  have  been  averaged  stoichiometrically  over  the  target-foil  atoms, 
excluding  the  hydrogen  atoms.  Theoretical  ranges  calculated  according  to 
Nielsen  are  given  in  the  last  column  of  Table  I.  Agreement  between  cal¬ 
culated-  and  theoretical-range  Values  is  rather  surprising  (within  a  factor 
of  two  for  most  ranges),  in  view  of  the  claim  by  Nielsen  that  the  theoretical- 
rat]ge  formvila  is  not  applicable  when  the  recoiling  atom  and  the  target  atoms 
are  common— i.e.,  having  the  same  Z,  A— because  the  recoi3  atoms  may  lose 
all  of  their  energy  in  a  head-on  collision  with  a  common  target  atom. 

Figure  1  shows  that  the  fraction  of  metal  atoms  that  recoil  out  of 
0.001-ln.  foils  varies  exponentially  with  l/A,  while  fltrarine,  chlorine,  and 
carbon  in  the  plastic  matrices  increase  slowly  with  l/A.  I'f'  it  is  assumed 
that  the  photoneutron  spectra  of  all  the  elements  are  centered  at  the  same 
energy  and  have  the  same  shape,  then  the  abscissa  of  Fig.  1  essentially 
represents  the  energy  of  the  recoil  atoms.  No  explanation  is  offered  for 
the  data  of  Fig.  1,  which  depend  on  parameters,  mass,  energy,  and  lattice 
structures.  By  performing  recoil  experiments  with  photoproton,  photodeuteron, 
and  photoalpha  transmutations,  it  will  be  possible  to  see  the  effect  on 
Nc/Nf.  by  varying  the  recoil  energy  only. 

It  should  be  noted  that  the  Nc/Nf  fractions  for  Cu^2  and  cu6^  agree  well 
within  experimental  error.  Since  the  photoneutron  thresholds  for  Cu°2  and 
Cu6i4^  are  nearly  equal  (IO.7  and  10.0  Mev,  respectively)  and  the  corresponding 
photoneutron  cross  sections  are  roughly  s\perimposable  after  normalization 
by  a  constant  factor,  the  kinetic  energies  of  the  evaporated  neutrons  are 
expected  to  be  identical.  By  measxiring  the  average  evaporated  photoneutron 
energy,  using  the  sinple  recoil  technique  outlined  in  this  paper,  it  is 
possible  to  compare  the  nuclear  tenperatures  of  pairs  of  excited  nuclei. 

The  conparlson  should  be  restricted  to  a  narrow  range  of  atomic  numbers  be¬ 
cause  of  the  variation  in  recoil  range  with  electronic  structure. 

Proton-energy  spectra  resulting  from  22-Mev  bremsstrahlixng  have  been 
measured  for  at  least  fourteen  photoproton  reactions  by  nuclear-emxilsion 
techniques.  In  general,  the  most  probable  proton  energy  is  about  6  Mev. 
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I/A  X  10"’ 

Fig.  1— Fraction  of  atoms,  N^/Nf,  recoiling  out  of  0.001-in.  foils 
as  a  function  of  Inverse  mass  of  recoiling  atom,  l/A 


The  most  probable  proton  energy  and  the  corresponding  energy  width  may  be 
converted  to  the  most  probable  recoil  energy  for  the  residual  atom  and  the 
corresponding  energy  width.  For  example,  recoil  energies  of  such  typical 
atoms  as  Lio,  Nb91,  an^  Hfl^  have  been  observed  to  be  about  200  i8o,^° 

60  +40,^'  and  40  il5  kev,^°  respectively.  About  805^  of  the  recoils  are 
within  the  respective  energy  limits.  Even  higher  recoil  energies  may  be 
expected  from  photodeuteron,  phototriton,  and  photoalpha  reactions  because 
the  emitted  particles  are  heavier  and,  in  the  case  of  photoalpha  emission, 
more  highly  charged.  By  choosing  the  proper  photoreaction,  measurements 
may  be  made  of  the  energies  of  other  recoil  atoms  in  the  range  10  to  100  kev 
in  lattices  composed  of  similar  materials  and  also  in  lattices  of  dissimilar 
materials,  such  as  alloys  and  crystals. 
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THE  EFFECTS  OF  NUCLEAR  IRRADIATION  ON 
METALLIC  AND  NONMETALLIC  MAGNETIC  MATERIALS 

by 

E.  I.  Salkovitz,  A  I.  Schindler,  G.  S.  Ansell 

U.  S.  Naval  Research  Laboratory 
Washington,  D.  C. 

Extensive  investigations  of  the  effects  of  nu¬ 
clear  environments  upon  magnetic  materials  have  been 
undertaken.  A  major  aim  of  the  program  is  to  obtain 
basic  information  concerning  the  mechanisms  produc¬ 
ing  the  observed  effects.  More  than  100  samples  have 
been  irradiated  in  the  Brookhaven  graphite  reactor 
at  an  integrated  flux  of  10^"^  NVT.  The  materials 
studied  have  been  in  the  form  of  toroids  or  rods  and 
have  consisted  mainly  of  various  ferrites  and  square 
loop  and  high  permeability  alloys.  In  addition, 
discs  of  permanent-magnet  type  ferrices  and  portions 
of  magnetic  devices  have  been  irradiated.  A  detailed 
discussion  will  be  given  of  the  method  of  canning 
temperature  control, and  the  means  by  which  pre-  and 
post-irradiation  magnetic  measurements  were  made. 

During  the  past  decade  there  have  appeared  many  papers 
on  irradiation  effects  in  many  materials.  Such  effects  on 
magnetic  materials,  however,  have  escaped  concerted  atten¬ 
tion.  A  notable  exception  has  been  the  work  of  Sery,  Fischell, 
and  Gordon'^)  dealing  primarily  with  several  alloys.  Conse¬ 
quently  at  the  U.  S.  Naval  Research  Laboratory  a  broad  pro¬ 
gram  is  underway  to  determine  the  changes  which  may  occur  in 
typical  magnetic  substances  due  to  nuclear  irradiation,  and 
to  investigate  mechanisms  which  may  account  for  these  changes. 
The  need  for  such  a  study  is  obvious.  From  the  theorist’s 
point  of  view,  new  Information  on  the  behavior  of  matter  will 
be  obtained.  And  from  the  practical  point  of  view,  there  is  a 
considerable  degree  of  urgency  for  this  type  of  data.  As  the 
application  of  nuclear  energy  increases  there  will  be  an 
accompanying  increase  in  the  use  of  magnetic  devices  operating 
in  radiation  fields.  For  example,  it  has  come  to  the  authors 
attention  that  an  accelerator  is  to  be  constructed  using  nearly 
a  ton  of  a  nickel  ferrite.  How  such  a  ferrite  will  behave  in 
a  billion  electron  volt  accelerator  at  the  operating  frequency 
is  an  unanswered  question  because  trhe  appropriate  investiga¬ 
tions  have  not  been  conducted  yet. 

At  the  outset  it  must  be  appreciated  that  an  investigation 
of  irradiation  effects  is  in  a  sense  an  investigation  of  environ 
mental  effects.  Consequently,  it  is  not  only  mandatory  to  know 
as  much  as  possible  about  the  material  under  study,  it  is 
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equally  necessary  to  be  able  tb  delineate  the  various  en¬ 
vironmental  parameters.  Unfortunately  this  latter  task  is  a 
very  formidable  one.  If  a  given  reactor  is  used  for  the  study 
it  is  desirable  to  state  the  power  level  (which  may  vary  dur¬ 
ing  the  test),  the  ambient  temperature  in  the  vicinity  of 
exposure,  the  energy  spectrum  of  radiation,  and  the  total  flux 
received  by  the  material  studied.  On  the  other  hand  the 
research  should  be  so  conducted  that  the  effects  may  be  dis¬ 
cussed  in  terms  of  the  influence  of  the  specific  variables 
already  mentioned  and  in  addition  to  these  such  parameters  as 
dose  rate  and  recovery  processes  which  may  occur  during  and 
after  irradiation.  It  should  be  clear  therefore  that  the  pro¬ 
gram  discussed  here  is  to  be  considered  still  in  the  exploratory 
stage.  Sufficient  results  have  been  uncovered  which  are  worth 
reporting,  but  which  indicate  that  further  study  is  required. 

To  complement  an  Investigation  of  the  behavior  of  irradia¬ 
ted  ferrites  in  the  microwave  region,  it  is  desirable  to 
examine  the  behavior  in  the  D.C..  and  60  cycle  region.  For  the 
sake  of  expediency  a  variety  of  commercially  available  ferrites 
have  been  studied.  Very  serious  disadvantages  exist  with  this 
approach.  Usually  the  exact  chemical  composition  and  manu¬ 
facturing  procedures  are  complicated  and  trade  secrets.  Study¬ 
ing  commercial  ferrites  however  furnishes  immediate  information 
on  materials  of  practical  importance.  Suffice  to  say  that 
experiments  on  more  controllable  ferrites  are  underway  and  will 
be  reported  in  the  near  future. 

Four  separate  sets  of  experiments  were  conducted  using  the 
Brookhaven  reactor  operating  between  15  and  22  megawatts.  Dur¬ 
ing  each  run,  ferrites  and  magnetic  alloys  (which  were  also 
studied)  were  distributed  in  aluminum  cans  over  a  3  or  4  foot 
portion  of  a  channel  in  the  reactor.  Consequently  the  total 
integrated  flux  (fast)  received  by  the  different  samples  varied 
between  the  limit  of  6  x  IQI®  NVT  and  3  x  lO^^  NVT.  By  means 
of  forced  air  cooling  the  temperatures  of  all  samples  were 
kept  below  75°C  and  some  as  low  as  30° C.  Approximately  100 
ferrite  samples  of  about  20  different  compositions  were  studied, 
as  well  as  35  metal  alloy  samples  of  9  different  compositions. 

To  obtain  60  cycle  hysteresis  loops,  toroids  were  wound 
.with  appropriate  primary  and  secondary  colls.  Groups  of  4  to 
7  such  toroids  were  then  moxmted  on  a  sample  holder.  Figure  1 
shows  such  a  sample  holder  patterned  after  that  of  the  early 
work  of  Sery,  Flschell  and  Gordon.  As  can  be  seen,  leads  from 
the  coils  were  connected  to  ceramic  octal  receptacles.  Electri¬ 
cal  con^ctlons  could  then  be  made  with  a  hysteresis  loop- 
tracer'^'  and  the  loops  presented  on  an  X-Y  recorder.  These 
measurements  were  made  remotely^  both  before  and  after  irradia¬ 
tion,  using  the  same  equipment  and  all  measurements  were  made 
in  the  hot  cell  facilities  of  the  U.  S.  Naval  Research  Labora¬ 
tory.  When  the  "before”  measurements  were  completed,  the 


2 


toroids  mounted  on  their  sample  holders  were  placed  in 
aluminum  Cans,  shipped  to  the  Brookhaven  National  Laboratory 
and  lnse]*ted  in  the  reactor  with  cdils  intact.  After  expo¬ 
sure  the  units  were  returned,  their  radioactivity  measured 
and  hysteresis  loops  again  obtained.  Figures  2  through  8 
are  representative  of  the  effects  of  the  irradiation  upon 
several  different  types  of  ferrites.  The  hysteresis  loops 
shown  are  direct  reproductions  of  the  loops  presented  on  the 
X-Y  recorder.  For  purposes  of  comparison  we  note  in  Table  I, 
changes  in  the  remanence  Br,  the  coercive  force  He  and  the 
ratio  Bn/Br  where  B^  is  the  maximum  flux  density.  If  a 
10  percent  change  in  any  of  the^e  parameters  arbitrarily  is 
considered  to  be  insignificant,  then  Samples  IV  and  II  which 
are  nickel,  zinc,  manganese  ferrites,  and  Sample  VI  which  is 
a  complex  ferrite  containing  mixtures  of  oxides  of  nickel , 
zinc-manganese,  magnesium,  and  copper  may  be  assumed  to  have 
been  unaffected.  All  other  ferrites  showed  larger  changes » 

The  kinks  which  appeared  in  the  hysteresis  loops  of  Sample 
VIII  (Mn,  Zn  ferrite)  and  Sample  IX  (Ni,  Zn,  Cu  ferrite)  after 
irradiation,  are  reproducible.  Such  kinks  appeared  in  more 
pronounced  fashion  in  several  irradiated  alloys.  Although  the 
work  on  alloys  will  be  covered  in  detail  in  a  paper  in  prepara¬ 
tion  it  is  germalne  to  indicate  a  similarity  between  ^e  "kink" 
phenomena  in  irradiated  ferrites  and  in  alloys.  Skich 

"kinks"  in  hysteresis  loops  have  been  found  in  unirradiated 
permalloys  and  have  been  studied  elsewhere.  These  kinks  may 
be  associated  with  the  order-disorder  processes  which  occur 
during  the  thermal  history  of  the  material.  Figure  9  shows 
loops  for  a  4-79  Mo  Permalloy  before  and  after  irradiation. 

The  kinks  in  the  "after"  curve  are  q^uite  apparent.  In  Figure 
10  are  the  loops  for  a  sample  of  Mumetal.  Figure  10a  shows 
the  "before"  curve.  Figure  10b  and  Figure  10c  show  "after" 
irradiation  loops  obtained  for  two  different  exposures.  Since 
these  were  exploratory  experiments,  the  radiation  flux  was  not 
carefully  monitored,  so  that  it  was  quite  possible  that  the 
radiation  environment  was  very  different  for  the  two  runs. 

The  Important  feature  is  that  the  two  loops  are  different 
(Figures  10b  and  10c)  but  are  the  types  of  loops  which  are 
found  when  the  ordering  treatment  is  interrupted  and  hysteresis 
loops  obtained  for  permalloy  type  of  alloys.  It  is  conceivable 
that  for  the  ferrites  a  similar  explanation  may  be  offered. 
However,  this  phenomena  is  under  investigation  and  should  not 
be  discussed  until  further  data  is  available.  Additional  re¬ 
sults  on  a  variety  of  commercial  alloys  are  shown  in  the 
succeeding  figures. 

As  an  exploratory  experiment  in  which  the  ratio  of  thermal, 
to  epl-thermal,  to  fast  neutrons  was  varied,  two  samples  of  a 
particular  ferrite  were  placed  in  close  proximity  in  the 
reactor.  One  of  these  however  was  wrapped  in  cadmium  foil. 
Table  II  indicates  that  the  energy  spectnm  of  the  neutrons 


3 


does  influence  the  effects  observed.  Obviously  more  work  of 
this  nature  is  desired. 

The  effect  of  irradiation  upon  the  Curie  temperature  of 
all  the  ferrites  were  investigated  by  Mr.  Glenn  C.  Bailey  of 
this  Laboratory.  For  this  experiment  the  samples  received  a 
total  flux  of  1.24  X  IQI^  nvt  (fast).  No  change  in  the  Curie 
temperature  was  found  for  any  of  the  ferrites. 

The  preliminary  results  of  still  another  experiment  may 
be  of  Interest.  The  effect  of  irradiation  (1.24  x  lO^'  NVT 
fast)  upon  the  magnetic  moment  of  d  Fe304  and  Fe203  has  been 
obtained  by  Dr.  W.  E,  Henry  and  one  of  the  authors  (EIS) . 
Whereas  the  moment  was  found  to  decrease  by  25  percent  in  the 
case  of  a  Fe304  a  decrease  of  only  2  percent  was  found  for 
Y  Fe203.  This  strongly  suggests,  as  would  be  expected,  that 
the  radiation  produces  a  redistribution  of  the  cations.  Neu¬ 
tron  diffraction  experiments  are  underway  to  determine  if  this 
is  the  case. 

From  the  foregoing  it  is  apparent  that  certain  magnetic 
properties  of  the  materials  investigated  may  be  altered  by 
neutron  irradiation.  These  are  primarily  the  structure  sensi¬ 
tive  properties.  The  degree  of  the  effect  varies  with  the 
chemical  composition  and  past  history  of  the  material.  Both 
factors  of  course  determine  the  cation  distribution.  Suf¬ 
ficient  variables  however  have  been  uncovered  that  until 
further  controlled  experiments  are  conducted  it  would  be  un¬ 
wise  to  dogmatically  advance  mechanisms  for  the  observed 
effects.  Having  uncovered  these  effects  in  the  very  low  fre¬ 
quency  region  it  is  now  desirable  to  investigate  the  behavior 
of  Irradiated  ferrites  at  higher  frequencies.  'Hils  will  be 
discussed  in  a  forthcoming  paper. 
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Sample  holder  used  for  Irradiation  experiments. 
Hysteresis  loops  before  and  after  Irradiation  for 
a  commercial  nickel  ferrite. 

B^prsteresls  loops  before  and  after  Irradiation  for 
w  commercial  nickel,  zinc,  manganese  ferrite. 
Hysteresis  loops  before  and  after  Irradiation  for 
another  commercial  nickel^  zinc,  manganese  ferrite. 
Hysteresis  loops  before  and  after  Irradiation  for 
a  commercial  nickel,  zinc,  manganese,  magnesium 
ferrite. 

Hysteresis  loops  before  and  after  Irradlatlop  for 
a  commercial  nickel,  zinc,  copper  ferrite. 
Hysteresis  loops  before  and  after  Irradiation  for 
another  commercial  nickel,  zinc,  copper  ferrite. 
Hysteresis  loops  before  and  after  Irradiation  for 
a  commercial  manganese,  zinc  ferrite. 

Hysteresis  loops  before  and  after  Irradiation  for 
a  commercial  4-79  Permalloy. 

Hysteresis  loops  before  and  after  Irradiation  for 
a  commercial  Mumetal,  (b)  and  (c)  represents  two 
different  samples  Irradiated  at  different  times. 

-Hysteresis  loops  before  and  after  irradiation  for 
commercial  alloys  as  indicated. 
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THE  EFFECT  OF  NUCLEAR  RADIATION  ON 
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Two  AN/ARC-34  UHF  Comnwand  Communications  Systems  were 
exposed  to  high  energy  neutron  and  gamma  radiation  from  the 
Ground  Test  Reactor.  When  set  number  C-3735  failed  due  to 
radiation  Induced  damage  to  a  1N69  crystal  diode  it  had  received 
an  Integrated  epicadmium  neutron  flux  of  4.68  x  10^3  neutrons/cm^ 
and  an  integrated  gamma  dose  of  4.71  x  10*^  Roentgen.  When 
set  number  C-3731  failed  due  to  radiation  Induced  damage  to  a 
1N69  crystal  diode  It  hod  received  an  Integrated  epicadmium 
neutron  flux  of  4.22  x  10^3  neutrons/cm^  and  an  integrated 
gamma  dose  of  7.03  x  10^  Roentgen.  A  second  failure  occurred 
In  set  number  C-3731  at  an  Integrated  epicadmium  neutron  flux 
of  2.82  X  10^^  neutrons/cm^  and  an  integrated  gamma  dose  of 
4.7  X  10^  Roentgen  when  a  0.22  microfarad  paper  capacitor  failed. 
No  transient  responses  were  noted  In  system  performance  while  sub¬ 
jected  to  radiation. 
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A  number  of  organic  engineering  materials  and  components  used 
by  the  General  Electric  Company  on  J-79  turbojet  engines  have 
been  and  are  presently  being  irradiation  tested  at  the  AN  PD  in  on 
effort  to  determine  the  applicability  of  these  materials  in  propulsion 
machinery  for  nuclear  powered  flight.  The  materials  and  components 
are  irradiated  in  anticipated  thermal  environments  while  being  sub¬ 
jected  to  functional  pressure  and  fluid  flow.  The  materials  and 
components  under  investigation  Include  a  turbojet  lubricant,  hydrau¬ 
lic  fluid,  fuel,  elastomer  seals,  gaskets,  and  hoses.  The  elastomer 
materials  consist  of  seven  commercial  products  in  the  form  of  fabri¬ 
cated  O-rIngs,  gaskets  and  flexible  hoses. 

This  paper  will  present  the  past,  present  and  planned  materials  and 
components  irradiation  programs,  post-irradiation  materials  evalu¬ 
ation,  the  necessity  for  information  of  this  kind,  and  other  problems 
and  considerations  associated  with  applications  proving  radiation 
tolerance  of  organic  engineering  materials. 
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